I'm trying to understand the NFW profile and how it causes a flat rotation curve. $$ \rho(r) = \frac{\rho_0}{\frac{r}{r_s}\left(1+\frac{r}{r_s}\right)^2} $$
I think I've got why it causes a flat rotation curve. The density seems to be inversely proportional to $r^3$ when $r >> r_s$ leading to a linearly increasing mass in proportion to radius since volume is proportional to $r^3$.
However I realized I can't really prove my statement above as I don't really understand what $r_s$ and $\rho_0$ are. Thus I was hoping to see a numerical example to get my head around the relationships.
While looking for one I came across this answer by Kyle Oman applying the NFW profile to calculate the dark matter density in the Solar System and was wondering where he got his $r_s$ and $\rho_0$ values or how one calculates them themselves.