p3 function

simulate Poisson point process