Abstract.We prove absolute convergence of the multi-body correlation functions as a power series in the density uniformly in their arguments. This is done by working in the context of the cluster expansion in the canonical ensemble and by expressing the correlation functions as the derivative of the logarithm of an appropriately extended partition function. In the thermodynamic limit, due to combinatorial cancellations, we show that the coefficients of the above series are expressed by sums over some class of two-connected graphs. Furthermore, we prove the convergence of the density expansion of the "direct correlation function" which is based on a completely different approach and it is valid only for some integral norm. Precisely, this integral norm is suitable to derive the OrnsteinZernike equation. As a further outcome, we obtain a rigorous quantification of the error in the Percus-Yevick approximation.