pub fn covariance_horizontal(x: &Matrix<f64>) -> Matrix<f64>
Covariance between rows (i.e. across columns)