USER REVIEWS
Since the eccentric anomaly cannot be directly calculated from certain given values, this function performs a simple Newton-Raphson iteration to solve for the eccentric anomaly within a given tolerance (default 10^-8 radians)
WRITE USER REVIEW