def METHOD_NAME(self, r, VP):
    r1 = r[0,:]
    r2 = r[1,:]
    pot = self.en_pot(r1, r2) + self.ee_pot(r1, r2)
    psi_val = self.psi_internal(r1, r2, VP)
    lapl = self.lap(r1, r2, VP)
    h = -0.5*lapl/psi_val + pot
    return h