昴問題（7体問題）をシンプレクティック8次で計算してGIFアニメにしてみる。
先日、Scratchでピタゴラスの3体問題を図示するものをアップロードしたら、こういう動きのあるものは結構海外の方からコメント付いたりするので気を良くして今度は昴問題をやってみよう。その前準備としてシンプレクティック8次で計算してGIFアニメ化。
昴(the Pleadis)問題は、”常微分方程式の数値解法I”に出ている例で、初期値として
(x,y) = (3,3),(3,-3),(-1,2),(-3,0),(2,0),(-2,-4),(2,4)に星があって、
(vx, vy) = (0,0),(0,0),(0,0),(0,-1.25),(0,1),(1.75,0),(-1.5,0)の初期速度であったときの運動。
ただし重さはmi=i（だんだん大きくなる）とする。
で計算結果はこちら。
星の間をすり抜けたり、なかなか面白い動きをする。
« え！スーパーマリオ64で遊ぶと脳が大きくなるという研究結果が！ | トップページ | 七五三の三角形をGeogebra4.2で描いて最大角度を確かめる。あと名古屋（758）と花さん（873）も。 »
「学問・資格」カテゴリの記事
- 高周波（RF・マイクロ波・ミリ波・5G）関連ニュース2020年9月17日 IEEE microwave magazine10月号はマイクロ波顕微鏡特集、 Microwave Journalはミリ波パッケージ。Vayyarの4Dイメージングレーダの仕組み、MIPIのA-PHY、など。(2020.09.17)
- 高周波（RF・マイクロ波・ミリ波・5G）関連ニュース2020年8月17日 Microwave Journal 8月号はドローンで衛星通信をテスト、220GHzの半導体テストシステム。DesignConで導体の粗さモデルと因果律について解説。Qualcommとワシントンポストの5Gクラス、Samsung Galaxy Watch3分解(2020.08.17)
- 高周波（RF・マイクロ波・ミリ波・5G）関連ニュース2020年8月10日 Microwave Magazine 9月号はシリコンフォトニクスをマイクロ波・ミリ波通信に使う、コロナと5G、IEEE RFICでimecが低消費電力60GHzレーダを発表、Qorvoの5G向けフィルタ解説、など。(2020.08.11)
- Python＋Scipyでルンゲクッタ8次のDOP853(Dormand＆Prince)を使う（その16) 時間依存のシュレーディンガー方程式のトンネル効果を計算してGIFアニメに。(2020.07.20)
- Python＋Scipyでルンゲクッタ8次のDOP853(Dormand＆Prince)を使う（その15) 蔵本モデルを計算してGIFアニメに。 (2020.07.18)
« え！スーパーマリオ64で遊ぶと脳が大きくなるという研究結果が！ | トップページ | 七五三の三角形をGeogebra4.2で描いて最大角度を確かめる。あと名古屋（758）と花さん（873）も。 »
コメント