using the algorithm by Y. Chung (see writeup.bib) for simulatin lower-truncated gammas