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

ind = b-a;

while ind > tol

    xx = (a+b)/;

    if f(a)*f(xx) <

        b = xx;

    else

        a = xx;

    end

    ind = b - a;

end

xc = xx;