Vectorized the CDF inverse function.  This speeds up computation by a factor of ~90 for 500k particles
Small changes...added a seed argument for reproducibly generating random positions