月亮近月点收敛性有些问题,按以下修改即可
把eph0.js中的//初算二次dt=1/36525;r1 = XL1_calc( 2, t-dt, 10 );改为//初算二次dt=2/36525;r1 = XL1_calc( 2, t-dt, 10 );这样就可以较可靠收敛了