« NHKスペシャル "いのちの記録"を未来へ~震災ビッグデータ~ をさっきから見てました。 | トップページ | LambertのW関数をkeisan.casio.jpの自作式として作ってみた。 »

2013年3月 4日 (月)

アレンストーフ軌道をGeoGebra4.2でGIFアニメにしてみる。

昨日のピタゴラスの三体問題に続いて、次は制限三体問題の例、アレンストーフ軌道を計算してみる。式は

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-μ

です。

Aren


計算はルンゲクッタ8次のDormand&Price法のDOP853のExcel VBA移植版を使用して、結果をGeoGebraにコピーしてます。

« NHKスペシャル "いのちの記録"を未来へ~震災ビッグデータ~ をさっきから見てました。 | トップページ | LambertのW関数をkeisan.casio.jpの自作式として作ってみた。 »

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

コメント

コメントを書く

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

トラックバック

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

この記事へのトラックバック一覧です: アレンストーフ軌道をGeoGebra4.2でGIFアニメにしてみる。:

« NHKスペシャル "いのちの記録"を未来へ~震災ビッグデータ~ をさっきから見てました。 | トップページ | LambertのW関数をkeisan.casio.jpの自作式として作ってみた。 »

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