rgtp function

Simulation of generalized Thomas process