pub fn covariance_matrix(x: &Matrix<f64>, axis: Axis) -> Matrix<f64>
Calculates the covariance matrix of the input data. Assumes input x is (n_samples, n_features).
x