from preliz import NegativeBinomial, style
style.use('preliz-doc')
mus = [1, 2, 8]
alphas = [0.9, 2, 4]
for mu, alpha in zip(mus, alphas):
    NegativeBinomial(mu, alpha).plot_pdf(support=(0, 20))