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