A single shell, identical objects, no avoidance, no disposal. The cascade equation dN/dt = L - N/τ + kN² has a stable and an unstable fixed point; they merge at the critical launch rate Lc = 1/(4kτ²). The test integrates the ODE numerically and finds the runaway threshold by bisection.
N(t) in units of Nc = 1/(2kτ) against t/τ. Curves that cross the unstable fixed point N+ run away.
Phase line dN/dt against N for the chosen launch rate: roots are the fixed points.