The following fact is quite standard and does not have a very long proof:


$(\\*)$ If $u$ is harmonic on $B\_1(0)\setminus \{0\}$ and uniformly bounded, then $u$ in fact extends to a harmonic function on the whole ball.


Some googling reveals that such statements are in fact true for large classes of elliptic operators with much more general singularity sets and growth conditions. There appears to be a vast literature on the subject.


I would like to know if the exact analogue of $(\\*)$ has a "simple" proof for linear elliptic operators. 
"Simple" is slightly ambiguous but is meant to mean tools present in standard elliptic theory textbooks, i.e. Gilbarg and Trudinger or Jost.



Removable Singularities for Elliptic Equations
The proof I had in mind was actually simpler and appears to work only in dimension greater than 2. Here is a sketch for a second order self-adjoint elliptic operator $Pu = \partial\_i(a^{ij}\partial\_ju)$. Suppose $u$ satisfies $Pu = 0$ on $B\backslash\{0\}$. We want to show that $u$ is a weak solution on $B$. It suffices to show that if $\phi$ is a smooth compactly supported function on $B$, then
$$
\int \phi Pu = 0.
$$
Let $\chi$ be a smooth nonnegative compactly supported function that is equal to $1$ on a neighborhood of $0$ and, given $\epsilon > 0$, let $\chi\_\epsilon(x) = \chi(x/\epsilon)$. Then
$$
\int \phi Pu = \int [\phi(1-\chi\_\epsilon) + \phi\chi\_\epsilon]Pu = \int \phi\chi\_\epsilon Pu
= \int P(\phi\chi\_\epsilon)u.
$$
If we now expand out $P(\phi\chi\_\epsilon)$, we see that it is $O(\epsilon^{-2})$. Therefore, if $u \in L\_\infty$, then
$$
\int P(\phi\chi\_\epsilon)u = O(\epsilon^{n-2})
$$
Therefore, taking the limit $\epsilon \rightarrow 0$ shows that $u$ is a weak solution on $B$. Elliptic regularity tells you that it is in fact a strong smooth solution.


