pub fn center_cov(pts: &[Point<f32>]) -> (Point<f32>, Matrix<f32>)
Computes the center and the covariance matrix of a set of points.