%admas二步內(nèi)插法求解微分方程的數(shù)值解.
%函數(shù)的輸入變量為：T為x的取值上限，n為把[0，T]區(qū)間等份的數(shù)目，u0為方程的初始值
function u=fadams(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*f(h*i,u(i));
    u(i+1)=u(i)+(h/12)*(5*f(h*(i+1),u(i+1))+8*f(h*i,u(i))-f(h*(i-1),u(i-1)));
end