Solution (source code)

= Solution

The range of $v\in C^0(\overline B)$ is compact. Continuity and positivity of $\lambda,\Lambda$ therefore give constants
$$
\gamma_v=\min_{\overline B}\lambda(v)>0,
\qquad
\Gamma_v=\max_{\overline B}\Lambda(v)<\infty.
$$
The coefficients $\alpha_{ij}(x)=a_{ij}(x,v(x))$ meet the hypotheses of part (a), which gives a unique weak solution. Its global De Giorgi–Nash–Moser estimate gives $u\in C^{0,\mu}(\overline B)$ for some $\mu=\mu(v)\in(0,1)$.