As you've not provided the detail, I suppose it's a system of nonlinear ODE with initial conditions, IVP. Therefore I'll do the followings:
- Convert the system of equations to state variable space
- Code and test a function to calculate the RHS of above equations getting input data
- Code the main script to solve the equations, print the numerical outputs (if required) and plot the required graphs
If it was not ODE system or IVP, then I will change the steps of my proposal.
I am the author of book "Numerical Recipes for Solving Engineering Problems using MATLAB," published on 2002.