Function colored_jacobian 

pub fn colored_jacobian<F, Func>(
    f: Func,
    x: ArrayView1<'_, F>,
    pattern: &SparsePattern,
    eps: F,
) -> IntegrateResult<SparseJacobian<F>>
where
    F: IntegrateFloat,
    Func: Fn(ArrayView1<'_, F>) -> IntegrateResult<Array1<F>>,

Compute sparse Jacobian using coloring