monte carlo simulation for roulette bphs