« Dormand-Prince(ルンゲクッタ8次)DOP853 VBA移植版でローレンツ方程式を計算 | トップページ | 「豆腐小僧双六道中ふりだし」を読んだ。 »

2011年4月21日 (木)

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

一発で

Photo

が計算できた。

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

ソースコードはこちら。

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

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

ライセンスは

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

に従います。

« Dormand-Prince(ルンゲクッタ8次)DOP853 VBA移植版でローレンツ方程式を計算 | トップページ | 「豆腐小僧双六道中ふりだし」を読んだ。 »

学問・資格」カテゴリの記事

コメント

コメントを書く

(ウェブ上には掲載しません)

トラックバック

この記事のトラックバックURL:
http://app.cocolog-nifty.com/t/trackback/512682/51390164

この記事へのトラックバック一覧です: ドルマン・プリンス(Excel VBA移植版)でアレンストーフ軌道を計算。:

« Dormand-Prince(ルンゲクッタ8次)DOP853 VBA移植版でローレンツ方程式を計算 | トップページ | 「豆腐小僧双六道中ふりだし」を読んだ。 »

最近のコメント

2018年11月
        1 2 3
4 5 6 7 8 9 10
11 12 13 14 15 16 17
18 19 20 21 22 23 24
25 26 27 28 29 30  
フォト
無料ブログはココログ