fn covariance_matrix(data: &[Vec<f64>], eps: f64) -> Vec<Vec<f64>>
Compute covariance matrix with regularization.