code: def create_grid_samples(order, dim=1):
      
      x_data = numpy.arange(1, order+1)/(order+1.)
      x_data = chaospy.quad.combine([x_data]*dim)
      return x_data.T
docs: Create samples from a regular grid.
  
      Args:
          order (int):
              The order of the grid. Defines the number of samples.
          dim (int):
              The number of dimensions in the grid
  
      Returns (numpy.ndarray):
          Regular grid with ``shape == (dim, order)``.