rk4.m

来自「自己编写的变步长龙格库塔方法解强非线性微分方程」· M 代码 · 共 6 行

function y=rk4(t,y,h)
k1=myfun(t,y);
k2=myfun(t+h/2,y+h/2*k1);
k3=myfun(t+h/2,y+h/2*k2);
k4=myfun(t+h,y+h*k3);
y=y+h*(k1+2*k2+2*k3+k4)/6;