<SYSTEM_TASK:>
Solve the following problem using Python, implementing the functions described below, one line at a time
<END_TASK>
<USER_TASK:>
Description:
 def create_grid_samples(order, dim=1):
 """ 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)``. """

    x_data = numpy.arange(1, order+1)/(order+1.)
    x_data = chaospy.quad.combine([x_data]*dim)
    return x_data.T