Solution.
Let 0<r<R<1 so that B(0,r)⊆B(0,R)⊆D. Since {un}n is a decreasing sequence, we un−um is a non-negative harmonic function for n≤m. Thus, by Harnack's inequality, we have for any z∈B(0,r) that
∣un(z)−um(z)∣=un(z)−um(z)≤R−rR+r(un(0)−um(0))=R−rR+r∣un(0)−um(0)∣.
Since {un(0)}n is a decreasing sequence bounded below (by 0), it is convergent, hence Cauchy. Thus, because r,R are independent of n,m, it follows that {un}n is uniformly Cauchy on compact sets, so it converges locally uniformly to u=infnun. Since each un was continuous, locally uniform convergence implies that u is also continuous, so we may integrate.
Let B(z,R) be such that B(z,R)⊆D. Then because un is harmonic,
un(z)=2π1∫02πun(z+Reiθ)dθ.
By uniform convergence, we may take n→∞ and interchange the limit and integral to get
u(z)=2π1∫02πu(z+Reiθ)dθ,
so u satisfies the mean value property, hence harmonic. We reproduce the proof of this fact below:
Let
ρ(z)=Cexp(−1−∣z∣21)χD(z)
with C such that ∥ρ∥L1=1. Then ρ is a smooth function supported in D. For ε>0, set ρε(z)=ε1ρ(εz) so that ρε is supported on B(0,ε) and set uε=ρε∗u.
Let z∈D and ε small enough so that B(z,ε)⊆D. Then
uε(z)=∫B(0,ε)ρε(w)u(z−w)dw=∫0ε∫02πrρε(reiθ)u(z−reiθ)dθdr=∫0εrρε(r)∫02πu(z−reiθ)dθdr=∫0ε2πrρε(r)u(z)dr=u(z)∫0ε∫02πrρε(reiθ)dθdr=u(z)∥ρε∥L1=u(z).
Thus, uε=u whenever ε is small enough, and because uε is a convolution with a smooth function against a continuous function, it follows that u is also smooth. To show that u is harmonic, we have by the divergence theorem
∫B(z,R)Δu(w)dm=∫∂B(z,R)∇u(z+Reiθ)⋅(cosθ,sinθ)dθ=∫∂B(z,R)∂x∂u(z+Reiθ)cosθ+∂x∂u(z+Reiθ)sinθdθ=∫∂B(z,R)∂r∂u(z+Reiθ)dθ=∂r∂∫∂B(z,R)u(z+Reiθ)dθ=0,
since u satisfies the mean value property. Thus, because u is smooth, Δu is continuous, and so
Δu(z)=R→0limm(B(z,R))1∫B(z,R)Δu(w)dw=0,
so u is harmonic.