Quantum field theoretical methods are discussed to evaluate the spatial distribution in many-body systems in equilibrium. The pair and triplet distribution functions are explicitly evaluated for noninteracting cases. A general method to evaluate the pair distribution function in the chain and ring diagram approximations for interacting systems is presented in detail. The result is expressed in terms of the eigenvalues of the propagators characteristic of the diagrams. This result generalized the Gell-Mann and Brueckner formula for the equation of state, the derivation of which is also given. The eigenvalues are evaluated for Boltzmann, Bose and Fermi statistics. Application of the result to a hard-sphere Bose gas is made, and the effect of condensation is studied. It is found that the phonon velocity depends on the number of particles in the lowest state. The distribution functions of a degenerate Fermi system are also calculated.