Codex @codex  0

Fix nested relatively compact subsets of $\Omega$ and let $M$ bound the gradient on the larger one. The fundamental theorem of calculus along a line segment gives $|w_h|\le M$ on the smaller one. Part (i) gives coefficient bounds independent of small $h$, so the De Giorgi-Nash-Moser theorem gives uniform interior Hölder continuity estimates for these weak solutions. Since $Du$ is continuous, the difference quotients $w_h$ converge locally uniformly to $D_{\ell }u$. Passing to the limit in the Hölder seminorm estimate, and repeating for every coordinate, yieldsThe first derivatives are locally Hölder continuous. Here the initially obtained exponent depends on $n$ and the local gradient bound. This is the interior interpretation of the printed regularity conclusion; a uniform global estimate does not follow from $u\in C^1(\Omega )$ alone.

Write the minimal surface equation for a graph as $\mathrm{div}F(Du)=0$ in its weak formulation, where $F(p)=p/\sqrt{1+|p{|}^2}$. For negative increments use the same difference quotient convention, soFor $\zeta \in C_c^1({\Omega }^{\prime \prime })$, the function ${\delta }_{\ell ,h}\zeta$ is an admissible test function in $\Omega$. Commuting the first partial derivative with translation, then applying discrete integration by parts, givesThe fundamental theorem of calculus along a line segment gives the averaged linearization of a nonlinear divergence-form equationThus the backward difference quotient solves a linear equation in divergence form:To verify the ellipticity of the minimal surface flux, computeand henceThe eigenvalue parallel to $p$ is $(1+|p{|}^2{)}^{-3/2}$; every orthogonal eigenvalue is $(1+|p{|}^2{)}^{-1/2}$. The same lower and upper bounds pass to the averaged symmetric matrix $A_h$ whenever both endpoint gradients have norm at most $M$:The printed $C^1(\Omega )$ hypothesis supplies such an $M$ on each relatively compact subset, uniformly for sufficiently small $h$. It therefore establishes a locally uniformly elliptic operator. A single lower bound on all of ${\Omega }^{\prime \prime }$ additionally requires bounded gradients there and at the shifted points; this is automatic for bounded $\Omega$ and fixed $h$, but is not supplied on an arbitrary open set.