df(r) / dr = - (f(r) / T)*rho*(V*r/R)^2 / r