code: def create_chebyshev_samples(order, dim=1):
      
      x_data = .5*numpy.cos(numpy.arange(order, 0, -1)*numpy.pi/(order+1)) + .5
      x_data = chaospy.quad.combine([x_data]*dim)
      return x_data.T
docs: Chebyshev sampling function.
  
      Args:
          order (int):
              The number of samples to create along each axis.
          dim (int):
              The number of dimensions to create samples for.
  
      Returns:
          samples following Chebyshev sampling scheme mapped to the
          ``[0, 1]^dim`` hyper-cube and ``shape == (dim, order)``.