It depends on the accuracy you need, and the order of your DEs. Note that, all of the numerical methods could be applied in MATLAB, such as "Leap Frog".
u can use MATLAB inbulit functions like ode45,ode23 etc.. some of these uses rungekutta family of methods.... But they can handle most types of ODEs stiff, nonstiff etc... check this link for more..
Basically I am trying to solve a problem related to plasma in sheath region. I also got my result using ode45 in matlab. But the fact is that I wanted to check it by using different kind of approach or method to check the accuracy and variation in the result. Thanks to everyone I got what I wanted....
Solving differential equation is a lot easier in Maple. You should give it a try. There is no need to define variables and all. You can directly write the equation and find the solution.