troot.m

数值分析各种公式的解法

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