def _hess_two_param(self, funct, p0, p1, dl=<NUM_LIT>, rts=False, **kwargs):
vals0 = self.get_values(p0)<EOL>vals1 = self.get_values(p1)<EOL>f00 = funct(**kwargs)<EOL>self.update(p0, vals0+dl)<EOL>f10 = funct(**kwargs)<EOL>self.update(p1, vals1+dl)<EOL>f11 = funct(**kwargs)<EOL>self.update(p0, vals0)<EOL>f01 = funct(**kwargs)<EOL>if rts:<EOL><INDENT>self.update(p0, vals0)<EOL>self.update(p1, vals1)<EOL><DEDENT>return (f11 - f10 - f01 + f00) / (dl**<NUM_LIT:2>)<EOL>
Hessian of `func` wrt two parameters `p0` and `p1`. (see _graddoc)