To-do's

1. GR 無しでの e_max vs a1 の計算もやってみる。

perturber が 15 AU で 0.5 木星質量のケースでは、摂動の一周期に 30Myr くらいかかっている。これでは、一度目の e_max に辿り着く前に migration が終わってしまう。

2. Derive octupole-order equations for dI/dt, d(node)/dt

3. GR も含めた e_max の式と、derivation を読む (Fabrycky & Tremaine)