Function finite_diff_jacobian 

pub fn finite_diff_jacobian<S: Scalar>(
    r: &dyn Fn(&[S], &mut [S]),
    x: &[S],
    m: usize,
    jac: &mut [S],
)

Compute a finite-difference Jacobian of residual r at x. jac is row-major m*n.