牛顿法和割线法.docxVIP

精品文档 作业十（第五章）：1. 在区间 (0,1.5) 上分别用 二分法、牛顿法和割线法编程求下面的函数的零点，精度要求 10-10 。 f ( x)=cos 2 (2 x) x2 二分法 function [X]=bisection(fx,xa,xb,n,delta) 二分法解方程 fx 是由方程转化的关于 x 的函数，有 fx=0 。 xa 解区间上限 xb 解区间下限 %解区间人为判断输入 n 最多循环步数，防止死循环。 %delta 为允许误差 x=xa;fa=eval(fx); x=xb;fb=eval(fx); for i=1:n 。 1欢迎下载 精品文档 xc=(xa+xb)/2;x=xc;fc=eval(fx); X=[i,xc,fc]; if fc*fa0 xb=xc; else xa=xc; end if (xb-xa)delta, break , end End 二分法结果： 迭代 34 次， xc=0.5149 牛顿法 function [X]=newton(fx,e,x0,m) x=x0;k=0; F=eval(fx); if abs(F)=e 。 2欢迎下载 精品文档 X=[x F]; disp(X); return end while k=m x=x0;g=eval(diff(fx)); x1=x0-F/g; x=x1;F=eval(fx);