Abstract
A radial distribution function for hard spheres is obtained, which improves at high densities on the Percus-Yevick function. This is achieved by solving the mean spherical integral equation for a Yukawa form, adjusted to satisfy the Carnahan-Starling equation of state. The solution to the related problem of the mean spherical model approximation for a model fluid which interacts via the hard sphere potential plus the Yukawa potential is reported.