function Xs = BisectionRoot_Schultz(Fun,a,b,TolMax)
if(nargin<4)
    tol = .000001;
else
    tol = TolMax;
end
error = 100;
x = (b-a)/2 + a;
count = 0;

while((error > tol) && (count < 1000))
    count = count + 1;
    value = Fun(x);
    if(Fun(a)*value < 0)
        b = x;
    else
        a = x;
    end
    plot(x,value,'rx-')
    x = (b-a)/2 + a;
    error = abs(value);
end
if count < 1000
    disp('Bisection Method-')
    disp('Tolerance Reached')
    disp('Error:')
    disp(error)
    disp('Solution:')
    disp(x)
else
    disp('Timeout!!')
end
Xs = x;
return;
end