from preliz import ZeroInflatedNegativeBinomial, style
style.use('preliz-doc')
psis = [0.7, 0.7]
mus = [2, 8]
alphas = [2, 4]
for psi, mu, alpha in zip(psis, mus, alphas):
    ZeroInflatedNegativeBinomial(psi, mu=mu, alpha=alpha).plot_pdf(support=(0,25))