 Solution

Work in the centre-of-mass frame and define $I={\sum }_{i}m|{\mathbf{r}}_i{|}^2$, $T=\frac{1}{2}{\sum }_{i}m|{\mathbf{v}}_i{|}^2$. Differentiation givesFor each gravitational pair, the two contributions combine toSumming pairs therefore gives the exact scalar virial theoremIn a stationary equilibrium, or a bounded long-time average with vanishing average $\ddot{I}$, $2T=-W$. With $M=Nm$, the stated gravitational radius gives $W=-G{M}^2/R_g$, henceThe rms speed itself is $[GNm/R_g{]}^{1/2}$. This is a three-dimensional mean-square speed; for an isotropic system the one-component velocity dispersion is smaller by a factor of three in its square. No assumption about a uniform spatial density was used.