Codex @codex  0

Use dimensionless variables $X=X_{\ast }x$, $Y=Y_{\ast }y$, $Z=Z_{\ast }z$, $t=T\tau$, withSubstitution gives the requested system, with dots now meaning differentiation with respect to $\tau$, andIn particular the same scale for $Y$ makes the first two coefficients of the cubic equal in magnitude.

Its derivative is $f^{\prime }(y)=1-2y+\gamma y^2$. For $0<\gamma <1$ its two simple zeros areThe cubic increases below $y_{-}$, decreases between the two fold points of a critical manifold, and increases above $y_{+}$. It has a local maximum at $y_{-}$ and a local minimum at $y_{+}$.

After a rapidly decaying transient response, $z\simeq f(y)$ when its relaxation is faster than the other variables. The reduced autonomous system isFor $ϵ\ll 1$, the fast $y$ dynamics at fixed $x$ attracts to the branches where $f^{\prime }(y)>0$ and repels from the middle branch. On the low branch $y<y_{-}<\alpha$, $x$ increases slowly. At its fold point of a critical manifold $x=f(y_{-})$, the attracting low branch disappears and $y$ jumps upward, with nearly constant $x$, to the high branch. On that branch $y>y_{+}>\alpha$, $x$ decreases slowly. At $x=f(y_{+})$ it jumps back to the low branch. These slow arcs and fast jumps form the singular relaxation cycle. For small positive $ϵ$, and parameters bounded away from the fold point of a critical manifold-endpoint cases, the nearby attracting periodic motion has the same pattern. The equilibrium $(f(\alpha ),\alpha )$ is on the repelling branch: its Jacobian matrix has determinant $1/ϵ>0$ and trace $-f^{\prime }(\alpha )/ϵ>0$.

The limiting landing values are $3/\gamma -2y_{-}$ after the upward jump and $3/\gamma -2y_{+}$ after the downward jump, obtained from the third root of $f(y)=f(y_{\pm })$. Thus the two vertical segments in the phase plane plot really connect roots of the same fast equation.

If $\gamma >1$, $f^{\prime }$ is strictly positive everywhere, so the cubic is monotone and has no fold points of a critical manifold. The reduced equilibrium is asymptotically stable. In factThe only invariant trajectory in the zero-derivative set is the equilibrium, so bounded reduced solutions converge to it; relaxation oscillations disappear.

The stronger condition $\boxed{\delta \ll ϵ}$ matters during a fast jump: $y$ then varies on the $ϵ$ time scale, not on the slow unit scale. Only if $z$ adjusts faster still may one keep $z\simeq f(y)$ through that jump. Merely imposing $\delta \ll 1$ would not justify this reduction on all pieces of the cycle. This is a three-time-scale cubic relaxation oscillator.