I have $\Re f(z) = e^x(xcosy -ysiny) = e^xcosy -e^xysiny$. Therefore, $$v_x= e^xxcosy+e^xcosy-e^xysiny = e^x(xcosy+cosy-ysiny)$$ $$v_y= -e^xxsiny-e^x(siny-ycosy) = e^x(-xsiny-siny+ycosy)$$
$u_x = v_y = e^x(-xsiny-siny+ycosy)$
$u_y = -v_x = -e^x(xcosy+cosy-ysiny)$
$$u = \int u_xdx = -\int e^x(xsiny+siny-ycosy)dx = ... = -e^x(xsiny-ycosy)$$ To find $u_y$, I do:
$$u_y = -e^x(xcosy+cosy-ysiny) = -e^x(xsiny-ycosy)$$ What do I do from here?
First of all, we observe that $$ \mathrm{e}^x(x\cos y-y\sin y)=\mathrm{Re}\,\big(\mathrm{e}^x(\cos y+i\sin y)\cdot (x+iy)\big)=\mathrm{Re}\,(\mathrm{e}^zz). $$ Next, use the fact that, if $g$ is analytic in some region and $\,\mathrm{Re}(g)=0$, then $g\equiv ic$, where $c$ real constant.
Hence the function you are looking for is $\mathrm{e}^zz+ic$.