steprk.m

来自「这是比例导引的MATLAB程序,其中龙格库塔积分的程序很简单且通用性强」· M 代码 · 共 13 行

%函数说明，t表示自变量（一般为时间），y表示变量的初值（是一个行向量），h表示步长
%运行的结果，也是一个行向量

function dy = stepRK(t,y,h)
K1=rightF(t,y);
K2=rightF(t+h/2,y'+h*K1/2);
K3=rightF(t+h/2,y'+h*K2/2);
K4=rightF(t+h,y'+h*K3);
dy=y'+(K1+2*K2+2*K3+K4)*h/6;
dy=dy';