cython
12 строк · 287.0 Байт
1import math2
3def great_circle(lon1, lat1, lon2, lat2):4radius = 3956 # miles5x = math.pi/180.06
7a = (90.0 - lat1)*x8b = (90.0 - lat2)*x9theta = (lon2 - lon1)*x10c = math.acos(math.cos(a)*math.cos(b) + math.sin(a)*math.sin(b)*math.cos(theta))11
12return radius*c13