% Function for the bisection method.
% Finds root for given function f, left and right boundaries a and b and tolerance TOL
% f(a).*f(b)<0 is a requirement for this function to work

function c = bisection(f, a, b, TOL)
  while(b−a)/2>TOL
    c = (a+b)./2;
    if f(c)=0
      return
    end %if
  
    if f(a)f(c)<0
      b=c;
    else
      a=c;
    end %If
  end %While  
end %Function