Monte Carlo simulations in the NVT ensemble have been performed for fluids of hard spheres plus attractive inverse power tails, to obtain the radial distribution function, the equation of state and the excess energy. The results obtained from the Monte Carlo simulations are compared with the predictions of the Tang-Lu [Y. Tang, B.C.-Y. Lu, Mol. Phys. 90 (1997) 215-224] integral equation theory. It is found that the agreement between theory and simulation is quite good for the radial distribution function and excellent for the equation of state and the excess energy. (c) 2007 Elsevier B.V. All rights reserved.