This function generates a systematic random pattern of points
  in the window win. The pattern consists of a rectangular grid
  of points with a random common displacement.  The grid spacing is determined by the distances dx, dy
  if they are present. If they are absent, then the grid spacing is
  determined so that there will be nx columns and ny rows
  of grid points in the bounding rectangle of win.
  The grid is then given a random displacement (the common displacement
  of the grid points is a uniformly distributed random vector in the
  tile of dimensions dx, dy).
  Some of the resulting grid points may lie outside the window win:
  if they do, they are deleted.
  The result is a point pattern inside the window win.
  This function is useful in creating dummy points for quadrature
  schemes (see quadscheme) as well as in simulating
  random point patterns.