Here you can solve the Kepler's equation M = E-eSin(E) for an elliptical orbit. Given
mean anomaly M and eccentricity e , you can solve for eccentric anomaly E.
The ranges for e and M are [0,1] and [0,PI]. If M<0 , you can solve for E
with |M| and associate a negative sign to the resulting solution.