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%以下.