get_velocities

cosmic.sample.cmc.elson.get_velocities(r, r_max_cluster, gamma)

The correct way to generate velocities: generate the distribution function from rho, then use rejection sampling.

returns (vr,vt) with same length as r