四阶龙格-库塔法.txt

数值计算方法的一些题

function y=my4_Runge_Kutta(fun,a,b,n,u)

%u为函数初值。

%n为区间等分数。

%a,b分别为上下界。

h=(b-a)/n;

t=a;

y=u;

for i=1:n

    K1=fun(t,y);

    K2=fun(t+h/2,y+K1*h/2);

    K3=fun(t+h/2,y+K2*h/2);

    K4=fun(t+h,y+h*K3);

    y=y+(K1+2*K2+2*K2+K4)/6;

    t=a+i*h;

end