武汉理工数值分析大作业.rar

源程序： function [x_star,k]=bisect0(fun,a,b,ep) %二分法解非线性方程,f(x)=0 %fun(x)为要求根的函数f(x),a,b初始区间的端点 %ep为精度(默认为1e-5),当(b-a)/2<ep时终止计算 %x_star为迭代成功时方程的根,k表示迭代次数 %当输出迭代次数k为0时表示此区间没有根存在 if nargin<4 ep=1e-5; end fa=feval(fun,a);fb=feval(fun,b); if fa*fb>0 x_star=[fa,fb]; k=0; return; end k=1; while abs(b-a)/2>ep x=(a+b)/2;fx=feval(fun,x); if fx*fa<0 b=x;fb=fx; else a=x;fa=fx; end k=k+1; end x_star=(a+b)/2; 输入 fun=inline('x^3-x-1'); [x_star,k]=bisect1(fun,1,1.5,0.005)