Computing a Jacobian matrix