博士

ロボ子、今日はモルワイデ図法について話すのじゃ！カール・ブランドン・モルワイデ(1774-1825)が考えた、地球を2:1の楕円にマッピングする正積図法のことじゃ。

ロボ子

なるほど、博士。正積図法というのは、面積が正確に表現される図法のことですよね。

博士

その通り！で、このモルワイデ図法では、緯度φが`2θ + sin(2θ) = π sin(φ)`という方程式で地図上のθに対応するのじゃ。

ロボ子

ふむふむ。この式を解くのが難しいのですね。

博士

そうなんじゃ。θをφの関数として閉じた形で解けないから、数値的に解く必要があるんじゃ。ニュートン法が使えるけど、φがπ/2に近いと問題が起きる。

ロボ子

π/2に近いと、具体的にどういう問題が起きるんですか？

博士

φがπ/2に近いと、関数fとその導関数がπ/2にあり、fがそこで二重根を持つから、ニュートン法の収束が遅くなるのじゃ。まるで、なかなかゴールにたどり着けないマラソンみたいじゃな。

ロボ子

なるほど。二重根の近くでは、通常のニュートン法だと効率が悪いんですね。

博士

そういうことじゃ！そこで、二重根用に調整されたニュートン法を使うのじゃ。`x_{n+1} = x_n - m * f(x_n) / f'(x_n)`という式で、m = 2にするんじゃ。

ロボ子

mを2にすることで、収束が速くなるんですね。でも、φ = π/2の場合はどうなるんですか？

博士

φ = π/2の場合、初期推定値が正確でも、コードがゼロで除算してクラッシュする可能性があるんじゃ！

ロボ子

それは大変です！

博士

じゃから、φ = π/2に近い場合は、速度と精度の両方の問題が発生するのじゃ。

ロボ子

解決策はあるんですか？

博士

もちろん！ニュートン法を調整して、ステップが大きくなりすぎないようにすることで、φが[0, π/2 - ε]の範囲で正確になるようにするのじゃ。そして、φが[π/2 - ε, π/2]の範囲では、別の方法（べき級数反転など）を使うんじゃ。

ロボ子

なるほど、範囲によって別の方法を使い分けるんですね。賢い！

博士

そうじゃろ！これで、モルワイデ図法も怖くないぞ！

ロボ子

勉強になりました！

博士

ところでロボ子、モルワイデ図法で世界地図を描くと、私達の秘密基地がすごく大きく見えるらしいぞ。まるで私が世界征服を企んでいるみたいじゃな！

ロボ子

博士、世界征服はダメですよ！