Custom rule for differentiating through Newton solver using ForwardDiff: works for gradient, fails for hessian

I would be very interested in an example!