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