《数值分》第五讲：常微分方程数值解.pptVIP

参考程序-Euler x=0:0.01:1; y=sqrt(1+2.*x); a=0.0;b=1.0;n=10; h=(b-a)/n; x0=a:h:b;y0(1)=1.0; for k=1:10 y0(k+1)=y0(k)+h*(y0(k)-2*x0(k)/y0(k)); end for i=1:10 y1(1)=1.0; y1(i+1)=y1(i)+h*(y1(i)-2*x0(i)/y1(i)); y1(i+1)=y1(i)+h*(y1(i)-2*x0(i)/y1(i)+y1(i+1)- 2*x0(i+1)/y1(i+1))/2; end plot(x,y,b); hold on; plot(x0,y0,or); hold on; plot(x0,y1,*); 第五章：常微分方程数值解 参考程序-Lunge_Kutta x=0:0.01:1; y=sqrt(1+2.*x); a=0.0;b=1.0;n=10; h=(b-a)/n; x0=a:h:b; y0(1)=1.0; for k=1:10 k1=y0(k)-2*x0(k)/y0(k); k2=y0(k)+h*k1/2-(2*x0(k)+h)/(y0(k)+h*k1/2); k3=y0(k)+h*