Python implementation of Frank van den Bosch's N-body simulation IC generation code, which uses standard inverse transform sampling to sample particle locations and then samples the NFW distribution function of Widrow (2000) using rejection sampling as described in Kuijken & Dubinski (1994).
