code: def in_hull(points, hull):
      
      from scipy.spatial import Delaunay, ConvexHull
      if isinstance(hull, ConvexHull):
          hull = hull.points
      hull = Delaunay(hull)
      return hull.find_simplex(points) >= 0
docs: Test if a list of coordinates are inside a given convex hull
  
      Parameters
      ----------
      points : array_like (N x ndims)
          The spatial coordinates of the points to check
  
      hull : scipy.spatial.ConvexHull object **OR** array_like
          Can be either a convex hull object as returned by
          ``scipy.spatial.ConvexHull`` or simply the coordinates of the points
          that define the convex hull.
  
      Returns
      -------
      result : 1D-array
          A 1D-array Boolean array of length *N* indicating whether or not the
          given points in ``points`` lies within the provided ``hull``.