« 「バガボンド33巻」を読んだ。「あきらめろ」というセリフに衝撃! | トップページ | van der Pol方程式を35段14次ルンゲクッタ法で計算(PARI/GP) »

2010年6月 2日 (水)

アレンストーフ軌道を35段14次ルンゲクッタ法(適応刻み幅機能付き)で計算(PARI/GP)

高次ルンゲクッタを導入したのはいいのだが、何も考えず刻み幅を決めるととんでもなく時間がかかってしまう。もともとの公式は埋め込み型で、刻み幅が調整できるのにサボっていてその機能使っていなかった。そこでようやく今頃になってその機能入れてみた。解いたのは、

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

計算結果はこちら。同じ計算結果を得るのに前は9000ポイント必要だったけど、1800ポイントまで減少した。図は両方を重ねているけど、全く重なっている。

Arendiff

ソースはこちら。↓

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

« 「バガボンド33巻」を読んだ。「あきらめろ」というセリフに衝撃! | トップページ | van der Pol方程式を35段14次ルンゲクッタ法で計算(PARI/GP) »

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

コメント

コメントを書く

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

トラックバック

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

この記事へのトラックバック一覧です: アレンストーフ軌道を35段14次ルンゲクッタ法(適応刻み幅機能付き)で計算(PARI/GP):

« 「バガボンド33巻」を読んだ。「あきらめろ」というセリフに衝撃! | トップページ | van der Pol方程式を35段14次ルンゲクッタ法で計算(PARI/GP) »

最近のコメント

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  
フォト
無料ブログはココログ