code: def gaussian_multi_polygons(points, n=10):
      
      polygons = gaussian_polygons(points, n*2)
      # Randomly stitch them together.
      polygon_pairs = [shapely.geometry.MultiPolygon(list(pair)) for pair in np.array_split(polygons.values, n)]
      return gpd.GeoSeries(polygon_pairs)
docs: Returns an array of approximately `n` `shapely.geometry.MultiPolygon` objects for an array of
      `shapely.geometry.Point` objects.