pub fn sample_cov(x: &[f64], y: &[f64]) -> Result<f64, String>
Sample covariance: sum((xi-mx)(yi-my)) / (n-1).