I have proven the converse of the MVT before, but I have not seen it done like this and I do not understand it.
Consider the limit:
$ \lim_{r \to 0} \frac{2n}{r^2}[\frac{1}{n \alpha (n)} \int_{\delta B(0,1)} u(x+r\omega)dS(\omega)-u(x)]$
In the above, $\alpha(n)$ represents the volume of the unit ball.
It says that you can use a Taylor expansion with remainder, and note that
$\int_{\delta B(0,1)}x_i=\int_{\delta B(0,1)}x_jx_k $ for $j\neq k$,
then it can be shown that this limit is equal to $ \Delta u(x)$, and then we conclude that the limit is $0$. I cannot begin to see how to reduce this limit to the laplacian of $u$. I am not very familiar with PDEs, but I am trying to learn. Can someone demonstrate how to manipulate this limit?
The converse to the mean value properties say if $$ u(x_0)= \frac1{\mathcal H^{n-1}(\partial B_\rho(x_0))}\int_{\partial B_\rho(x_0)} u(x) \, d\mathcal H^{n-1}_x $$ for all $B_r(x_0) \subset \subset \Omega$ then $\Delta u=0$ in $\Omega$. This directly follows from the Proposition since $$\frac1{\mathcal H^{n-1}(\partial B_\rho(x_0))}\int_{\partial B_\rho(x_0)} u(x) \, d\mathcal H^{n-1}_x -u(x_0)=0 $$for all $B_r(x_0) \subset \subset \Omega$.