They ask me to find the analytic function $f(z)=u(x,y)+iv(x,y)$ if we know that $$u(x,y)=\frac{e^x( (x^2-y^2)\cos y+ 2xy\sin y)}{(x^2+y^2)^2}.$$
Am I supposed to find $v(x,y)$ satisfying the Cauchy-Riemann conditions in cartesian coords.? Isn't it too cumbersome to integrate and differentiate $u$? I don't think polars make it easier since there'd be $\sin(r\sin\theta)$ type functions and a quotient...
It's cumbersome, but not too cumbersome. On the other hand, it is a mistake to ask for the analytical function whose real part is the given function; there are infinitely many of them.
Here's a way of getting one of them. Take\begin{align}v(x,y)&=\int_0^y\frac{\partial u}{\partial x}(x,t)\,\mathrm dt\\&=\frac{e^x\left(\left(x^2-y^2\right) \sin (y)-2 x y \cos (y)\right)}{\left(x^2+y^2\right)^2}.\end{align}It is then clear that $\frac{\partial u}{\partial x}=\frac{\partial v}{\partial y}$. And a direct computation shows that $\frac{\partial u}{\partial y}=-\frac{\partial v}{\partial x}$. So, you can take\begin{align}f(x+yi)&=u(x,y)+v(x,y)i\\&=\frac{e^x\left(\left(x^2-y^2\right)\cos(y)+2xy\sin(y)\right)}{\left(x^2+y^2\right)^2}+\frac{e^x\left(\left(x^2-y^2\right)\sin(y)-2xy\cos(y)\right)}{\left(x^2+y^2\right)^2}i\\&=\frac{e^x((x^2-y^2)e^{yi}-2xyie^{yi})}{(x^2+y^2)^2}\\&=\frac{e^{x+yi}(x-yi)^2}{(x+yi)^2(x-yi)^2}\\&=\frac{e^{x+yi}}{(x+yi)^2.}\end{align}In other words, you can take $f(z)=\dfrac{e^z}{z^2}$.