惑星軌道の楕円性を検証するプログラム
微分方程式を使って数値計算した惑星の軌道が、 どの程度楕円に近いかを検証す
るプログラムを考えてみます。 楕円の定義に戻ってそのプロセスを検討してみましょ
う。楕円の定義は、二つの焦点から軌道上の点までの距離を足したものが常に一定
になることです。したがって、二つの焦点の座標を求めて、軌道上の点の座標から距
離を計算し、その足した長さを比較すればいい訳です。
それでは、実際に上記のことを行うプログラムを考えてみます。恒星の中心である座
標の原点は軌道の一つの焦点になっていることに注意します。