I'm trying to work out a problem in John Lee's Introduction to Smooth Manifold:
A harmonic function on a Riemannian manifold is a function $u$ satisfies
$-\Delta u = -\frac{1}{\sqrt{\mathrm{det} g}} \frac{\partial}{\partial x^i}(g^{ij}\sqrt{\mathrm{det}g} \frac{\partial u}{\partial x^j})=0$
(a) If $M$ is a compact, connected Riemannian manifold without boundary, then the only harmonic functions on $M$ are constants.
(b) If $M$ is a compact, connected Riemannian manifold with boundary, and $u,v$ are harmonic functions that agree on the boundary, then $u \equiv v$.
In Why must harmonic functions on compact Riemannian manifolds be constant? Ivo Terek gives a solution to (a). My idea is to use the mean value property of Laplacian equation, but when we pullback $u$ to a chart of $M$, $u$ does not satisfy Laplacian equation unless $g_{ij}=\delta_{ij}$. Does this have something to do with elliptic equation? I've heard of this, but I haven't taken advanced PDE course. How to deal with problem (b) if we don't use the theory of PDE? Thanks in advance.