Solution (source code)

= Solution

First repeat the <Caccioppoli inequality> calculation on two concentric <balls>. If $h$ is a <harmonic function> on $B_R$, choose a <cutoff function> equal to one on $B_r$, supported in $B_R$, and with <gradient> bounded by $C/(R-r)$. The argument of (ii), with $A=I$ and $k=0$, gives
$$
\int_{B_r}|\nabla h|^2\,dx\leq\frac{C_*}{(R-r)^2}\int_{B_R}|h|^2\,dx
\qquad(0<r<R\leq1).
$$
Here $C_*$ depends at most on the dimension. When $h$ is only known to be smooth inside $B_R$ and in $L^2(B_R)$, the <weak formulation> and all <integrals> involving the <cutoff function> are still legitimate, since its <support> is strictly inside $B_R$.

Every classical <partial derivative> of the given smooth <harmonic function> $u$ is again a <harmonic function>, since constant-coefficient <partial derivatives> commute with the <Laplacian>. Choose indices $i_1,\ldots,i_m$ representing the <multi-index> $\alpha$, and define
$$
h_0=u,\qquad h_j=\partial_{i_j}\cdots\partial_{i_1}u,\qquad
R_j=1-\frac{j}{2m}\quad(0\leq j\leq m).
$$
For each $j\geq1$, the preceding <Caccioppoli inequality> yields
$$
\int_{B_{R_j}}|h_j|^2
\leq\int_{B_{R_j}}|\nabla h_{j-1}|^2
\leq C_*(2m)^2\int_{B_{R_{j-1}}}|h_{j-1}|^2.
$$
Starting with $u\in L^2(B_1)$, these estimates also establish the integrability needed at each subsequent step. Iterating $m$ times gives
$$
\boxed{\int_{B_{1/2}}|\partial^\alpha u|^2\,dx
\leq\bigl(4C_*m^2\bigr)^m\int_{B_1}|u|^2\,dx.}
$$
This is the desired <interior derivative estimate for a harmonic function>, with a constant depending only on $n$ and $m$, uniformly over every <multi-index> of size $m$.