Psim function

Simulate Binomial Spatial Point Process