adams方法探讨(精品·公开课件).pptVIP

精确解matlab编程： function g=f(x,y) g=1/(1+x^2)-2*y^2; end h=0.1； h=0.1;y(101)=0;x=0:h:10; for i=1:3 k1=h*f(x(i),y(i)); k2=h*f(x(i)+h/2,y(i)+k1/2); k3=h*f(x(i)+h/2,y(i)+k2/2); k4=h*f(x(i)+h,y(i)+k3); y(i+1)=y(i)+(k1+2*k2+2*k3+k4)/6; end for i=1:97 y(i+4)=y(i+3)+h*(55*f(x(i+3),y(i+3))-59*f(x(i+2)... ,y(i+2))+37*f(x(i+1),y(i+1))-9*f(x(i),y(i)))/24; end g=x./(1+x.^2);plot(x,y,+,x,g,-) 四步四阶adams显式matlab程序（龙格库塔法解出前四个点）： 用单步欧拉法达到h的四次方精度（即用欧拉法解出前四个点）： h=0.1;y(101)=0;x=0:h:10;t=h^4; x(2)=t;x(3)=2*t;x(4)=3*t; for i=1:3 y(i+1)=y(i)+t*f(x(i),y(i)); end for i=1