Metropolis-Hastings, or how to sample without the integral
Bayesian work hands you densities you cannot integrate or invert. Metropolis-Hastings samples them anyway, with a random walk and an acceptance ratio that cancels the normalization constant. Includes the NumPy implementation and the three ergodicity conditions.