« #せやねん メチャ売れのコーナー、9/23の特集はエースコックの新食感焼そば モッチッチ。 | トップページ | 味仙監修の名古屋名物台湾ラーメンのカップ麺を食べた。なかなかの辛さでした。 »

2017年9月23日 (土)

1223個の新しい三体問題の周期解が見つかったということで、私もExcel VBAにルンゲクッタ8次のDOP853を移植したもので計算してみた。なかなか面白い動きをする。

だいぶ前に、

という記事を書いた。
このときも感心していたのだが、先日こんな論文があるのを知った。

なんじゃこれ!力技で1223個も見つけたのか!

とりあえず検証してみよう。ルンゲクッタ8次のDormand&Princeルーチン、DOP853を使っていると論文には書いてあるが、私もちょうどそのフォートランルーチンをExcelのVBAに移植しているのだ(この記事の後にリンクつけます)
では計算しよう。どれがいいかな、、、1223個もあるし、、、
とりあえずなんとなくみかんのあみあみ袋みたいなI.A68(0.5)というのにしよう。
3body_ia68_05
論文の結果と同じになった。GIFアニメにもしてみよう。

うーん、なんちゅう動きをしてるんだ(ちょっと表示の刻みが粗すぎるか、、、)

次はIID2(2)ってやつ。

3body_iid2_2

いろいろ計算しようかと思ったけど1223個か、、、萎える、、、
---
関連リンク:

« #せやねん メチャ売れのコーナー、9/23の特集はエースコックの新食感焼そば モッチッチ。 | トップページ | 味仙監修の名古屋名物台湾ラーメンのカップ麺を食べた。なかなかの辛さでした。 »

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

日記・コラム・つぶやき」カテゴリの記事

コメント

コメントを書く

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

トラックバック

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

この記事へのトラックバック一覧です: 1223個の新しい三体問題の周期解が見つかったということで、私もExcel VBAにルンゲクッタ8次のDOP853を移植したもので計算してみた。なかなか面白い動きをする。:

« #せやねん メチャ売れのコーナー、9/23の特集はエースコックの新食感焼そば モッチッチ。 | トップページ | 味仙監修の名古屋名物台湾ラーメンのカップ麺を食べた。なかなかの辛さでした。 »

最近のコメント

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