This general coupled differential equation solver uses an adaptive Runge-Kutta method. The time step size used to solve the coupled differential equations is calculated according to a real-time accuracy estimate. - View it on GitHub