组号 1北 方 民 族 大 学学生上机报告课程名称 数值计算方法 上机名称 牛顿迭代法求方程的根 院系 数学与信息科学学院 专业 数学与应用数学 班级 12 级应数 2 班 第 组 组长 周海龙 同组成员 周海龙 上机日期 2014.6 成绩 组号 2第 次上机 一、 上机目的进一步了解掌握牛顿迭代法求方程的根算法的设计思路和算法流程，培养动手实践能力和分析能力。 二、 上机内容（问题描述）第一部分： 牛顿迭代法求方程的根 10（B 类） 、利用迭代法求解 的全部根，要求绝3210x对误差限小于 (先确定含根区间，然后构造迭代公式8102进行求解)组号 3三、 算法描述（流程图）开始输入 x0，e，N,1=kf(x0)=0? 是否X0-f(x0)/f(x0)=x1x1-xokx1=x0组号 4四、 上机程序function k,x,wuca,yx=newton(x0,tol)k=1;x1=x0-fun(x0)/fun1(x0);disp(x=)while abs(x1-x0)tolx0=x1;k=k+1;x1=x0-fun(x0)/fun1(x0);disp(x1)endk;x=x1;wuca=abs(x1-x0)/2;yx=fun(x);end%分程序 1：function y1=fun(x)y1=x3+3*x2-1;end组号 5%分程序 2：function y2=fun1(x)%函数 fun(x)的导数y2=3*x2+6*x;end%调用函数输入数据x0=1.0;xr=-2.5;xp=-0.5;tol=0.5*10-8;k,x,wuca,yx=newton(x0,tol)k,x,wuca,yx=newton(xr,tol)k,x,wuca,yx=newton(xp,tol)五、 运行结果x=0.54860.53240.53210.53210.5321组号 6k =6x =0.5321wuca =5.9952e-015yx =2.2204e-016x=-2.9009-2.8797-2.8794-2.8794-2.8794k =6x =-2.8794wuca =2.6645e-015组号 7yx =-3.5527e-015x=-0.6528-0.6527-0.6527k =4x =-0.6527wuca =1.0850e-009yx =-1.1102e-016