function a =secante(fun,x0,x1,tol,maxiter)
% Aproxima por el método de la secante una raiz de la ecuacion fun(x)=0
%cercana a x0, tomando como criterio de parada abs(fun(x))<tol o la cota sobre
%el numero de iteraciones dada  por maxiter. 
%
% Variables de entrada:
%     fun: funcion a calcular la raiz, se introduce en modo simbolico 'fun'   
%     x0, x1: estimaciones iniciales para el proceso de iteración
%     tol: tolerancia en error absoluto para la raiz
%     maxiter: maximo numero de iteraciones permitidas
%
% Variables de salida:
%     a: valor aproximado de la raiz
fprintf(1, 'Metodo de la secante \n');
f0=subs(fun,x0); 
f1=subs(fun,x1);
iter=1;
while(abs(f1)>tol) & (iter<maxiter)
	a = x1-f1*((x1-x0)/(f1-f0));		        	% formula de iteracion
   f0=f1; f1=subs(fun,a); 		%Actualiza f0 y f1 
   fprintf(1, 'iter= %i, a= %x0,f= %e  \n', iter,a ,f1)
	iter = iter + 1;	                           	% Cuenta los pasos
	x0=x1;	x1=a;			% actualiza x
end
% Salida