I need to study the Hamiltonian differential system $$ \begin{align} \dot{x} &= -2ye^{-x^2}\\ \dot{y} &= 2xe^{-x^2}(1-y^2) \end{align}$$ with Hamiltonian function $$ \begin{align} H &\colon \mathbb{R}^2 \to \mathbb{R} \\ (x,y) &\mapsto e^{-x^2}(1-y^2) \end{align}$$
Thanks to the explicit independence of $H$ from the time, $H$ is conserved along the solutions. So I need to study the level sets of $H$. In $(0,0)$ there is a local maxima, so there is a neighborhood of the constant solution $0,0$ in which the other solutions are cycles. I wanted to prove that this neighborhood is in fact the region (strictly) inside the the two lines $(1-y^2)=0$. But I don't know how to prove that in such region, the hamiltonian $H$ has "cyclic" (closed curves) level sets.
Here is the Plot and the contour lines computed by Wolfram Alpha.
When $H\gt0$, $1-y^2\ge H$, or $$|y|\le\sqrt{1-H}.$$
The curve has the explicit equation $$x=\pm\sqrt{\ln{H(1-y^2)}},$$ formed of two symmetric arcs (it's an oval), and $$|x|\le\sqrt{\ln H}.$$