function [valpr,vectpr]=vectprop(A)
// renvoie les vecteurs propres de A
  // classés par valeurs propres corresp décroissantes
  printf("calcul en cours ...\n");
  [Diag,Base]=bdiag(A);
  // diagonalisation
  if Base(:,1)==-abs(Base(:,1)) then Base=-1*Base; end;
  // afin que les coef du 1er vect prop ne soient pas tous < 0.
  [valpr,k]=sort(diag(Diag));
  // classement décroissant des valeurs propres
  vectpr = Base(:,k);
  printf("calcul termine ...\n");
endfunction