Simulation for gamma