This is a VBA routine for Runge Kutta 4th order integration.
Function rhs(x, t, dxdt, p)
'==============================================’
' Input your right hand sides in the form dxdt(1) = expression
' in t, x(1), x(2), etc. The right hand side may also include
' several...