get_velocities_plummer

cosmic.sample.cmc.elson.get_velocities_plummer(r, r_max_cluster)

The correct way to generate velocities

this is a special case for gamma=4, which is the plummer sphere (and has an analytic distribution function)

returns (vr,vt) with same length as r