Tengo la siguiente ODA:
$$\dot x=-y(x^2+y^2), \dot y=x(x^2+y^2)$$
Quiero esbozar el retrato de fase (manualmente) y quiero encontrar el flujo $\phi_t$, la órbita $O(x_0)$ y el límite establecido $\omega(x_0)$
Empiezo tomando coordenadas polares y cambio el sistema a $\dot r=-r^3\sin\theta, \dot\theta=r^3\cos\theta$
El retrato de fase se parece al de un centro estable, ¿verdad?
¿Cómo puedo seguir encontrando el flujo de la función, es decir, la solución de la ecuación diferencial?