function [x,n]=two(a,b,e)
%x=two(f,a,b,e)
if f(a)*f(b)>=0
    return
end
n=0;
while abs(a-b)>e
    x=(a+b)/2;
    y=f(x);
    if abs(y)<e
        return
    end
    if f(a)*f(x)<0
        b=x;
    end
    if f(x)*f(b)<0
        a=x;
    end
    n=n+1;
end
n