We're given the following dynamical system:
$$ \begin{aligned} \dot x &= -x + y + x (x^2 + y^2)\\ \dot y &= -y -2x + y (x^2 + y^2) \end{aligned} $$
What's the largest constant $r_0$ such that the circle $x^2 + y^2 < r_0^2$ lies in the origin's basin of attraction?
So far, with relatively easy algebra, I've got:
$$ \begin{aligned} \dot r &= \frac{r}{2}(-2-\sin(2 \phi)+2r^2) \\ \dot \phi &= -(1+\cos^2(\phi)) \end{aligned} $$
Which immediately shows $r_0 \geq \sqrt{\frac12}$. How to show that there is no better bound?
Best Answer
(This is only a "numerical" answer to the question)
A basin of attraction can be found numerically using reverse-time (or backward) integration. If we choose some initial point and integrate backwards for a sufficient time, we can obtain the set of the system states that it had before it reached the selected initial point. So, if we choose sufficiently many initial points that are contained in the basin of attraction (i.e. sufficiently close to the attracting steady state) and integrate backwards, we can get some idea of how the basin of attraction looks like.
The following picture demonstrates the aprroximation of the basin of attraction of the system: The blue curves fill the basin of attraction. The red circle is the largest circle that fits in the basin of attraction. I have choosen 36 initial points on the circle of radius 0.05. Here is the Matlab code: