アレンストーフ軌道をGeoGebra4.2でGIFアニメにしてみる。

昨日のピタゴラスの三体問題に続いて、次は制限三体問題の例、アレンストーフ軌道を計算してみる。式は

y1" = y1 + 2 y2' - μ' (y1 + μ）/D1　- μ(y1 - μ') / D2

y2" = y2 - 2 y1' - μ' y2 /D1　- μ y2 - μ' / D2

ここで

D1 = ((y1 + μ)^2 + y2^2) ^ (3/2)

D2 = ((y1 - μ')^2 + y2^2) ^ (3/2)

また

μ= 0.012277471

μ' = 1-μ

です。

計算はルンゲクッタ8次のDormand&Price法のDOP853のExcel VBA移植版を使用して、結果をGeoGebraにコピーしてます。