数值分析程序纪要(MATLAB)

一、求方程的根

    (1)二分法

数值分析程序纪要(MATLAB)

---

function xc=bisect(f,a,b,tol)

if sign(f(a))*sign(f(b)) ＞＝0

    error('f(a)f(b)<0 not satisfied')

end

fa=f(a);

fb=f(b);

while (b-a)/2＞tol

    c=(a+b)/2;

    fc=f(c);

    if fc==0

    break

end

if sign(fc)*sign(fa)<0

  b=c;fb=fc;

else

  a=c;fa=fc;

    end

end

xc=(a+b)/2;

定义函数f=@(x) 表达式，

xc=bisect(f，a，b，tool)

＊注记: