 Quadratic accretion in a central Keplerian potential

In a fixed central Newtonian gravitational potential, pressureless free fall has speed $u=a{R}^{-1/2}$, and randomly moving stars have typical gas-relative speed $V=b{R}^{-1/2}$. A steady, undepleted spherical mass flux $A$ gives $\rho =A{R}^{-3/2}/(4\pi a)$. The Bondi--Hoyle--Lyttleton accretion rate therefore becomes $\dot{M}=\lambda M^2$ with $\lambda =G^{2}A/(a{b}^3)$, independent of radius. Integrating $d(1/M)/dt=-\lambda$ gives the displayed solution before $\lambda M_{0}t=1$. Its finite-time divergence signals failure of the fixed potential, undepleted gas supply and negligible feedback assumptions; it is not a physically unlimited mass supply.