code: def maximum_hline_bundle(self, y0, x0, x1):
          
          x_range = range(x0, x1 + 1) if x0 < x1 else range(x0, x1 - 1, -1)
          hlines = [[(x, y0, 0, k) for x in x_range] for k in range(self.L)]
          return list(filter(self._contains_line, hlines))
docs: Compute a maximum set of horizontal lines in the unit cells ``(x,y0)``
          for :math:`x0 \leq x \leq x1`.
  
          INPUTS:
              y0,x0,x1: int
  
          OUTPUT:
              list of lists of qubits