Interior derivative estimate for a harmonic function

If $u$ is harmonic on $B_r(x_0)$, then for every integer $m\ge 1$,The first-order estimate follows by differentiating the ball mean-value formula and applying the divergence theorem; iteration gives higher orders.