<s>[INST] converting a differential equation to polar coordinates I have the following family of autonomous systems, I'm having trouble with part b):
$$x'=x(1-\sqrt {x^2+y^2})-y-\epsilon y$$
$$y'=y(1-\sqrt {x^2+y^2})+x+\epsilon(x+x^2+y^2)$$
a) Convert the system to polar coordinates $(r(t;r_0,\theta_0,\epsilon),\theta(t;r_0,\theta_0,\epsilon))$:
$$x=r\cos\theta \qquad y=r\sin\theta \qquad r^2=x^2+y^2$$
$$x'=r'\cos\theta - r\sin\theta\theta' \quad y'=r'\sin\theta + r\cos\theta\theta'$$
$$\left ( \begin{array}\\ \cos{\theta} & -r \sin{\theta} \\ \sin{\theta} & r \cos{\theta} \end{array} \right ) \left ( \begin{array}\\ r' \\ \theta' \end{array} \right ) = \left ( \begin{array}\\ r\cos\theta(1-r)-r\sin\theta-\epsilon r \sin\theta\\r\sin\theta(1-r) + r\cos\epsilon + \epsilon(r\cos\theta+r^2) \end{array} \right )$$
I found the inverse matrix:
$$\left ( \begin{array}\\ \cos{\theta} & -r \sin{\theta} \\ \sin{\theta} & r \cos{\theta} \end{array} \right )^{-1}=\left ( \begin{array}\\ \cos{\theta} &  \sin{\theta} \\ \frac{-1}{r}\sin{\theta} & \frac{1}{r} \cos{\theta} \end{array} \right )$$
Multiplied on the left side, simplified and obtained:
$$\left ( \begin{array}\\ r' \\ \theta' \end{array} \right )=\left ( \begin{array}\\ r(1+\epsilon) \\1 + \epsilon r \end{array} \right )$$
b)See that for $\epsilon=0$ there is a periodic circular orbit with radius 1, period $2\pi$, and is attractive. In particular $(r(t;1,0,0),\theta(t;1,0,0))=(1,t)$.
For $\epsilon=0$ what I get is:
$$\left ( \begin{array}\\ r' \\ \theta' \end{array} \right )=\left ( \begin{array}\\ r \\ 1 \end{array} \right )$$
$$r'(t)=r(t) \quad \theta'(t)=1$$
Therefore $\theta (t)=t$ matches the initial value problem $\theta(t;1,0,0)=t$. However $r(t)=r'(t)$ can't match the initial value problem $r(t;1,0,0)=(1)$ because:
$$r(t)=1 \Rightarrow r'(t)=0 \Rightarrow r(t) \neq r'(t)$$
Where did I go wrong?
 [/INST] $$\left ( \begin{array}\\ \cos{\theta} & -r \sin{\theta} \\ \sin{\theta} & r \cos{\theta} \end{array} \right )^{-1}= \left ( \begin{array}\\ \cos{\theta} &  \sin{\theta} \\ \frac{-1}{r}\sin{\theta} & \frac{1}{r} \cos{\theta} \end{array} \right )$$
Multiplied on the left side, simplified and obtained:
$$\left ( \begin{array}\\ \cos{\theta} &  \sin{\theta} \\ \frac{-1}{r}\sin{\theta} & \frac{1}{r} \cos{\theta} \end{array} \right ) * \left ( \begin{array}\\ r\cos\theta(1-r)-r\sin\theta-\epsilon r \sin\theta\\r\sin\theta(1-r) + r\cos\theta + \epsilon(r\cos\theta+r^2) \end{array} \right )$$
$$\left ( \begin{array}
*r\cos^2\theta (1-r)-r\sin\theta\cos\theta-\epsilon r \sin\theta\cos\theta+r\sin^2\theta(1-r) + r\cos\theta\sin\theta + \epsilon\sin\theta (r\cos\theta+r^2) 
\\ -\sin\theta \cos\theta(1-r)+\sin^2\theta+\epsilon \sin^2 \theta+\sin\theta\cos\theta(1-r)+\cos^2\theta+\epsilon(\cos^2 \theta+r\cos\theta)
\end{array} \right )$$
Cancel out:
$$\left ( \begin{array}  +r(1-r) + \epsilon\sin\theta r^2
\\ +1+\epsilon+\epsilon r\cos\theta
\end{array} \right )$$
Doesn't equal:
$$\left ( \begin{array}\\ r' \\ \theta' \end{array} \right )=\left ( \begin{array} \\ r(1+\epsilon) \\1 + \epsilon r \end{array} \right )$$
</s>