pub fn covariance_py( x: &Bound<'_, PyArray1<f64>>, y: &Bound<'_, PyArray1<f64>>, ddof: usize, ) -> PyResult<f64>
Calculate covariance