The only equation used is F= -GMm/r[sup]2 [/sup]
I do not know what you would like me to explain?
Do you mean the numerical method? (I used Runge-Kutta 4th order method).