from math import radians, cos, sin, asin, sqrt
def distance(lat1, lat2, lon1, lon2):
# math contains a unit
# radians on the . function
# which converts angles from the degree system to the radius system
lon1 = radians(lon1)
lon2 = radians(lon2)
lat1 = radians(lat1)
lat2 = radians(lat2)
# Haversine formula
dlon = lon2 - lon1
dlat = lat2 - lat1
a = sin(dlat / 2)**2 + cos(lat1) * cos(lat2) * sin(dlon / 2)**2
c = 2 * asin(sqrt(a))
# Earth's radius in kilometres.
r = 6371
# Earth's radius in kilometres.
return(c * r)
# Test
lat1 = 53.32055555555556
lat2 = 53.31861111111111
lon1 = -1.7297222222222221
lon2 = -1.6997222222222223
print(distance(lat1, lat2, lon1, lon2), "K.M")