untitled4.m

来自「这是我从一本基于MATLAB实现的数值分析算法的书中自带的光盘程序」· M 代码 · 共 14 行

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