code: def gaussian_points(loc=(0, 0), scale=(10, 10), n=100):
      
  
      arr = np.random.normal(loc, scale, (n, 2))
      return gpd.GeoSeries([shapely.geometry.Point(x, y) for (x, y) in arr])
docs: Generates and returns `n` normally distributed points centered at `loc` with `scale` x and y directionality.