Question

急！！！寻数学专业的MATLAB高手帮写一道题！

要求学过计算方法这门课，懂得用MATLAB语句写出四阶Runge-Kutta方法求解微分方程组（两个方程）的程序，不能直接用数值解的语句ode来写，注意！！！一定是用四阶Runge-Kutta方法。
题目是蛮简单的，主要是我MATLAB学得不好，而且最近忙考试，没时间写。只有一道题，觉得没问题的朋友请留下邮箱，我会将题目发到你的邮箱，完成以后将程序代码发给我，回答满意可追加分，万分感谢！！！

Answer 1

function y = DELGKT4_lungkuta(f, h,a,b,y0,varvec)
format long;
N = (b-a)/h;
y = zeros(N+1,1);
y(1) = y0;
x = a:h:b;
var = findsym(f);
for i=2:N+1
K1 = Funval(f,varvec,[x(i-1) y(i-1)]);
K2 = Funval(f,varvec,[x(i-1)+h/2 y(i-1)+K1*h/2]);
K3 = Funval(f,varvec,[x(i-1)+h/2 y(i-1)+K2*h/2]);
K4 = Funval(f,varvec,[x(i-1)+h y(i-1)+h*K3]);
y(i) = y(i-1)+h*(K1+2*K2+2*K3+K4)/6;
end
format short;

以前下过，代码挺简单的～～