pca: compute new eigenvectors (but it does not converge)
pca: rewrite and fix computation of covariance matrix