The thermodynamic and statistical properties of square-well hexamers, octamers, and hexadecamers have been studied using the Born-Green-Yvon method. We have obtained concentration- and temperature-dependent site-site correlation functions for each system and have used these results to predict the equation of state and the phase diagram. Where possible we have compared our results with those obtained using Monte Carlo simulations and other theoretical methods. (C) 2003 American Institute of Physics.