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").