The phonon Monte Carlo simulation