%To use Runge Kutta is easy.  Here are the Lorenz equations.

[t,y]=solveRK4(@lorenzeq,[0,20],randn(3,1),0.001);
plot3(y(:,1),y(:,2),y(:,3))