<SYSTEM_TASK:>
Solve the following problem using Python, implementing the functions described below, one line at a time
<END_TASK>
<USER_TASK:>
Description:
 def gaussian_points(loc=(0, 0), scale=(10, 10), n=100):
 """ Generates and returns `n` normally distributed points centered at `loc` with `scale` x and y directionality. """


    arr = np.random.normal(loc, scale, (n, 2))
    return gpd.GeoSeries([shapely.geometry.Point(x, y) for (x, y) in arr])