import spharpy
coords = spharpy.samplings.equidistant_cuboid(3)
spharpy.plot.scatter(coords)