《数学计算方法实验报告》由会员分享,可在线阅读,更多相关《数学计算方法实验报告(5页珍藏版)》请在金锄头文库上搜索。
1、数学计算方法实验报告习题二2.估计用二分法求方程f(x)=x3+4x2-10=0 在区间 1,2内根的近似值,为使方程不超过10 时所需的二分次数。f(xk) 程序过程:function two (tolerance) a=1;b=2;counter=0; while (abs(b-a)tolerance) c=(a+b)/2; fa=a3+4*a2-10; fb=b3+4*b2-10; fc=c3+4*c2-10; if (fa=0|fb=0) disp(counter); elseif (fa*fcg) a=b; fa=a3+4*a2-10; ga=3*a2+8*a; b=a-fa/ga;
2、 k=k+1; end format long; disp(a); disp(k); 实验结果:程序结果计算结果8.用弦割法求方程f(x)=x3-3x2-x+9=0 在区间 -2,-1内的一个实根近似值 xk,|f(xk)|10-5. 程序过程:function eight (t) a=-2; b=-1; fa=a3-3*a2-a+9; fb=b3-3*b2-b+9; c=b-fb*(b-a)/(fb-fa); k=1; while(abs(c-b)t) a=b; b=c; fa=a3-3*a2-a+9; fb=b3-3*b2-b+9; c=b-fb*(b-a)/(fb-fa); k=k+1;
3、 end format long; disp(k); disp(b); 实验结果:计算结果9.用艾特肯算法求方程x3+4x2-10=0 在区间 1,2内的根的近似值(取 x0=1.5,g(x)= 410x,精确到|xk+1-xk|10-5,并与 2,3,6 结果比较。程序过程:function nine (tolerance) x0=1.5; y0=sqrt(10/(x0+4); z0=sqrt(10/(y0+4); x1=x0-(y0-x0)2)/(z0-2*y0+x0); k=0; while (abs(x1-x0)tolerance) x0=x1; y0=sqrt(10/(x0+4); z0=sqrt(10/(y0+4); x1=x0-(y0-x0)2)/(z0-2*y0+x0); k=k+1; end disp(x1); disp(k); 实验结果:计算结果