from preliz import ZeroInflatedBinomial, style
style.use('preliz-doc')
ns = [10, 20]
ps = [0.5, 0.7]
psis = [0.7, 0.4]
for psi, n, p in zip(ns, ps, psis):
    ZeroInflatedBinomial(psi, n, p).plot_pdf(support=(0,25))