Added a subroutine to calculate the true anomaly using Newton's method