Python+Scipyでルンゲクッタ8次のDOP853(Dormand Prince)を使う（その9） 三体問題の周期解いろいろ（3） 1223個見つかった解のうち、I.A68(0.5),I.A2(4),I.B59(0.75)

今回は、さらに見つかった1223個の解について：

1223個の新しい三体問題の周期解が見つかったということで、私もExcel VBAにルンゲクッタ8次のDOP853を移植したもので計算してみた。なかなか面白い動きをする。

初期値は
x0=[-1,0.2138410831, 0,0.0542938396,1,0.2138410831, 0,0.0542938396,0,-2*0.2138410831/0.5, 0,-2*0.0542938396/0.5]

で、m1=m2=1, m3=1/2とする。

計算結果はこちら：クリックするとGIFアニメが始まります。

先日よりさらにすごい動き。

次は、I.A2(4)

v1=0.9911981217
v2=0.7119472124
m1=1
m2=1
m3=4
x0=[-1,v1, 0,v2,1,v1, 0,v2,0,-2*v1/m3, 0,-2*v2/m3]

 

つぎは？I.B59(0.75)

v1=0.4101378717
v2=0.1341894173
m1=1
m2=1
m3=0.75
x0=[-1,v1, 0,v2,1,v1, 0,v2,0,-2*v1/m3, 0,-2*v2/m3]