Vamos a exigir la sustitución de $z=x+iy$ e $\bar z=x-iy$. Entonces, tenemos
$$\begin{align}
\frac{\partial \phi}{\partial x}&=\frac{\partial \phi}{\partial z}\frac{\partial z}{\partial x}+\frac{\partial \phi}{\partial \bar z}\frac{\partial \bar z}{\partial x}\\\\
&=\frac{\partial \phi}{\partial z}+\frac{\partial \phi}{\partial \bar z}
\end{align}$$
y
$$\frac{\partial^2 \phi}{\partial x^2}=\frac{\partial^2 \phi}{\partial z^2}+2\frac{\partial^2 \phi}{\partial z \partial \bar z}+\frac{\partial^2 \phi}{\partial \bar z^2} \tag 1$$
Del mismo modo, para la derivada parcial con respecto a $y$, tenemos
$$\frac{\partial^2 \phi}{\partial y^2}=-\frac{\partial^2 \phi}{\partial z^2}+2\frac{\partial^2 \phi}{\partial z \partial \bar z}-\frac{\partial^2 \phi}{\partial \bar z^2} \tag 2$$
La adición de $(1)$ e $(2)$, obtenemos
$$\frac{\partial^2 \phi}{\partial x^2}+\frac{\partial^2 \phi}{\partial y^2}=4\frac{\partial^2 \phi}{\partial z \partial \bar z}=0 \tag 3$$
con lo cual la solución de $(3)$ para $\phi$ revela que $\phi(z,\bar z) = f(z)+g(\bar z)$. Por último, la sustitución de la espalda para $z$ e $\bar z$ rendimientos
$$\phi(x,y)=f(x+iy)+g(x-iy)$$
como iba a ser mostrado.