Computational Methods for Kepler's Law

The purpose is to find solutions to Kepler's equation via python to find the position, velocity, or mass of the two celestial bodies given as a function of time. The first part employs the relaxation method to solve and plot Kepler's equation to find a value of eccentric anomaly (E) for a given value of mean anomaly (M) when the orbital eccentricity (e) is less than 1. The second portion involved an orbital eccentricity greater than 1. The third part of the program found the orbit of the Great Comet of 1680, employing the methods from the first case where the orbital eccentricity was less than one, where the orbital eccentricity and total period to orbit the sun are given. Results shown via graphs on the poster.