a = sin²(Δφ/2) + cos(φ1) · cos(φ2) · sin²(Δλ/2)
c = 2 · atan2(√a, √(1−a))
d = R · c
R = 6371 km(地球の平均半径)、φは緯度、λは経度、Δは差分。この公式は地球を完全な球体と仮定しており、誤差は約0.3%です。