ode solver code snippets in matlab

create ode solver for 2 dof in matlab