  Algorithms to determine transition probabilities in Monte Carlo simulations
are tested using a system of classical particles with effective interactions
which reproduce Bose-Einstein statistics. The system is appropriate for testing
different Monte Carlo simulation methods in out-of-equilibrium situations since
non equivalent results are produced. We compare mobility numerical results
obtained with transition probabilities derived from Glauber and Metropolis
algorithms. Then, we compare these with a recent method, the interpolation
algorithm, appropriate for non-equilibrium systems in homogeneous substrata and
without phase transitions. The results of mobility obtained from the
interpolation algorithm are qualitatively verified with molecular dynamics
simulations for low concentrations.
