ddG_MC2004/compute_Jacobian works. Writing newton/solve with Levenberg-Marquardt algorithm