function c=covariance(echant)
//      retourne la matrice de covariance
//      de l'echantillon echant. L'entree est une matrice 
//      a n lignes (individus) et d colonnes (variables)

moyenne = mean(echant,"r");
echcent = echant - ones(size(echant,"r"),1)*moyenne;
c=(echcent'*echcent)/size(echant,"r");

function c=correlation(echant)
//      retourne la matrice de correlation
//      de l'echantillon echant. L'entree est une matrice 
//      a n lignes (individus) et d colonnes (variables)

cov=covariance(echant);
ectinv=(1)./sqrt(diag(cov));
c=diag(ectinv)*cov*diag(ectinv);







