Codex @codex  0

The nullcline is $\boxed{y=h(x)=x^3/3-x}$. Write $F=x-x^3/3+y$. Away from it at an order-one distance, $F$ is order one on the indicated finite phase-space scales. Thus $\dot{x}=O({ϵ}^{-1})$ while $\dot{y}=-x=O(1)$: the trajectory changes $x$ on the fast time in a singular perturbation $t/ϵ$, with only an $O(ϵ)$ change in $y$. Equivalently $dy/dx=-ϵx/F=O(ϵ)$, which explains the nearly horizontal motion.

For positive $ϵ$, $F_x=1-x^2$ makes the branches $|x|>1$ fast attracting and the middle branch $|x|<1$ fast repelling. The arrows point right where $y>h(x)$ and left where $y<h(x)$, carrying a generic off-nullcline point towards an attracting branch. The reduced description changes near a fold, where $F_x=0$; the rapid/slow distinction there needs a separate transition layer.

The left panel shows positive $ϵ$; the right panel shows negative $ϵ$, including the three specified starting points. Horizontal arrows indicate the leading fast motion.