function  x = regula_falsi(fun,a,b,maxiter)
% Aproxima por el método  de la regula falsi una raíz de la ecuación fun(x)=x

fprintf(1, 'Método de la regula falsi\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;

for n=1:1:maxiter
   c=a-(u*(b-a)/(v-u)); 
   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