Solution (source code)

= Solution

Differentiate the <Cauchy-Riemann equations> and use equality of mixed <partial derivatives>:
$$
\boxed{\Delta u=u_{xx}+u_{yy}=(v_y)_x-(v_x)_y=0.}
$$
Thus the real part is a <harmonic function>. The same calculation gives $\Delta v=0$. Every <holomorphic function> has the required smoothness.