Пусть частица массой $m$ движется в неподвижной среде, в которой на неё действует сила вязкого трения $\vec{F}_\text{тр} = -\varepsilon m \vec{v}$, где $\vec v$ – скорость частицы. Ускорение свободного падения $\vec g$, изначально частица находилась в начале координат и имела скорость $\vec v_0$.
Рассмотрим теперь движение частицы в отсутствии силы тяжести, но теперь на неё действует постоянная по модулю вращающаяся сила $F_x = F_0 \cos(\omega t)$, $F_y = F_0 \sin(\omega t)$.
Известно, что у полученных вами уравнений есть решение вида
$$
\left(x(t),y(t)\right)=\left(A e^{at} \cos(\Omega t + \varphi_0), A e^{at} \sin(\Omega t+ \varphi_0)\right).
$$Эти функции являются решениями при произвольных значениях $\varphi_0$, в частности при $\varphi_0 = 0$ и $\varphi_0 = \pi/2$ получим $$\left(x_1(t),y_1(t)\right)=\left(A e^{at} \cos(\Omega t), A e^{at} \sin(\Omega t)\right),$$ $$\left(x_2(t),y_2(t)\right)=\left(-A e^{at} \sin(\Omega t), A e^{at} \cos(\Omega t)\right).$$ При этом существует две пары возможных значений параметров $(a, \Omega)$.
В дальнейшем будем использовать обозначения
$$A\vec{r}_{1i}=\left(Ae^{a_it} \cos(\Omega_i t), A e^{a_it} \sin(\Omega_i t)\right), \;A\vec{r}_{2i}=\left(-Ae^{a_it} \sin(\Omega_i t), A e^{a_it} \cos(\Omega_i t)\right).$$ Здесь индекс $i$ может принимать два значения $1, 2$, отвечающие двум возможным значениям $(a, \Omega)$.