Harnack inequality for harmonic functions

If $u\ge 0$ is a harmonic function on a neighbourhood of $\overline{B_{4r}(y)}$, thenFor $x,z\in B_r(y)$, use $B_r(x)\subset B_{3r}(z)\subset B_{4r}(y)$ and the mean value property for harmonic functions. Nonnegativity lets us enlarge the integral, and the ratio of volumes is $3^n$. This compares values without requiring derivative bounds.