Run principal component analysis on the dataset X
[U, S, X] = pca(X) computes eigenvectors of the covariance matrix of X
Returns the eigenvectors U, the eigenvalues (on diagonal) in S