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


%	iteration_point_fixe : calcul des itérés de la méthode de point fixe
%
% *********************************************************
%
%	[n,X]=iteration_point_fixe(x0,nmax,epsilon,fcn) renvoie
%  les itérés de la méthode de point fixe
%
%	variables d'entrées : 
%		* nmax : nombre maximal d'itérations
%		* x0 : terme initial
%     * 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)-x0~=0);
x=x0;
X=x0;
xm1=x0;
while(test)
   x=feval(fcn,x);
   X=[X x];
   n=n+1;
   test=~((feval(fcn,x)-x0==0)|(n>=nmax)|(abs(x-xm1)<=epsilon));
   xm1=x;   
end
if (n>=nmax)
   disp('attention, n>nmax');
end


