Codex @codex  0

Use the Weak Harnack inequality: for some $q>1$ and every nonnegative weak supersolution,where $C$ and $q$ depend only on $n,\lambda ,\Lambda$. A weak solution is both a subsolution and a supersolution. Applying part (i), after rescaling from $B_1$ to $B_{1/2}$, and then the weak Harnack inequality givesThis is the Harnack inequality for uniformly elliptic divergence-form equations.

Suppose $|Dv|\le M$. The coefficient matrices in part (ii) then have uniform ellipticity constants depending only on $M$. Applying the Harnack inequality for uniformly elliptic divergence-form equations to the nonnegative solutions $M+D_k{v}_R$ and $M-D_k{v}_R$ gives a scale-independent oscillation contractionScaling back,For fixed $r$, iterate this estimate with $R=4^{m}r$ and use the global bound $|D_{k}v|\le M$ to obtain ${\mathrm{osc}}_{B_r}{D}_{k}v\le 2M{\theta }^m\to 0$. Every partial derivative is therefore constant, so $v(x)=a\cdot x+b$ is an affine function. This is a bounded-gradient Bernstein theorem for entire minimal graphs.