rk3.txt

来自「龙哥库塔4阶函数的计算」· 文本 代码 · 共 17 行

function R=rk4(f,t,y0,h)

M=t/h;
T=zeros(1,M+1);
Y=zeros(1,M+1);
T=0:h:t;
Y(1)=y0;
for j=1:M
   k1=h*f(T(j),Y(j));
   k2=h*f(T(j)+h/2,Y(j)+k1/2);
   k3=h*f(T(j)+h,Y(j)-k1+2*k2);
   Y(j+1)=Y(j)+(k1+4*k2+k3)/6;
end

R=Y(M+1);

return