计算方法试验10

1、9实验十一、编程并在计算机上调试修改运行1. 利用公式计算2. 利用公式计算oa=zndn xn!sin x 八ndD2n 1x(2n 1)!二、阅读理解下列程序，并在计算机上运行1. jzhexp.m%利用&的展开式计算矩阵函数eA的值fun cti on E,k=jzhexp(A);%A是方阵E=zeros(size(A);F=eye(size(A);k=1;while norm(F,1)>0E=E+F;F=A*F/k;k=k+1;end2. xjfdd .m % x0,x为初值，w为精度，n为最大迭代次数function a,k=xjfdd(x0,x,w,n)a(1:2)=

2、x0,x;for k=1:nif f(x)-f(x0)=0disp('分母为零')breakend y=x-(x-x0)*f(x)/(f(x)-f(x0);a(k+2)=y;if abs(y-x)>wx=y;elsebreakendendif k=ndisp('迭代失败')endfunction y=f(x) % 子函数y=x*exp(x)-1;3. kxjfdd .m % xO,x 为初值，w为精度，n为最大迭代次数fun cti ona,k=kxjfdd(x0,x1,w ,n)a(1:2)=x0,x1;for k=1: nif f(x1)-f(x0)=

3、0disp('分母为零')breakendy=x1-(x1-x0)*f(x1)/(f(x1)-f(x0);a(k+2)=y;if abs(y-x1)>wx0=x1;x1=y;elsebreakendendif k=ndisp(' 迭代失败')endfun cti ony=f(x)y二x*exp(x)-1;%子函数三、选做题：修改单点弦截法程序dxjf.m，并在计算机上运行function x2=dxjf(x0,x1,w,N) k=1;while k=Nif g(x1)=0x2='分母为零endx2=x1-f(x1)/g(x1);if abs(x2-

4、x1)vw breakendk=k+1;x仁 x2;endif k=Nx2='迭代失败endfunction y=f(x)% 子函数y=4*cos(x)-exp(x);function y=g(x)% 子函数y=(f(x) -f(x0)"(x -x0);、1、functiony=f(x, n)y=0;for k=0:n y=y+(xAk)/(jieche ng(k); endfun ctio ny仁jieche ng(n)y1=1; if n=0y1=1; else for k=1:ny仁 y1*k;endend>> f(1,10) ans =2.71828180

5、1146385>> f(2,10) ans =7.388994708994708>> f(1,100) ans =2.718281828459046>> f(0,10)ans =12、fun ctio ny=f(x ,n)y=0;for k=1:ny=y+(-1)A(k-1)*xA(2*k-1)/(jieche ng(2*k-1);endfun ctio ny1=jieche ng(a)y1=1;if a=0y1=1; else for k=1:ay仁 y1*k; end end>> f(1,10) ans =0.841470984807897&

6、gt;> f(2,10) ans =0.909297426825641 >> f(pi/2,100) ans =1.0000000000000001、>> A=1 2;2 1;>> E,k=jzhexp(A)E =10.226708182179552 9.8588287410081099.858828741008109 10.226708182179552 k =225>> A=1 0;0 1;>> E,k=jzhexp(A)E =2.71828182845904600 2.718281828459046k =1792、>

7、> a,k=xjfdd(0,1,1e-4,100)a =Colu mns 1 through 30 1.000000000000000 0.367879441171442Colu mns 4 through 60.6062435350855970.692200627555347 0.500473500563637Columns 7 through 90.545395785975027 0.579612335503379 0.560115461361089Colu mns 10 through 120.571143115080177 0.564879347391049 0.56842872

8、5029061Colu mns 13 through 150.566414733146883 0.567556637328284 0.566908911921495Colu mns 16 through 180.567276232175570 0.567067898390788 0.567186050099357Colu mn 190.56711904005721517>> a,k=xjfdd(0.5,1,1e-4,100)a =Colu mns 1 through 30.500000000000000 1.000000000000000Colu mns 4 through 60.

9、568322402555088 0.567076908522704k =0.5463692378777420.56714702931860843、>> a,k=kxjfdd(0.5,1,1e-4,100)a =Colu mns 1 through 30.500000000000000 1.000000000000000Colu mns 4 through 60.560794677581843 0.567252363024186Colu mn 70.567143290359005k =0.5463692378777420.567142721988679>> a,k=kxjfdd(0.5,6,1e-4,100) a =Colu mns 1 through 30.500000000000000Colu mns 4 thr