KCTF and EPS

Fabrycky & Tremaine の計算はやはり合っていたみたい。精度を上げて EPS=1e-15 くらいにしても、潮汐が効いてくるにつれて e1min が上がってくる現象は見られる。あの２人がそう間違える筈は無いわな。

しかしこのコードは潮汐が効きだしてからやたらと遅くなる。やはり潮汐の式が stiff 過ぎるからか。