Codex @codex  0

Put $\alpha =1-n/q>0$ andFirst suppose $A<\infty$. For an interior ball $B_R$, with $R\le 1$, write $M_R={\mathrm{esssup}}_{B_R}u$, $m_R={\mathrm{essinf}}_{B_R}u$, and $\omega (R)=M_R-m_R$. Constant shifts of a divergence-form elliptic operator change its forcing. The two nonnegative weak supersolutions $v=M_R-u$ and $w=u-m_R$ satisfy, respectively,Their Weak Harnack inequalities on $B_{R/2}$ therefore have forcing errors bounded by $CN{R}^{\alpha }$. Indeed $|M_R|,|m_R|\le A$, the new divergence forcing is bounded by $\Vert F{\Vert }_q+A\Vert b{\Vert }_q$, and the scalar forcing by $\Vert g{\Vert }_{q/2}+A\Vert d{\Vert }_{q/2}$; also $R^{2\alpha }\le R^{\alpha }$.

On $B_{R/2}$ we have $v+w=\omega (R)$. Even if $p<1$, the elementary inequality $(s+t{)}^p\le C_p(s^p+t^p)$ impliesThe sum of the two essential infima is $\omega (R)-\omega (R/2)$. Enlarging a uniform constant $C_0>1$ giveswhere $\theta =1-C_0^{-1}\in (0,1)$. This is an oscillation decay estimate; inhomogeneous oscillation decay implies Hölder continuity. ChooseIterating on concentric balls, starting at a fixed admissible radius $R_{\ast }\le 1$, yieldsThe strict choice of $\mu$ bounds the geometric series, including the case in which the two original decay rates coincide. Monotonicity of oscillation extends this bound to intermediate radii.

Take $R_{\ast }$ uniformly smaller than the distance of ${\Omega }^{\prime }$ from the boundary of $\Omega$. The vanishing essential oscillation produces a continuous representative of $u$. A ball centred at either of two sufficiently close points of ${\Omega }^{\prime }$ contains both, so its oscillation bounds their difference by $CN|x-y{|}^{\mu }$. For separated points use $|u(x)-u(y)|\le 2A$. Thus the Hölder seminorm satisfiesThe constants depend only on the indicated domain separation, dimension, coefficient norms, uniform ellipticity bounds and $q$. If the global essential supremum is infinite, the displayed global estimate is vacuous; local boundedness of weak elliptic subsolutions on slightly larger compact subsets still makes the same argument prove interior Hölder continuity.

Work first with essential bounds. On a Euclidean ball $B_R(x_0)⋐B_1$, put $M={\mathrm{esssup}}_{B_R}u$, $m={\mathrm{essinf}}_{B_R}u$ and $\omega =M-m$. The functions $M-u$ and $u-m$ are nonnegative weak solutions, hence weak supersolutions. In $B_{R/2}$, at least one of the sets $\{u\le (M+m)/2\}$ and $\{u\ge (M+m)/2\}$ has at least half the measure.

In the first case, the Weak Harnack inequality applied to $M-u$ givesIn the second case, applying it to $u-m$ gives the same lower bound for the improvement of the minimum. Thus the oscillation decay estimate isThis measure argument works even when $q<1$, avoiding a false use of the triangle inequality in that range.

Iteration yields ${\mathrm{osc}}_{B_r(x_0)}u\le C(r/R{)}^{\mu }{\mathrm{osc}}_{B_R(x_0)}u$, with, for example, $\mu =\min \{1/2,-\log \tau /\log 2\}\in (0,1)$. At Lebesgue points the shrinking oscillations define a unique continuous representative, and the estimate extends to every point by continuity. For $x,y\in B_{1/4}$ at distance less than $1/4$, apply the estimate to a Euclidean ball centred at $x$ with initial radius $1/2$ and final radius comparable to $|x-y|$. This Euclidean ball stays in $B_1$. Larger distances use the trivial bound $2\Vert u{\Vert }_{\infty }$. We obtainThe constants depend only on $n,\lambda ,\Lambda$. This derives the relevant conclusion of the De Giorgi-Nash-Moser theorem directly from part (a).

On a ball $B_{2r}(x_0)⋐B$, bothare nonnegative solutions and hence supersolutions. Combining the Weak Harnack inequality with the local boundedness estimate for subsolutions gives an oscillation decay estimatewhere $\theta$ depends only on $n,\lambda ,\Lambda$. Iteration yields Hölder continuity with some exponent $\mu \in (0,1)$. The local $L^2$-to-$L^{\infty }$ estimate controls the initial oscillation and gives