We report for the first time a formulation for finding the path integral of many boson systems using hyperspherical harmonics. The accuracy of the method is established by comparison with analytical results of simple systems. The partition function, and hence several thermodynamic properties are also calculated using interactions that are relevant to Bose-Einstein condensation. (C) 2005 Elsevier B.V. All rights reserved.