« マツダの水木一郎アニキが出てるCMのBGMはスペクトルマン? | トップページ | スターバックスダブルショットのCMのsmoke on the waterはSenor Coconut »

2009年11月22日 (日)

「昴問題」をsymplectic8次、Excelで計算してみる(VBAソース付き)。

久々にシンプレクティック法のプログラムをいじってみた。「常微分方程式の数値解法I」に載っていた7体問題で「昴問題」というらしい。

pxi '= ΣGmimj(xj - xi)/rij3

pyi '= ΣGmimj(yj - yi)/rij3

で、

(m1, ... , m7)=(1,2,3,4,5,6,7)

初期値は

(x1, ..., x7) = (3,3,-1,-3,2,-2,2)

(y1, ..., y7) = (3,-3,2,0,0,-4,4)

(x1', ..., x7') = (0,0,0,0,0,1.75,-1.5)

(y1', ..., y7') = (0,0,0, -1.25,1,0,0)

でt=3まで計算。本では自動ステップ幅調整が必要な例として書いてあるが(衝突があるから)、無理やりシンプレクティックで最初からステップ幅めちゃくちゃ小さくして計算。でも1分くらいですよ。計算は。

で結果は本の結果と同等なものが得られた。

Subaru

色分けしたほうがわかりやすいな。

ソースはこちら。

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

« マツダの水木一郎アニキが出てるCMのBGMはスペクトルマン? | トップページ | スターバックスダブルショットのCMのsmoke on the waterはSenor Coconut »

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

コメント

コメントを書く

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

トラックバック

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

この記事へのトラックバック一覧です: 「昴問題」をsymplectic8次、Excelで計算してみる(VBAソース付き)。:

« マツダの水木一郎アニキが出てるCMのBGMはスペクトルマン? | トップページ | スターバックスダブルショットのCMのsmoke on the waterはSenor Coconut »

最近のコメント

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