rk.m

来自「仿真 三点法 matlab 数值仿真 机动目标弹道仿真」· M 代码 · 共 28 行

%龙格库塔子函数，n:微分方程的个数；h:积分步长

function rk(n,h)

global y;

a(1)=h/2;
a(2)=a(1);
a(3)=h;
a(4)=h;

dy=grkt2f(y);
 
for i=1:n;
	old_y(i)=y(i);
end

for j=1:3;
	for i=1:n;
        y1(i)=old_y(i)+a(j)*dy(i);
		y(i)=y(i)+a(j+1)*dy(i)/3.;
    end
    dy=grkt2f(y1);
end
     
for i=1:n;
    y(i)=y(i)+h*dy(i)/6.;
end