docs: def get_polygon_constraints(self,
                                  range_polygones=range(3, 5),
                                  print_out=False):
          """
          :param range_polygones: list of numbers of polygones to test.
          
          :return A, b: the constraints on the theta-vector of the form A*theta = b
          """
          rows_A = []
          rows_b = []
          for m in range_polygones:
              if (print_out):
                  print('checking {}-polygones'.format(m))
              polygons = self.get_convex_polygons(m)
              row_A, row_b = self.get_polygon_constraints_m(polygons, print_out)
              rows_A.append(row_A)
              rows_b.append(row_b)
          self.A = np.vstack(rows_A)
          self.b = np.hstack(rows_b)
          return self.A, self.b