Functions for Bayesian unfolding
Example 6.4.1 from arXiv:1201.4612v4