monte carlo simulation for roulette bohc