pp = beta_prior(p, 2, 2)
post = posterior(pp, l)
plot_post(pp, l, post, p)