function r=rayonspectral (A,X,k); 

% r= rayonspectral(A,X,k); 
% Calcule (de manière approchée) le rayon spectral
% d'une matrice A (réelle ou complexe) de taille (N,N), diagonalisable et
% présentant une unique valeur propre de module maximal, à partir de
% l'algorithme itératif initié au vecteur X et itéré k fois.

x=X; 
for i=1:k
    y=A*x;
    x=y/norm(y); 
end
r=norm(y); 