function [root1,ea1,iter1]=bisect(func,xl,xu,es,maxit,varargin)

while (1)

xrold = xr

    xr = (xl + xu)/2;

    iter1 = iter1 + 1;

    

    if xr ~= 0,ea1 = abs((xr - xrold)/xr)*100;

    end

    

    test = func(xl,varargin{:})*func(xr,varargin{:});

    if test < 0

        xu = xr;

    elseif test > 0

        xl = xr;

    else

    ea1 = 0;

    end

    

    if ea1 <= es | iter1 >= maxit,break,end