monte carlo simulation for roulette buph