Monte Carlo calculation