Excelで三体問題をシンプレクティック8次で解いてみる。
この前、8次のSymplectic公式の数表を書いたので、せっかくなんで三体問題でも解いてみよう。
重力ポテンシャルはU= -G Mi Mj / |Ri-Rj|
でめんどくさいんでG=M=1だ。で半径1の円の上に120°ずつ星を並べて計算。すると…
きれいに知恵の輪のようになった。五体問題にすればオリンピックになるのかな(初期条件適当に与えればできるのかもしれないがめんどくさい。)
次はフェルミ・パスタ・ウラムの問題かな。あの時代のアルゴリズムって何使ってたんだろ。エネルギーのことがからむし。多分leap frogじゃないかなあ。フェルミがいてエネルギーが保存しないスキーム使うはずないしな。
« リアルキアイダーを飲んでみた。 | トップページ | これは面白いのか(シェークスピアと言う名のプログラミング言語) »
「学問・資格」カテゴリの記事
- 高周波(RF・マイクロ波・ミリ波・5G)関連ニュース2021年2月16日 IEEE Microwave Magazineの特集はオールデジタルのRFID、Microwave JournalはEバンド ミリ波通信に衛星や気球を使う話、アメリカの半導体企業がバイデンに投資を迫る、(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)
- フィッツヒュー・南雲 (FitzHugh-Nagumo) 方程式をPython+Scipyでルンゲクッタ8次のDOP853(Dormand Prince)で計算。(2021.02.23)
- 「水晶振動子の等価回路計算」をカシオの高精度計算サイトkeisan.casio.jpの自作式としてUP! インピーダンスの大きさと位相がグラフ化できる。(2021.02.12)
« リアルキアイダーを飲んでみた。 | トップページ | これは面白いのか(シェークスピアと言う名のプログラミング言語) »
コメント