pub fn linear_est( x: &VectorF64, c: &VectorF64, cov: &MatrixF64 ) -> Result<(f64, f64), Value>
Returns (y, y_err).
(y, y_err)