B. Runge-Kutta 法で解く

今さら数値計算する気がなくなる…まあ、気が向いたら。

と書いたのですが、 歳は取りたくないもので、 ずっと昔、授業 (『振り子の運動』) でやったことがありました。 その計算でも、 確かに振幅 $1$ ラジアンで $6.63$ % 周期が長くなっていました。 問題として、完全楕円積分を級数展開で計算しろ、なんて書いてある (笑)。 ともあれ、全然違う方法で計算して得た数値が良く一致したので、とりあえず安心。