Matlab - Sistema de ecuaciones diferenciales

Hola, he tratado durante dias crear un programa con Rungge Kutta para solucioner un sistema de 3 ecuaciones y de segundo orden...Alguien tendra uno redactado parecido? Les coloco el que pude realizar que es solo de dos ecuaciones y de primer grado, espero poder partir de ahi.. Si tienen uno parecido bien sea de otro orden les agradeceria. Saludos

function [X1,Y1,T]=RK4Dosvariables
h=0.2;
n=2;
x=-1;
y=6;
%z=5
%xprima=1
%yprima=2
%zprima=3
Fx=inline('2*x+4*y','x','y');
Fy=inline('-x+6*y','x','y');
%Fz=inline(' ,x,y,z )

m1=Fx(x,y);
k1=Fy(x,y);;
m2=Fx(x+((m1*h)/2),y+((k1*h)/2))
k2=Fy(x+(m1*h)/2,y+(k1*h)/2)
m3=Fx(x+(m2*h)/2,y+(k2*h)/2)
k3=Fy(x+(m2*h)/2,y+(k2*h)/2)
m4=Fx(x+h*m3,y+h*k3)
k4=Fy(x+h*m3,y+h*k3)

%dx/dt
x1=x+(h/6)*(m1+2*m2+2*m3+m4)

%dy/dt
y1=y+(h/6)*(k1+2*k2+2*k3+k4)

x=x1;
y=y1;
T(n)=i+h;
n=n+1;
end

una sugerencia.
Primero haz para una sola ecuación el método de Runge Kutta
luego para dos ecuaciones y observa que modificaciones haz hecho
luego para tres ecuaciones siguiendo el mismo análisis para dos ecuaciones.
...
luego generalizas para n ecuaciones.