For vc=220kms−1, this is about 0.0140(r/8kpc)−2M⊙pc−3. It is the spherical model's total density; it is singular if extrapolated to the very centre and is intended only over the flat-rotation region.
For a sharp outer truncation at r∗, the exterior mass is fixed at M∗=vc2r∗/G. Integrating the gravitational work from r to the boundary and then to infinity gives
Without those extra assumptions, the measurement gives only a bound: a bound star's total speed is at least its radial component, and the sample maximum need not reach escape, so r∗≳19.9kpc within this model. An unbound interloper would invalidate even that bound. This is not adetermination of the real Milky Way halo edge.