逆二乗の中心力のもとで物体がどのように運動するかを調べるケプラー問題。微分方程式を解く能力がなければ楽しめないかと言えば、そういうわけでもない（物理を勉強する人にはもちろん必須）。
最近の強力なプログラミング環境を使えば、物体の軌道を簡単に視覚化できる。インタラクティブなものにするのも簡単（さすがに、微分方程式は立てられないといけない）。
Mathematicaならこんな感じ（インタラクティブに操作できるのは、重力の強さを決めるkと時間の上限tmax、初速度v0、初期座標pos）。
実際に動かしてみたい人は、Mathematica Player用のノートブックを使って欲しい。Mathematica Playerはフリーソフトウェアではないが無料で使える。ブラウザ上で動くようになるとすごいと思う。
Mathematicaの微分方程式ソルバーNDSolveはかなり強力で、このエントリのような単純なものだけではなく、摂動があるような問題も扱える。
解析的に解きたい人は、ランダウ＝リフシッツ『力学・場の理論』のような教科書を参照して欲しい（福島登志夫『天体の位置と運動』の方がやさしいかな ）。ジンマーマン, オルネス『物理学のためのMathematica』では、ハミルトン＝ヤコビ方程式を使う方法が紹介されている。
関連：プリンキピアを読む
Related posts:
