rk4.m

该文件中包含有常微分方程初值问题的四种数值解法

function RK4(h,x0,xn,y0)
N=(xn-x0)/h;
N=round(N);
x=zeros(1,N+1);
y=zeros(1,N+1);

for i=0:N
    x(i+1)=x0+h*i;
end
y(1)=y0;

for i=1:N
    k1=f(x(i),y(i));
    k2=f(x(i)+h/2,y(i)+h*k1/2);
    k3=f(x(i)+h/2,y(i)+h*k2/2);
    k4=f(x(i)+h,y(i)+h*k3);
    y(i+1)=y(i)+(k1+2*k2+2*k3+k4)*h/6;
end
    

fprintf('\n用RK4方法求解常微分方程初值问题：\n');
for i=0:N
    fprintf('X(%d)=%f,  Y(x%d)=%f\n',i,x(i+1),i,y(i+1));
end