yugu.m

来自「Adams四阶预测的源码」· M 代码 · 共 24 行

function yu1=yugu(f,a,b,ya,m);

h=(b-a)/m;

t=zeros(1,m+1);

y=zeros(1,m+1);

t=a:h:b;

y(1)=1.0309
y(2)= 1.0608
y(3)=1.0896
y(4)= 1.1174

for j=1:m-4

y(j+4)=y(j+3)+h/24*(55*feval(f,t(j+3),y(j+3))-59*feval(f,t(j+2),y(j+2))+37*feval(f,t(j+1),y(j+1))-9*feval(f,t(j),y(j)));    
y(j+4)=y(j+3)+h/24*(19*feval(f,t(j+3),y(j+3))-5*feval(f,t(j+2),y(j+2))+feval(f,t(j+1),y(j+1))+9*feval(f,t(j+4),y(j+4)));

end

yu1=[t',y'];