アレンストーフ軌道を14次ルンゲクッタ法で計算（PARI/GP)

PARI/GPで35段14次ルンゲクッタを計算するシリーズ。アレンストーフ軌道も高精度に計算できるかな？と思ってやってみた。

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

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

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

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

μ=0.012277471

μ' = 1-μ

y1(0) = 0.994, y1'(0)=0,

y2(0)=0,y2'(0)=-2.00158510637908252240537862224

計算結果はこちら。

ソースはこちら。↓

「aren.gp」をダウンロード