Euler方法以及Euler改进.doc

Euler方法以及Euler改进 方法的主程序： format long fun x+y; a -0.2;b 0.2; h 0.02; n b-a /h;X a:h:b; Y zeros 1,n+1 ; X 1 a;Y1 1 1;Y2 1 1; Y3 1 1; %Euler方法 for i 1:n x X i ;y Y1 i ; Y1 i+1 Y1 i +h*eval fun ; end %Euler 预测-校正方法 X 1 a;Y2 1 1; for i 1:n x X i ;y Y2 i ; YP Y2 i +h*eval fun ; YC Y2 i +h*eval fun ; Y2 i+1 YP+YC /2 end %计算精确值 fun1 2*exp x -x-1; for i 1:n x X i+1 ; Y3 i+1 eval fun1 end [X,Y1,Y2,Y3] %显示结果 %绘制结果图形 fplot fun1,[a,b] hold on plot X,Y1,*r,X,Y2,+m hold on plot X,Y3,ok title a,b取值为[-0.2，0.2] return 当a，b取不同的值的时候： [0,0.1] ans 0 1.000000000000000 1.000000000000000 1.000000000000000 0.0座机电话号码000000 1.0座机电话号码000000 1.0座机电话号码000000 1.0座机电话号码053512 0.0座机电话号码000000 1.0座机电话号码000000 1.0座机电话号码000000 1.0座机电话号码384776 0.0座机电话号码000000 1.0座机电话号码000000 1.0座机电话号码816000 1.0座机电话号码090719 0.0座机电话号码000000 1.0座机电话号码000000 1.0座机电话号码257283 1.0座机电话号码349917 0.1XXXXXXXXXX0000 1.10座机电话号码00000 1.110座机电话号码2880 1.110座机电话号码1295 [-0.1,0.1] ans -0.1XXXXXXXXXX0000 1.000000000000000 1.000000000000000 1.000000000000000 -0.0座机电话号码000000 1.01座机电话号码00000 1.01座机电话号码00000 0.座机电话号码 -0.0座机电话号码000000 1.0座机电话号码000000 1.0座机电话号码000000 0.座机电话号码 -0.0座机电话号码000000 1.0座机电话号码000000 1.0座机电话号码575200 0.座机电话号码 -0.0座机电话号码000000 1.0座机电话号码000000 1.0座机电话号码644420 0.座机电话号码 0 1.0座机电话号码080000 1.0座机电话号码860237 1.000000000000000 0.0座机电话号码000000 1.11座机电话号码01600 1.1座机电话号码128214 1.0座机电话号码053512 0.0座机电话号码000000 1.1座机电话号码533632 1.1座机电话号码257804 1.0座机电话号码384776 0.0座机电话号码000000 1.1座机电话号码904305 1.1座机电话号码747611 1.0座机电话号码090719 0.0座机电话号码000000 1.1座机电话号码382391 1.1座机电话号码707713 1.0座机电话号码349917 0.1XXXXXXXXXX0000 1.座机电话号码 1.座机电话号码 1.110座机电话号码1295 [0,0.2] ans 0 1.000000000000000 1.000000000000000 1.000000000000000 0.0座机电话号码000000 1.0座机电话号码000000 1.0座机电话号码000000 1.0座机电话号码053512 0.0座机电话号码000000 1.0座机电话号码000000 1.0座机电话号码000000 1.0座机电话号码384776 0.0座机电话号码000000 1.0座机电话号码000000 1.0座机电话号码816000 1.0座机电话号码