public static double GetEccentricAnomaly(double M, double e)
{
if (e < 1.0)
{
return(Kepler.NewtonElliptical(M, e, 1));
}
if (e > 1.0)
{
return(Kepler.NewtonHyperbolic(M, e));
}
if (e == 1.0)
{
return(Kepler.NewtonParabolic(M, e, 1));
}
return(M);
}