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の特集はエースコックの新食感焼そば モッチッチ。 | トップページ | 味仙監修の名古屋名物台湾ラーメンのカップ麺を食べた。なかなかの辛さでした。 »
「学問・資格」カテゴリの記事
- iPhoneで使える数学ソフトMaple Calculatorで手書き入力でMathematicaの常微分方程式の例を試す（その１）(2021.03.29)
- 高周波（RF・マイクロ波・ミリ波・5G）関連ニュース2021年3月16日 IEEE Microwave magazineの特集は中国のマイクロ波技術、Microwave Journalの特集はCATRでミリ波ビームスキャンを測定、SIJの特集はなんと450ドルのネットアナ、フリーの高周波シミュレータ（Qucsなど）。(2021.03.16)
- 高周波（RF・マイクロ波・ミリ波・5G）関連ニュース2021年2月16日 IEEE Microwave Magazineの特集はオールデジタルのRFID、Microwave JournalはＥバンド ミリ波通信に衛星や気球を使う話、アメリカの半導体企業がバイデンに投資を迫る、(2021.02.17)
- カオスを生じる電気回路、Chua’s circuitをLTspiceで回路シミュレーションしてみる。(2021.02.19)
- Labyrinth Chaos(迷宮カオス）を生むThomas-Rössler方程式のパラメータbを色々変えて、Python＋Scipyでルンゲクッタ8次のDOP853(Dormand＆Prince)を使って計算してGIFアニメ(2021.02.16)
「日記・コラム・つぶやき」カテゴリの記事
- 新型コロナウイルス、中国、日本、韓国、アメリカ、ドイツ、フランス、イギリスでの感染者数を指数関数＆ロジスティック関数＆Log-Logプロットでべき関数フィッティングした（4/11更新）日本は陽性者50万人を超え第四波の始まり、というかイギリス以外は全部増えだしてる。(2021.04.12)
- ３乗すると（奇数回べき乗すると）元に戻る数59161727622001114846846461792218008213239954784512519836425781249（２乗すると元に戻る数はペレリマン数列）(2021.04.11)
- 円と接線を使った1枚の絵でsin/cos/tan/sec/csc/cotを表すもの、複素数を使って接線のパラメータ表示(z=exp(iθ)+i exp(iθ) t )で計算してみる。(2021.04.05)
- 新型コロナウイルス、中国、日本、韓国、アメリカ、ドイツ、フランス、イギリスでの感染者数を指数関数＆ロジスティック関数＆Log-Logプロットでべき関数フィッティングした（4/4更新） フランス・ドイツの増え方がエグイが日本・韓国もやはり増加中。イギリスの減速がすごいな。(2021.04.05)
- 分数階微分(Fractional Calculus)の計算をPythonのdifferint（Riemann-Liouvilleの定義で計算するもの）でやってみる。(2021.04.02)
トラックバック
この記事へのトラックバック一覧です： 1223個の新しい三体問題の周期解が見つかったということで、私もExcel VBAにルンゲクッタ8次のDOP853を移植したもので計算してみた。なかなか面白い動きをする。:
« #せやねん メチャ売れのコーナー、9/23の特集はエースコックの新食感焼そば モッチッチ。 | トップページ | 味仙監修の名古屋名物台湾ラーメンのカップ麺を食べた。なかなかの辛さでした。 »
コメント