Busque ayuda de expertos en programación y utilice la programación de Matlab para resolver ecuaciones diferenciales utilizando el método Longo Kuta.
Función: Utiliza el método de Runge-Kutta de cuarto orden para resolver ecuaciones diferenciales ordinarias
-------------------- --- -----------------------
Función R=Rungkuta4(f, a, b, n, ya) p>
%.f: el identificador de función en el extremo derecho de la ecuación diferencial
% a, b: los dos puntos finales del intervalo de valor de la variable independiente
% n : el número de partes iguales del intervalo
p>
% ya: función y(a)
% R=[x',y']: matriz compuesta por la suma de la variable independiente X y la variable independiente Y
h=(b-a)/n;
x=zeros(,n+1);
y =ceros(1,n+1);
x=a:h:b;
y(1)=ya;
para i=1 :n
k1=h*feval(f , x(i),y(i));
k2=h*feval(f, x(i)+h/ 2,k2=h x(i)+h/2,y(i )+k1/2
k3=h*feval(f,x(i)+h/2,y(i); )+k2/2);
k4= h*feval(f,x(i)+h,y(i)+k3);
y(i+1) =y(i)+(k1+2*k2+2*k3+ k4)/6;
fin
R=[x',y'];
y(i+1)=y(i)+( k1+2*k2+2*k3+ k4)/6