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