 Solution

The three-dimensional fundamental solution of the Laplace equation, with $\Delta F=+\delta$, isIt is harmonic away from the source, and its gradient has total outward flux one through a small source-centred sphere.

For $s=|{\mathbf{r}}^{\prime }|\in (0,1)$, place the image at ${\mathbf{r}}^{\ast }={\mathbf{r}}^{\prime }/s^2$, outside the unit ball. On $|\mathbf{r}|=1$, expansion of squared distances gives $|\mathbf{r}-{\mathbf{r}}^{\ast }|=|\mathbf{r}-{\mathbf{r}}^{\prime }|/s$. The method of images therefore gives the Dirichlet Green function inside a sphereThe boundary value is zero by the distance identity, and the image term is harmonic throughout the physical ball. Thus $\Delta G_0=\delta (\mathbf{r}-{\mathbf{r}}^{\prime })$ there. At the central source, take the finite limit of the weighted image term: