If $u(x, y)$ is harmonic in the unit disk $\mathbb{D}=\{z=x+i y:|z|<1\}$, define
$$
f(z)=2 u\left(\frac{z}{2}, \frac{z}{2 i}\right)-u(0,0) .
$$
In other words, formally replace the real variables $x$ and $y$ in the two-variable Taylor series for $u$ by the complex variables $\alpha$ and $\beta$ then set $\alpha=z / 2$ and $\beta=z / 2 i$. Where will the resulting series converge? Prove that $f$ is well defined (i.e., that this substitution makes sense) and analytic in $\mathbb{D}$ and $\operatorname{Re} f=u$. Note that this procedure does not work for the functions $u=\log \left(x^2+y^2\right)$ and $\operatorname{Re} \frac{1}{z}$ which are harmonic in $\mathbb{D} \backslash\{0\}$. However, if $u$ is harmonic in a neighborhood of $\left(x_0, y_0\right)$ then, for $z_0=x_0+i y_0$,
$$
f(z)=2 u\left(\left(z-z_0\right) / 2+x_0,\left(z-z_0\right) /(2 i)+y_0\right)-u\left(x_0, y_0\right)
$$
satisfies $\operatorname{Re} f=u$. Try this procedure for a few familiar functions such as $e^x \cos y$ and $\operatorname{Re}\left(z^n\right)$