This function generates a realisation
of the homogeneous Poisson process in three dimensions,
with intensity lambda
(points per unit volume).
The realisation is generated inside the three-dimensional region
domain
which currently must be a rectangular box (object of
class "box3"
).