There seems to be a removable singularity theorem for harmonic functions that somewhat goes like:
Let $\Omega\subset \Bbb R^n$ be an open region and $x_0\in\Omega$, $u$ is harmonic in $\Omega-\{x_0\}$. Now if we have $u(x)=o(\Phi(x-x_0))$ around $x_0$ then there exists $\tilde u$ harmonic in all of $\Omega$ and $u=\tilde u$ in $\Omega-\{x_0\}$.
Here $\Phi$ is the fundamental solution to $\Delta u=0$ in $\Bbb R^n$, and, since dropping the normalising constant doesn't affect our statement above, you may just take $\Phi=-\log|x-x_0|$ in 2D and $=1/|x-x_0|^{n-2}$ in 3D or higher.
I have two questions concerning this theorem:
1). I can't find a satisfactory proof. So far I only found this one, which however doesn't look okay to me because there seems to be an abuse of $\epsilon-\delta$ language: when the author concluded the proof by letting $\epsilon\to 0$, he seemed to forget the previous part of his proof depended on $\delta$ which depended on $\epsilon$. So is there any other available proof?
2). I think the statement of the theorem is unnecessarily complicated. If $x_0$ is really removable for $u$, then I think it's necessary $u$ is bounded around $x_0$. And then boundedness clearly implies the condition that $u(x)=o(\Phi(x-x_0))$. So why not just say "if $u$ is bounded around $x_0$ then it's a removable singularity"? Why all the fuss about comparing $u$ against $\Phi$?