popstock
sample_powerlaw()
This is the actual spectral index x^alpha, not the negative of the spectral index