Sample from distributions commonly used in astronomy.
pip install astrosampler
numpy
from astrosampler import SalpeterIMF
# Sample 1000 stellar masses from the Salpeter (1955) IMF
imf = SalpeterIMF(mass_range=(0.08, 100))
masses = imf.sample(1000)
A package for sampling typical distributions in astronomy
pip install astrosampler==0.0.2