def combination_2_n(l): n, r = len(l), [] for i in range(0, n, 1): s = i + 1 for ii in range(s, n, 1): r.append([l[i], l[ii]]) return r # l1, l2 = [23, 123], [24, 124] def rad(d): return d * np.pi / 180.0 def compute_diff(l1, l2): lat1, lng1 = l1 lat…