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

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

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

コメント

コメントを書く

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

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

最近の記事

最近のコメント

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