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