function [VP,Niter,errY]=AppVPdom1bis(A,X,Y,epsilon,k); 

% [VP,Niter,errY] = AppVPdom1bis(A,X,Y,epsilon,k); 
%Calcule le vecteur propre approchant un vecteur propre associé au rayon
%spectral r d'une matrice A de taille (N,N)) (correspondant à une valeur
%propre réelle). L'algorithme est initié au vecteur X (de taille (N,1)). Le
%nombre d'itérations autorisées est k. La sortie Niter correspond au nombre
%d'itérations nécessaires avant que : 
% min(norm(VP(p+1)-VP(p)),norm(VP(p+1)+VP(p))) <=epsilon.
%La liste errY figure la liste des erreurs :
% min(norm(VP(k)-Y,inf),norm(VP(k)+Y,inf)), k=1,...,Niter. 

VP=X/norm(X,inf); %changement de norme 
err=2*epsilon; 
Niter=0; 
errY=[]; %ligne ajoutee (cf. AppVPdom)
while (err>epsilon) && (Niter<k) 
    VPold=VP; 
    y=A*VPold;
    r=norm(y,inf); %changement de norme 
    VP=y/r; 
    err=min(norm(VP-VPold,inf),norm(VP+VPold,inf)); %changement de norme 
    errY=[errY,min(norm(VP-Y,inf),norm(VP+Y,inf))]; %ligne ajoutee (cf. AppVPdom)
    %+ changement de norme  
    Niter=Niter+1; 
end 