Covariance Matrix computation