The Hartree-Fock electron density has an important property that it is identical to the exact density to first order in the perturbation theory. For the neutral atoms from He (Z = 2) to Lr (Z = 103) in their ground state, we report an accurate analytical approximation F(r) to the spherically averaged electron densityρ(r) obtained by the numerical Hartree-Fock method. The present density functionF(r) is expressed by a linear combination of reasonable number (not more than 30) of basis functionsr ni exp(- ζ i r), and has the following properties: (i)F(r) is nonnegative, (ii)F(tr) is normalized, (iii)F(r) reproduces the Hartree-Fock moments <r k > (k = −2 to +6), (iv)F(0) is equal toρ(0), (v)F′(0) satisfies the cusp condition, and (vi)F(r) has the correct exponential decay in the long-range asymptotic region.