adams.m

来自「把微分方程的一维经典初值问题的有限插分法：Euler法、改进的Euler法、Ad」· M 代码 · 共 11 行

%admas二步外插法求解微分方程的数值解.
%函数的输入变量为：T为x的取值上限，n为把[0，T]区间等份的数目，u0为方程的初始值
function u=adams(T,n,u0)
h=T/n;
u=[];
u(1)=u0;
u(2)=u(1)+h*f(h,u(1));
for i=2:n-1
    u(i+1)=u(i)+(h/2)*(3*f(h*i,u(i))-f(h*(i-1),u(i-1)));
end