For this choice, we obtain the solution \(x(t)\)  and \(y(t)\)  in the form [1]}
\(x(t)=\sin {(\omega t+\delta )}{\Bigg [I_{1}-q(b_1(t)+b_2(t))\Bigg ]}^\frac{-1}{2m+1},\nonumber \\y(t)=\cos {(\omega t+\delta )}{\Bigg [I_{1}-q(b_1(t)+b_2(t))\Bigg ]}^\frac{-1}{2m+1},\) 
