Just to let others know: this method appears to be from "Numerical Recipes in C", which is now available, full-form, online. See c7-2.pdf.

this also features recipes for other distributions, plus lots of other good stuff ....

- j