matlab root finding function
% 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