docs: def i4_sobol_generate(dim_num, n, skip=1):
      """
      i4_sobol_generate generates a Sobol dataset.
  
      Parameters:
        Input, integer dim_num, the spatial dimension.
        Input, integer N, the number of points to generate.
        Input, integer SKIP, the number of initial points to skip.
  
        Output, real R(M,N), the points.
      """
      r = np.full((n, dim_num), np.nan)
      for j in range(n):
          seed = j + skip
          r[j, 0:dim_num], next_seed = i4_sobol(dim_num, seed)
  
      return r