数学软件习题.doc

数学软件习题 1、用前进Euler公式求解并和精确解比较： （1）， 由题可知：， M文件： ： function y = a1(h) a=0; b=2; y0=1; N = (b-a)/h; y = zeros(N+1,1); x = a:h:b; y(1) = y0; for i=2:N+1 y(i) = y(i-1)+2*h*y(i-1)*((i-1)*h); end ： x = 0:0.1:2; z=zeros(21,1); for i=1:21 z(i) = exp(x(i)^2); end y1=a1(0.1); y2=a1(0.01); y3=a1(0.001); y22=zeros(21,1); for i=1:21 y22(i) = y2(10*(i-1)+1); end y33=zeros(21,1); for i=1:21 y33(i) = y3(100*(i-1)+1); end plot(x,y1,b-o,x,y22,g-x,x,y33,k-+,x,z,r); legend(h=0.1,h=0.01,h=0.001,exact); xlabel(x) ylabel(exact solution and approximate solutions) ： e1=zeros(21,1); for i=1:21 e1(i) = abs(z(i)-y2(10*(i-1)+1)); end e2=zeros(21,1); for i=1:21 e2(i) = abs(z(i)-y2(10*(i-1)+1)); end e3=zeros(21,1); for i=1:21 e3(i) = abs(z(i)-y3(100*(i-1)+1)); end plot(x,e1,b-o,x,e2,g-x,x,e3,k-+); legend(h=0.1,h=0.01,h=0.001); xlabel(x) ylabel(绝对误差) 命令窗口： a2 a3 X Z Y1 E1 Y2 E2 Y3 E3 0.0 1.000000 1.000000 0.000000 1.000000 0.000000 1.000000 0.000000 0.1 1.010050 1.020000 0.009950 1.011053 0.001003 1.010150 0.000100 0.2 1.040811 1.060800 0.019989 1.042835 0.002024 1.041013 0.000203 0.3 1.094174 1.124448 0.030274 1.097255 0.003081 1.094483 0.000309 0.4 1.173511 1.214404 0.040893 1.177695 0.004184 1.173930 0.000419 0.5 1.284025 1.335844 0.051819 1.289360 0.005334 1.284560 0.000535 0.6 1.433329 1.496146 0.062816 1.439841 0.006512 1.433983 0.000653 0.7 1.632316 1.705606 0.073290 1.639974 0.007658 1.633085 0.000769 0.8 1.896481 1.978503 0.082022 1.905130 0.008649 1.897350 0.000869 0.9 2.247908 2.334633 0.086725 2.257150 0.009242 2.248838 0.000930 1.0 2.718282 2.801560 0.083278 2.727268 0.008986 2.719187 0.000905 1.1 3.353485 3.417903 0.064419 3.360546 0.007061 3.354197 0.000712 1.2 4.220696 4.238200 0.017504 4.222698 0.002002 4.220898 0.000202 1.3 5.419481 5.340132 0.079349 5.410677 0.008803 5.418590 0.000891 1.4 7.099327 6.835369 0.263958 7.069330 0.029997