pub fn covar(J: &MatrixF64, epsrel: f64, covar: &mut MatrixF64) -> Value
Compute the covariance matrix cov = inv (J^T J) by QRP^T decomposition of J