Abstract
Numerical methods are demonstrated for the generation of the radial coefficients in the angular expansion of the pair distribution function in the mean spherical approximation (MSA). These techniques are applied to the MSA closure for mixtures of point dipoles embedded in hard spheres of additive diameter.