Runge-Kutta solver (matrix form)

"run_rk4" uses the fourth order Runge-Kutta numerical method to solve nth order nonlinear and linear ODE systems expressed

This function is a custom implementation fourth order runge-kutta numerical method to solve nth order nonlinear and linear ordinary differential equations (ODE) systems which can be expressed as follows :

where F(x) represents the system's equation in a matrix form. e.g:

dxdt = [

10*(x(2)-x(1));

28*x(1)-x(2)-x(1)*x(3);

x(1)*x(2)-(8/3)*x(3)

];