f=inline('816*x^3-3835*x^2+6000*x-3125');

a=1;
b=2;

format long;
k=0;
while abs(b-a) > eps*abs(b)
    m=(a+b)/2;
    if sign(f(m))==sign(f(b))
        b=m;
    else
        a=m;
    end
    k=k+1;
    disp([k a b]);
end