In complex analysis, Liouville's theorem is that every bounded entire function is constant. To prove it, Cauchy intergral formula is used$$f(z) = \frac{1}{2\pi i}\int_C\frac{f(s)}{s-z}ds$$ where C is circle with radius of R, and $|f'(z)|\le \frac{M}{2\pi R}$ by using bounded condition($|f(z)|\le M)$ and Cauchy's differentiation formula. Then, $f'(z)$ goes to $0 $ as R goes to infinity. So, f(z) is constant. The same method can be used in proving Liouville's theorem for harmonic function?
In Partial differential equation, Liouville's theorem is that $u(\mathbf x)$ is constant if $u(\mathbf x)$ is harmonic function in $\Bbb R^N (N=2,3) $ and $u(\mathbf x)$ is bounded by $M>0$.
By mean value principle, $$u(\mathbf x) = \frac{1}{|\partial B_R(\mathbf x)|}\int_{\partial B_R(\mathbf x)}u(s)ds$$ Then, how can I differentiate $u(\mathbf x)$?
Yes the same method can be used. A Cauchy-type estimate on the derivative of a harmonic function holds, and the proof is very much like that for the holomorphic case.
Let $B$ denote the open unit ball in $\mathbb R ^n.$ Claim: There is a constant $C$ such that
$$|\nabla u (0)| \le C\sup_{\partial B} |u|$$
for any function harmonic on a neighborhood of $\overline B.$
Proof: This comes from differentiating through the integral sign in the Poisson integral representation of $u$ on the unit sphere. Doing this will give an appropriate bound on $|\partial u/\partial x_k (0)|,$ $k=1,\dots ,n.$ This proves the claim.
Now suppose $u$ is harmonic on $\mathbb R ^n$ and $|u|\le M$ everywhere. For any $r>0,$ we can set $u_r(x) = u(rx).$ Note each $u_r$ is harmonic on $\mathbb R^n$ and $|u_r|\le M.$ Using the above, we get
$$|\nabla u_r (0)| \le C\sup_{\partial B} |u_r|\le C\cdot M.$$
But $\nabla u_r (0) = r\nabla u (0).$ This shows $|\nabla u (0)| \le CM/r.$ Now let $r\to \infty$ to see $|\nabla u (0)| =0.$ Finally, apply this to any translate of $u$ to see $|\nabla u | =0$ everywhere. This implies $u$ is constant.