def METHOD_NAME(self):
    mc = mcscf.approx_hessian(mcscf.CASSCF(mf, 9, 8), auxbasis='weigend')
    mc.conv_tol = 1e-8
    mo = mf.mo_coeff.copy()
    mo[:,[15,16,17,18]] = mf.mo_coeff[:,[17,18,15,16]]
    emc = mc.mc1step(mo)[0]
    self.assertAlmostEqual(emc, -230.72211519779304, 6)
    #?self.assertAlmostEqual(emc, -230.72681659920534, 6)