code: def create_point_set(order, nsd):
      h = Rational(1, order)
      set = []
  
      if nsd == 1:
          for i in range(0, order + 1):
              x = i * h
              if x <= 1:
                  set.append((x, y))
  
      if nsd == 2:
          for i in range(0, order + 1):
              x = i * h
              for j in range(0, order + 1):
                  y = j * h
                  if x + y <= 1:
                      set.append((x, y))
  
      if nsd == 3:
          for i in range(0, order + 1):
              x = i * h
              for j in range(0, order + 1):
                  y = j * h
                  for k in range(0, order + 1):
                      z = k * h
                      if x + y + z <= 1:
                          set.append((x, y, z))
  
      return set
  
file_path: sympy/examples/advanced/fem.py
imports: from sympy import symbols, Symbol, factorial, Rational, zeros, eye, integrate, diff, pprint, reduced
label_desc: Use image node_id 2 for calling a global function with example usage: create_point_set(order, nsd) and returns: set