Finite difference Jacobian

Implement the computation of the Jacobian with impl_jacobian and finite differences.