feuler.m

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

%改进的euler法解一维微分方程的通用子函数
%函数的输入变量为：T为x的取值上限，n为把[0，T]区间等份的数目，u0为方程的初始值
function u=feuler(T,n,u0)
h=T/n;
u=[];
u(1)=u0;
for i=1:n-1
    u(i+1)=u(i)+h*f(h*i,u(i));
    u(i+1)=u(i)+(h/2)*(f(h*(i+1),u(i+1))+f(h*i,u(i)));
end