monte carlo simulation for roulette gphm