Newtonian gravity simulation using Dormand-Prince Runge-Kutta ODE solver in Python

The video shows a simulation of Newtonian gravity, produced using a Runge-Kutta Ordinary Differential Equation (ODE) solver. An adaptive step size Dormand-Prince method was used in the code, which I wrote myself using Python. The code is freely available at: I am new to Python, so sorry to all those experienced Pythoners, but the code reads a bit like MATLAB. The simulation itself just uses a point particle approximation, the finite siz