1223個の新しい三体問題の周期解が見つかったということで、私もExcel VBAにルンゲクッタ8次のDOP853を移植したもので計算してみた。なかなか面白い動きをする。
だいぶ前に、
という記事を書いた。
このときも感心していたのだが、先日こんな論文があるのを知った。
なんじゃこれ！力技で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個か、、、萎える、、、
---
関連リンク：
« #せやねん メチャ売れのコーナー、9/23の特集はエースコックの新食感焼そば モッチッチ。 | トップページ | 味仙監修の名古屋名物台湾ラーメンのカップ麺を食べた。なかなかの辛さでした。 »
「学問・資格」カテゴリの記事
- Wordleを100回やってみて正解率の分布を見た。1回目で当たる確率をpとして、2回目からはk倍ずつ確率が上がっていく(kp,k²p,k³p…）とすると、私の結果はp=0.0089, k=2.4になった。(2022.05.27)
- 高周波回路シミュレータQucsStudioを使ってみる（その５）マルチポート（４ポートだとs4pみたいな）のSパラメータデータをTouchstoneファイルに出力しようとするとExport to SnPでは2ポートになってしまう。Project→Import Dataからポート数を手で入れないと。(2022.05.24)
- 高周波（RF・マイクロ波・ミリ波・5G）関連ニュース2022年5月19日 IEEE Microwave Magazineで磁性体を使わない非可逆デバイス、Microwave Journalは5Gスモールセル、THz、Movanoがミリ波血圧＆血糖値計発表、三菱の3Dプリンタ衛星アンテナ、など。(2022.05.19)
- シン・ウルトラマンをIMAXで観てきた。面白かった！物理学者の有岡君がラグランジアン,AdS,プランクブレーンとかしゃべってる！Randall–Sundrumも。マグカップは標準理論＋重力で、物理監修は予想通り橋本幸士さんでした！ブラックホールの構造もタイムリー！(2022.05.14)
- 個人の正解確率pとして多数決をとったときの正解確率が、pが0.5に近づくとどのくらいの人数で0.9を超えるか計算してみた。p=0.50005だと日本の人口くらいの多数決。（神とさざなみの密室読んで正当性確率が気になったので）(2022.05.10)
「日記・コラム・つぶやき」カテゴリの記事
- Wordleを100回やってみて正解率の分布を見た。1回目で当たる確率をpとして、2回目からはk倍ずつ確率が上がっていく(kp,k²p,k³p…）とすると、私の結果はp=0.0089, k=2.4になった。(2022.05.27)
- 高周波回路シミュレータQucsStudioを使ってみる（その５）マルチポート（４ポートだとs4pみたいな）のSパラメータデータをTouchstoneファイルに出力しようとするとExport to SnPでは2ポートになってしまう。Project→Import Dataからポート数を手で入れないと。(2022.05.24)
- 新型コロナウイルス、日本の陽性者数＆ワクチン接種者数総計をプロット＆中国、韓国、アメリカ、ドイツ、フランス、イギリスの陽性者数もプロット（5/22更新）どの国もじわじわと増加は続いている。マスクをしなくなったり油断してるからかな…(2022.05.22)
- 高周波（RF・マイクロ波・ミリ波・5G）関連ニュース2022年5月19日 IEEE Microwave Magazineで磁性体を使わない非可逆デバイス、Microwave Journalは5Gスモールセル、THz、Movanoがミリ波血圧＆血糖値計発表、三菱の3Dプリンタ衛星アンテナ、など。(2022.05.19)
- 新型コロナウイルス、日本の陽性者数＆ワクチン接種者数総計をプロット＆中国、韓国、アメリカ、ドイツ、フランス、イギリスの陽性者数もプロット（5/15更新）韓国も中国も減速だが日本はじわじわと増え続けてる。アメリカも増え始めた。(2022.05.15)
トラックバック
この記事へのトラックバック一覧です： 1223個の新しい三体問題の周期解が見つかったということで、私もExcel VBAにルンゲクッタ8次のDOP853を移植したもので計算してみた。なかなか面白い動きをする。:
« #せやねん メチャ売れのコーナー、9/23の特集はエースコックの新食感焼そば モッチッチ。 | トップページ | 味仙監修の名古屋名物台湾ラーメンのカップ麺を食べた。なかなかの辛さでした。 »
コメント