昴問題（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）も。 »
「学問・資格」カテゴリの記事
- #仮面ライダービルド 第40話の黒板の数式は、、、一瞬しか映らずわからなかった、、、ハミルトニアンと、L,Rが見えたので量子もつれ関連の何か？→全然違いました。密度行列繰り込み群か！(2018.06.17)
- #仮面ライダービルド 第40話の話数を表す数式はペル方程式x^2-77y^2=1の最小解がx=351,y=40ということ。(2018.06.17)
- #AKB総選挙 2018年の結果がジップの法則（というかべき乗則）に乗るかを検証。一位は松井珠理奈さんでした。3位、7位、15位で関数形が変わる。(2018.06.16)
- Powerpointで数式入力をLaTeX形式で行う。WordはデフォルトでできるけどPowerpointはちょっと初期設定が必要。(2018.06.13)
- #仮面ライダービルド 第39話の黒板の数式はワインバーグ・サラム理論（＆ヒッグス機構）？(2018.06.10)
トラックバック
この記事のトラックバックURL:
http://app.cocolog-nifty.com/t/trackback/512682/58527360
この記事へのトラックバック一覧です： 昴問題（7体問題）をシンプレクティック8次で計算してGIFアニメにしてみる。:
« え！スーパーマリオ64で遊ぶと脳が大きくなるという研究結果が！ | トップページ | 七五三の三角形をGeogebra4.2で描いて最大角度を確かめる。あと名古屋（758）と花さん（873）も。 »
コメント