地球を真球と仮定して球面余弦定理に当てはめれば近似値が出る。 ググるとhaversine fomula（半正矢関数）という関数を使った計算式が出てくる。これは三角関数表を引く回数を減らすために球面余弦定理を変形したもの。PCではsinやcosを簡単に計算できるので…