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
');
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
', iter,a ,f1)
iter = iter + 1; % Cuenta los pasos
x0=x1; x1=a; % actualiza x
end
% Salida
