真近点角νの求める解の精度を、ここではEps＝1e-10とします。 ケプラーの方程式は、S=tan ν/2 と置けば、M=S(１＋pow(S,２)/３) になります。 ここで、漸近法により、初期値So＝Mとして、Sn+1＝(M＋２pow( Sn ,３)/３)/(１＋pow(Sn ,２)) を、Sn+1−Sn＜Eps …