ドルマン・プリンス（Excel VBA移植版）でアレンストーフ軌道を計算。

Dormand & Prince (ルンゲ・クッタ8次）　DOP853　VBA移植版で制限三体問題の例、アレンストーフ軌道を計算してみる。式は

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次のドルマン・プリンス法で計算してみよう。

一発で

が計算できた。

（オリジナルのDOP853は　http://www.unige.ch/~hairer/software.html　です）

ソースコードはこちら。

「DOP853.bas」をダウンロード

「Aren.bas」をダウンロード

ライセンスは

「licence.txt」をダウンロード

に従います。