こんにちは。 楕円弧長は、下記の第二種不完全楕円積分で表されます。与えられた緯度に対し赤道からの子午線弧長です。 プログラムでは小さい扁平率を仮定し、倍精度浮動小数点数に現れない$n^7$(第三扁平率の冪乗)以上の級数展開項は近似で切り捨てています。このベッセルの級数式は優れた収束ですね。これの級数係数の計算については分数の表示(リスト)。参考:https://en.wikipedia.org/wiki/Meridian_arc#Expansions_in_the_third_flattening_.28n.29 \textrm{An elliptic arc length formula from a modification of Bessel (1825)}\\ \begin{aligned} m(\varphi) &= \int_0^\varphi M(\varphi) d\var
![楕円弧長計算 - Qiita](https://cdn-ak-scissors.b.st-hatena.com/image/square/2cb7930a0befd7b58832da1a9d73389b6ebe4192/height=288;version=1;width=512/https%3A%2F%2Fqiita-user-contents.imgix.net%2Fhttps%253A%252F%252Fcdn.qiita.com%252Fassets%252Fpublic%252Farticle-ogp-background-9f5428127621718a910c8b63951390ad.png%3Fixlib%3Drb-4.0.0%26w%3D1200%26mark64%3DaHR0cHM6Ly9xaWl0YS11c2VyLWNvbnRlbnRzLmltZ2l4Lm5ldC9-dGV4dD9peGxpYj1yYi00LjAuMCZ3PTkxNiZoPTMzNiZ0eHQ9JUU2JUE1JTk1JUU1JTg2JTg2JUU1JUJDJUE3JUU5JTk1JUI3JUU4JUE4JTg4JUU3JUFFJTk3JnR4dC1jb2xvcj0lMjMyMTIxMjEmdHh0LWZvbnQ9SGlyYWdpbm8lMjBTYW5zJTIwVzYmdHh0LXNpemU9NTYmdHh0LWNsaXA9ZWxsaXBzaXMmdHh0LWFsaWduPWxlZnQlMkN0b3Amcz05ZWIxZjQzMWRmZGJlMmU4YWExYmUwODcxYmNhNDUzNw%26mark-x%3D142%26mark-y%3D112%26blend64%3DaHR0cHM6Ly9xaWl0YS11c2VyLWNvbnRlbnRzLmltZ2l4Lm5ldC9-dGV4dD9peGxpYj1yYi00LjAuMCZ3PTYxNiZ0eHQ9JTQwa2tkZCZ0eHQtY29sb3I9JTIzMjEyMTIxJnR4dC1mb250PUhpcmFnaW5vJTIwU2FucyUyMFc2JnR4dC1zaXplPTM2JnR4dC1hbGlnbj1sZWZ0JTJDdG9wJnM9YzRlZDdiMGE5MjhiZjU5MjBjNTZjZTliYzAwMmIwYWE%26blend-x%3D142%26blend-y%3D491%26blend-mode%3Dnormal%26s%3Ddf61fc4a9f0c307b2aaed9cffd7dc827)