Saturating autocatalytic switch (source code)

= Saturating autocatalytic switch

For
$$
\dot g=s+k\frac{g^2}{1+g^2}-g,
\qquad k>2,
$$
the zero-input stable states are $0$ and $(k+\sqrt{k^2-4})/2$. The switching threshold is the <saddle-node bifurcation> determined by $f(g_c,s_c)=\partial_gf(g_c,s_c)=0$ on the low-concentration branch.