code: def approx_contains(self, point, atol):
          
          try:
              # Duck-typed check of type
              point = np.array(point, dtype=np.float, copy=False, ndmin=1)
          except (ValueError, TypeError):
              return False
  
          if point.size == 0:
              return True
          elif point.shape != (self.ndim,):
              return False
  
          return self.dist(point, exponent=np.inf) <= atol
docs: Return ``True`` if ``point`` is "almost" contained in this set.
  
          Parameters
          ----------
          point : `array-like` or float
              Point to be tested. Its length must be equal to `ndim`.
              In the 1d case, ``point`` can be given as a float.
          atol : float
              Maximum allowed distance in maximum norm from ``point``
              to ``self``.
  
          Examples
          --------
          >>> min_pt, max_pt = [-1, 0, 2], [-0.5, 0, 3]
          >>> rbox = IntervalProd(min_pt, max_pt)
          >>> # Numerical error
          >>> rbox.approx_contains([-1 + np.sqrt(0.5)**2, 0., 2.9], atol=0)
          False
          >>> rbox.approx_contains([-1 + np.sqrt(0.5)**2, 0., 2.9], atol=1e-9)
          True