untitled4.m

利用MATLAB进行求解常微分方程组的四阶R-K算法法

clear;
%程序C8-2  四阶R-K法
t(1)=0;y(1)=1;h=1.0;n=1;
while t<=6
    K1=h*f_tm(y(n),t(n));
    K2=h*f_tm(y(n)+0.5*K1,t(n)+0.5*h);
    K3=h*f_tm(y(n)+0.5*K2,t(n)+0.5*h);
    K4=h*f_tm(y(n)+K3,t(n)+h);
    y(n+1)=y(n)+(K1+2*K2+2*K3+K4)/6;
    t(n+1)=n*h;
    n=n+1;
end
plot(t,y)
y