Integrate a system of ODEs using the Second Order Runge-Kutta (Midpoint) method
github.com/scijs/integrate-adaptive-simpson
scijs/integrate-adaptive-simpson