  We construct analytic formulas that represent the coupling-constant-averaged
pair distribution function $\gxcav(r_s,\zeta, k_Fu)$ of a uniform electron gas
with density parameter $r_s =(9\pi/4)^{1/3}/k_F$ and relative spin polarization
$\zeta$ over the whole range $0<r_s<\infty$ and $-1<\zeta<1$, with
energetically-unimportant long range ($u\to \infty$) oscillations averaged out.
The pair distribution function $g_{xc}$ at the physical coupling constant is
then given by differentiation with respect to $r_s$. Our formulas are
constructed using {\em only} known theoretical constraints plus the correlation
energy $\ec(r_s,\zeta)$, and accurately reproduce the $g_{xc}$ of the Quantum
Monte Carlo method and of the fluctuation-dissipation theorem with the
Richardson-Ashcroft dynamical local-field factor. Our $g_{xc}$ seems to be
correct even in the high-density ($r_s\to 0$) and low-density ($r_s \to
\infty$) limits. When the spin resolution of $\ec$ into $\uu$, $\dd$, and $\ud$
contributions is known, as it is in the high- and low-density limits, our
formulas also yield the spin resolution of $g_{xc}$. We also analyze the
kinetic energy of correlation into contributions from density fluctuations of
various wavevectors.
