To convert a single second-order ODE to a system of two couples first-order ODEs, we create a new function equal to the derivative of our unknown function: we define . Then . So then our system is , or in matrix notation,

From here, you then perform the fourth-order Runge-Kutta on this system.

--Kevin C.