a = sin²(Δφ/2) + cos(φ1) · cos(φ2) · sin²(Δλ/2)
c = 2 · atan2(√a, √(1−a))
d = R · c
其中 R = 6371 km（地球平均半径），φ 为纬度，λ 为经度，Δ 为差值。该公式假设地球为完美球体，误差约 0.3%。