You first turn the system into a possibly larger system of first order ordinary differential equations.

Then you use an ODE solver (like ode45) for the resulting vector first order system.

(there are some examples in this part of MHF showing this so just search for the examples)

.