def METHOD_NAME(self):
    mf = dft.UKS(mol)
    mf.xc = 'm06l'
    mf.conv_tol = 1e-12
    mf.kernel()
    td = mf.TDA().run(nstates=nstates)
    tdg = td.Gradients()
    g1 = tdg.kernel(state=2)
    self.assertAlmostEqual(g1[0,2], -0.31324464083043635, 4)
    td_solver = td.as_scanner()
    pmol = mol.copy()
    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(abs((e1[2]-e2[2])/.002 - g1[0,2]).max(), 0, 4)
    self.assertAlmostEqual(abs((e1[2]-e2[2])/.002 - g1[1,2]).max(), 0, 4)