数值计算方法 方程求根.doc

计算方法实习报告数值计算方法实验报告实验内容： 方程求根 实 验 室： 专业班级： 学 号： 姓 名： 2.用MATBAB软件，用二分法求方程f（x）=x3+4*x2-10=0在区间1,2内根的近似值，为使误差不超过10-5时所需要的二分次数。function bisection_time(tolerance)a=1;b=2;k=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)|(fc=0) disp(k); elseif(fa*fc0) b=c;k=k+1; elseif(fc*fbtolerance) x0=x1; k=k+1; a=x03+4*x02-10; b=3*x02+8*x0; x1=x0-a/b;enddisp(x1);disp(k);运行结果1.3652338.弦割法求方程f（x）=x3-3*x2-x+9=0在区间-2，-1内的一个实根近似值Xk，使|f(x) |k)x0=x1;x1=x2;a=x13-3*x12-x1+9;b=x03-3*x02-x0+9;x2=x1-a*(x1-x0)/(a-b);t=t+1;enddisp(x1);disp(t)运行结果-1.52510269.用艾特肯算法求方程f（x）=x3+4*x2+10=0在区间1,2内的根的近似值（取X0=1.5，g（x）=，精确到|Xk+1-Xk|tolerance);x1=x0-(y0-x0)2/(z0-2*y0+x0);y1=sqrt(10/(x1+4);z1=sqrt(10/(y1+4);z0=z1;y0=y1;k=k+1;enddisp(z0);disp(k);运行结果1.365231结果比较方法次数结果二分法171.36523弦割法3-1.525102牛顿迭代61.36523艾特肯11.365236