eccentricity:とでるので、計算したい軌道の離心率をいれる。epsinit:とでる。初期には低次の公式から出発する。分小さいステップにする必要がある。epsinit/|加速度の導関数の絶対値|としているので、それで大丈夫な程度に小さくする。testinで実行できる。位置、速度は現在のところ上のような感じに画面にでる。そのへんは適宜ソースファイルをいじるか、この出力を加工すること。違う初期条件にするならそれもファイルから読むようにするとか改造すること。2017/10/19。