Transonic spherical accretion rate in a power-law potential (source code)

= Transonic spherical accretion rate in a power-law potential
{title2=$\dot M$}

With reservoir density $\rho_\infty$ and <adiabatic sound speed> $c_\infty$, put $q=2/\beta-1/2$ and $h=1-q(\gamma-1)>0$. The selected <mass accretion rate> is $4\pi\rho_\infty(A\beta/2)^{2/\beta}c_\infty^{1-4/\beta}h^{q-1/(\gamma-1)}$. It follows by evaluating the conserved spherical flux at the <sonic point>.