3body2019927red2.nb[setting] time measured by (G Msun)/AU^3 i.e. Year, distance by AU, mass by Mjupiter radius \[Epsilon][i] for each planet i \[Epsilon][i,j] mutual min. distance r0 initial revolution radius rd0 ratio of the initial separation to hill radius for the most heavy planet whenevent distance > 10 r0, stop validity check by energy conservation
いくつかのラン: スキームはクランクニコルソン(mathematica default) 1000年走らせて,途中で初期半径より10倍大きくなるものがあれば飛んだとみなし計算を止める. エネルギー誤差(2番目のグラフ)は大体0.1%以下.