p3_var_rate function

simulate variable rate Poisson point process