A Monte Carlo algorithm was constructed for calculating statistical physics
quantities characterizing an ideal bosonic lattice gas. To demonstrate the
possibilities of the algorithm the Bose-Einstein condensation was simulate
d within a three-dimensional isotropic harmonic oscillator where, in the fr
ame of a canonical ensemble, average occupancy of the ground state and spec
ific heat was calculated versus temperature and number of bosons. The algor
ithm can be further applied for studying both the static properties of idea
l bosons within other trapping potentials and the relaxation of the system
to the condensate. (C) 1999 Elsevier Science B.V. All rights reserved.