rk.m

把微分方程的一维经典初值问题的有限插分法：Euler法、改进的Euler法、Admas二步外插法、Admas二步内插法编写成通用子程序

%4级runge-kutta法编写成子函数
%函数的输入变量为：T为x的取值上限，n为把[0，T]区间等份的数目，u0为方程的初始值
function u=rk(T,n,u0)
h=T/n;
u=[];
u(1)=u0;
for i=1:n-1
    w1=f(i*h,u(i));
    w2=f(i*h+0.5*h,u(i)+0.5*w1*h);
    w3=f(i*h+0.5*h,u(i)+0.5*w2*h);
    w4=f(i*h+0.5*h,u(i)+0.5*w3*h);
    u(i+1)=u(i)+(1/6)*h*(w1+2*w2+2*w3+w4);
end