Hermite Corrector II.

修正を一度だけかけると、a と e が少しずつ下がっていく。Kokubo, Yoshinaga & Makino (1998) の結果と同じ。

しかし修正を二度かけると、a も e も上がっていく。これは修正のどこかで計算ミスしている筈だ。

もう一度修正のかけ方を復習。

新しい加速度と、加速度の一次微分を使って、エルミート補完式から
a(2) と a(3) を求める。これが修正子。
これを、修正式を使って、予測子に足してやる。これで予測子 x_{c,i} と v_{c,i} が
求まる。



この予測子を使って、もう一度中心星からの加速度と一次微分を求め直す。