troot.m

来自「数值分析各种公式的解法」· M 代码 · 共 17 行

function r=troot(a,b,eps)
k=0
while abs(b-a)>eps
ab2=(a+b)/2
switch sign(f(ab2)*f(a))
    case 0
        r=(a+b)/2
    case -1
        b=(a+b)/2
    case 1
        a=(a+b)/2
end
k=k+1
end
r=ab2