Scratchで新しい三体問題の周期解（1223個）を計算してみる、、、計算できたが遅すぎる。8の字解を計算したものを代わりに。。。
この前、三体問題の周期解が1223個も見つかったという論文を見て、、、
ExcelのVBAにルンゲクッタ8次のDOP853 (Dormand＆Prince)を移植したもので計算したりしていたが！
今度はScratchで計算してみよう。II.D2（2）というやつだが、、、
計算はできたんですよ、計算は。
でもDOP853はさすがに移植できず、シンプレクティックオイラー法を使ったら
刻み幅をとんでもなく小さくしないとすぐに発散する。
なので、ターボモードですら遅すぎる、、、
そのためちょっとアップロードは保留。
代わりに8の字解を計算したものを再掲。こっちは速い！
« 9/24 20時過ぎからJRの三ノ宮駅に人が立ち入り、JR神戸線と京都線が遅れてる！ | トップページ | 35段14次のルンゲクッタ法をPARI/GPに実装、ローレンツ方程式を計算し、通常の4次、そして8次、オイラー法と比較してみる。（再掲） »
「パソコン・インターネット」カテゴリの記事
- 高周波回路シミュレータQucsStudioを使ってみる（その３）Mixed Mode S parameterを計算(2018.10.12)
- 高周波回路シミュレータQucsStudioを使ってみる（その２）SパラメータのTouchStoneフォーマットで出力するには？(2018.10.11)
- 高周波回路シミュレータQucsStudioを使ってみる（その１）まずは何をさておきμの文字化けだけには注意。(2018.10.10)
- 円の弧長,弦長,矢高,半径のどれか2つを与えて残りを計算（カシオの高精度計算サイト自作式）で180°以上、複数解に対応。(2018.10.09)
- macbook proをmacOS Mojaveにアップデート。せっかくなんでダークモードにしてみる。(2018.09.26)
「学問・資格」カテゴリの記事
- 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)
「日記・コラム・つぶやき」カテゴリの記事
- 円と接線を使った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)
- 新型コロナウイルス、中国、日本、韓国、アメリカ、ドイツ、フランス、イギリスでの感染者数を指数関数＆ロジスティック関数＆Log-Logプロットでべき関数フィッティングした（3/28更新） イギリスが目に見えて減速だがフランスドイツが上昇。日本も微増傾向。(2021.03.29)
- 3/14(日本時間は昨日）は円周率Πの日で、Google検索の電卓機能でπを何桁覚えておけるかのゲームができる！(2021.03.15)
« 9/24 20時過ぎからJRの三ノ宮駅に人が立ち入り、JR神戸線と京都線が遅れてる！ | トップページ | 35段14次のルンゲクッタ法をPARI/GPに実装、ローレンツ方程式を計算し、通常の4次、そして8次、オイラー法と比較してみる。（再掲） »
コメント