deletions | additions       diff --git a/Method.md b/Method.md  index e6133cc..e88c9c7 100644  --- a/Method.md  +++ b/Method.md 
...
Every 50 time steps, the program output the kinetic, potential, and total energy of the system, as a function of time, calculated via a separate subroutine and called upon in the main Do Loop.   ##Data Analysis  To study the evolution of the proto-solar system model, I analyzed how the orbital parameters of the particles and Neptune changed over the course of the simulation. I particularly focused on the orbit parameters of Neptune, which could be studied quantitatively given the data output by the FORTRAN code. The FORTRAN code output calculated raw data, which was plotted and further analyzed using two IDL codes. The first IDL code, orbit_parameters.pro (see attached code), took Neptune orbital radius data as its input. In addition to simply plotting the orbital radius as a function of time, this program used the apoapsis and periapsis of every orbit to calculate semi-major axis $a$, orbital period $p$, and eccentricity $e$ as a function of time, via the following formulae:  \begin {equation}  a=(r_apoapsis+r_periapsis)/2.  \end {equation}  \begin {equation}  p=2*(t_apoapsis-t_periapsis)  \end {equation}  \begin {equation}  e=(r_apoapsis-r_periapsis)/(r_apoapsis+r_periapsis)  \end {equation}