Add a function to test that the analytic Jacobian is actually the one used by the root finder