Parker wind equation (source code)

= Parker wind equation
{c}

For radial speed $u$, sound speed $c_s$, and central mass $M$, steady spherical mass and momentum conservation give
$$
\left(u-\frac{c_s^2}{u}\right)\frac{du}{dr}
=\frac{2c_s^2}{r}-\frac{GM}{r^2}.
$$
At a regular sonic point, $u=c_s$ and $r=GM/(2c_s^2)$ simultaneously.