If f(z)=u(x,y)+iv(x,y) is analytic in some region of the z plane, then at every point in the region the Cauchy-Riemann conditions are satisfied:
\frac{ \partial u }{ \partial x } = \frac{ \partial v }{ \partial y } \quad \text{and} \quad \frac{ \partial u }{ \partial y } = -\frac{ \partial v }{ \partial x }
and therefor:
\frac{ \partial^2 u }{ \partial x^2 } = \frac{ \partial^2 v }{ \partial x \partial y } \quad \text{and} \quad \frac{ \partial^2 u }{ \partial y^2 } = -\frac{ \partial^2 v }{ \partial y \partial x }
provided these second derivatives exist. In fact,one can show that if f(z) is analytic in some region R, all its derivative exist and are continuous in R. Equating the two cross terms, we obtain
\frac{ \partial^2 u }{ \partial x^2 } + \frac{ \partial^2 u }{ \partial y^2 } = 0
throughout the region R.
Similarly, by differentiating the first of the Cauchy-Riemann equations with respect to y, the second with respect to x, and subtracting we obtain
\frac{ \partial^2 v }{ \partial x^2 } + \frac{ \partial^2 v }{ \partial y^2 } = 0
Equations (6.12a) and (6.12b) are Laplace's partial differential equations in two independent variables x and y. Any function that has continuous partial derivatives of second order and that satisfies Laplace's equation is called a harmonic function.
We have shown that if f(z)=u(x,y)+iv(x,y) is analytic, then both u and v are harmonic functions. They are called conjugate harmonic functions.