function  x = biseccion(fun,a,b,tol)
% Aproxima por el método  de la bisección una raíz de la ecuación fun(x)=x
fprintf(1, 'Método de la bisección\n');
fprintf(1,'\n'); 

n=1;

u=subs(fun,a);
v=subs(fun,b);
if sign(u)==sign(v)
disp('ERROR:la funcion debe cambiar signo en a,b');
break;
end;

while ((b-a)*0.5>tol)
c=(b+a)/2; w=subs(fun,c);
fprintf(1, 'n= %i, c= %f, f(c)= %e  \n',n, c,w);
if sign(u)==sign(w)
	a = c; u=w;
else
	b=c; v=w;
end
	n=n+1;
end;
x=c