import spharpy
coords = spharpy.samplings.gaussian(n_max=5)
spharpy.plot.scatter(coords)