The Fermi contact term is closely related to the atomic hyperfine constants. It often dominates the hyperfine splittings. The quality of the wave function near the origin and the correlation effect between electrons are two main factors which affect the numerical accuracy of the Fermi contact term. It is not an easy task to compute the Fermi contact term with high precision for a general atom. In the present paper, the Schrödinger equations of the \text1s2s3s\;^\text4\textS , \text1s2s4s\;^\text4\textS and \text1s2s2p\;^\text4\textP states of lithium atom and lithium-like ions (Z = 4–10) are solved by using Rayleigh-Ritz variational method in Hylleraas coordinates. The variational energiesenergy converges to an accuracy of 10^–13. Then the Fermi contact terms for these states are calculated based on the high precision variation wave functions. In particular, the Drachman global method are adopted in order to improve the convergence of the Fermi contact term. The effect of finite nuclear mass on Fermi contact term, i.e. the first-order mass polarization coefficient is also calculated. The Fermi contact term converges to an accuracy of 10^–10, which is the most accurate result at present. Our results can be used as a reference for other theoretical methods or relevant experimental studies.