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