function [n,X]=iteration_secante(xm1,x0,nmax,epsilon,fcn)


%	iteration_secante : calcul d'une racine la méthode de la sécante
%
% *********************************************************
%
%	[n,X]=iteration_secante(x0,xm1,nmax,epsilon,fcn) renvoie
%  les itérés de la méthode de la sécante
%
%	variables d'entrées : 
%		* nmax : nombre maximal d'itérations
%		* xm1 et x0 : les deux premiers termes (x_{-1} et x_0)
%     * epsilon : précision souhaitée
%     * fcn est une chaîne de caractère représentant la fonction 
%      (de type inline, builtin ou par fichier M-file);
%
%	variables de sortie :
%     * n est l'indice correspondant à x_n avec f(x_n)=0 ou |x_n-x_{n-1}|<epsilon 
%                                         ou n>=nmax
%     * X est une vecteur qui contient les valeurs de x_0 à x_n
%
% 
% 
% ************ Fonctions auxiliaires utilisées ************
%
%       aucune
%
% *********************************************************
%




% Contrôles d'entrée
qs=exist(fcn);
if (qs~=2) & (qs~=5) & (qs~=1)
   error('fcn doit être le nom d''une fonction (built-in, M-file ou inline)');
end

% corps d'algorithme
n=0;
test=(feval(fcn,x0)~=0);
x=x0;
xmm1=xm1;
fxmm1=feval(fcn,xm1);
fx=feval(fcn,x);
X=x0;
while(test)
   xnp1=x-((x-xmm1)/(fx-fxmm1))*fx;
   n=n+1;
   fxnp1=feval(fcn,xnp1);
   test=~((fxnp1==0)|(n>=nmax)|(abs(x-xmm1)<=epsilon));
   xmm1=x;
   fxmm1=fx;
   x=xnp1;
   fx=feval(fcn,x);
   X=[X x];
end   
if (n>=nmax)
   disp('attention, n>nmax');
end


