\frac{d N}{d t} = n_{H} \cdot 4 π {r_{i}}^{2} \frac{d r_{i}}{d t} = Q (1 - {(r_{i} / r_{s})}^{3})