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

だいぶ前に、

なんと三体問題の新しい解が13個も見つかっていた！さっそく計算して図示＆GIFアニメ化

という記事を書いた。

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

The 1223 new periodic orbits of planar three-body problem with unequal mass and zero angular momentum

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

とりあえず検証してみよう。ルンゲクッタ8次のDormand&Princeルーチン、DOP853を使っていると論文には書いてあるが、私もちょうどそのフォートランルーチンをExcelのVBAに移植しているのだ（この記事の後にリンクつけます）

では計算しよう。どれがいいかな、、、1223個もあるし、、、

とりあえずなんとなくみかんのあみあみ袋みたいなI.A68(0.5)というのにしよう。

論文の結果と同じになった。GIFアニメにもしてみよう。

1223個の新しい三体問題の周期解が見つかったということで、私もExcel VBAにDOP853移植ルーチンを使って計算、GIFアニメにしてみた。今回はI.A68(0.5)というもの。えぐい動きをしてるな、、、計算の詳細は続くツイートにて。 pic.twitter.com/BjCu7uJFIQ

— tomo (@tonagai) 2017年9月23日

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

次はIID2(2)ってやつ。

1223個の新しい三体問題の周期解が見つかったということで、私もExcel VBAへのDOP853移植ルーチンを使って計算、GIFアニメにしてみた。今回はII.D2(2)というもの。 pic.twitter.com/npOPHjfBSe

— tomo (@tonagai) 2017年9月23日

いろいろ計算しようかと思ったけど1223個か、、、萎える、、、

---

関連リンク：

Dormand-Prince(ルンゲ・クッタ8次）の有名なルーチンDOP853をExcel VBAに移植

Dormand-Prince(ルンゲ・クッタ8次）DOP853をExcel VBAに移植(説明編)