Celestial body simulation with Runge-Kutta 4