def METHOD_NAME(self):
    td = tdscf.TDA(mf_lda).run(nstates=nstates)
    tdg = td.nuc_grad_method()
    g1 = tdg.kernel(td.xy[2])
    self.assertAlmostEqual(g1[0,2], -0.7944872119457362, 6)
    td_solver = td.as_scanner()
    e1 = td_solver(pmol.set_geom_('H 0 0 1.805; F 0 0 0', unit='B'))
    e2 = td_solver(pmol.set_geom_('H 0 0 1.803; F 0 0 0', unit='B'))
    self.assertAlmostEqual((e1[2]-e2[2])/.002, g1[0,2], 4)