Given a sequence of harmonic functions u1, u2, ... on an open connected subset G of the Euclidean space Rn, which are pointwise monotonically nondecreasing in the sense that for every point x of G, then the limit automatically exists in the extended real number line for every x. Harnack's theorem says that the limit either is infinite at every point of G or it is finite at every point of G. In the latter case, the convergence is uniform on compact sets and the limit is a harmonic function on G.[1] The theorem is a corollary of Harnack's inequality.
If un(y) is a Cauchy sequence for any particular value of y, then the Harnack inequality applied to the harmonic function um − un implies, for an arbitrary compact set D containing y, that supD |um − un| is arbitrarily small for sufficiently large m and n. This is exactly the definition of uniform convergence on compact sets.
Having established uniform convergence on compact sets, the harmonicity of the limit is an immediate corollary of the fact that the mean value property (automatically preserved by uniform convergence) fully characterizes harmonic functions among continuous functions.
[2] The proof of uniform convergence on compact sets holds equally well for any linear second-order elliptic partial differential equation, provided that it is linear so that um − un solves the same equation.
The only difference is that the more general Harnack inequality holding for solutions of second-order elliptic PDE must be used, rather than that only for harmonic functions.