Now, du is the extra mass density required to support the "weight" of the r to r+dr section.
Therefore du * s = u * dr * (GM/(r^2) - r*w^2)
or du/dr = (u/s) * (GM/(r^2) - r*w^2)
Now to solve this differential equation
u(r) = A * exp(-GM/(r*s) - (0.5s^-1) * (r^2) * (w^2))
with A the undetermined constant.
Set u(R0) = u0
Get A = m/s*(GM/R0^2 - R0*w^2) * exp(GM/(R0*s) + (0.5s^-1) * (R0^2) * (w^2))
EDIT: dropped a factor of 1/s
Therefore du * s = u * dr * (GM/(r^2) - r*w^2)
or du/dr = (u/s) * (GM/(r^2) - r*w^2)
Now to solve this differential equation
u(r) = A * exp(-GM/(r*s) - (0.5s^-1) * (r^2) * (w^2))
with A the undetermined constant.
Set u(R0) = u0
Get A = m/s*(GM/R0^2 - R0*w^2) * exp(GM/(R0*s) + (0.5s^-1) * (R0^2) * (w^2))
EDIT: dropped a factor of 1/s
Comment