The nonequilibrium behavior of neutral atoms confined in a harmonic potential is simulated by a Monte Carlo method. Atoms obey Maxwell-Boltzmann statistics. Initially, the probability distribution function (pdf) for each of the three energy components, Ex, Ey and Ez, is taken to be a truncated...