  As shown recently [Phys. Rev. E 95, 033307 (2017)], spheroidal harmonics
expansions are well suited for the external solution of Laplace's equation for
a point source outside a spherical object. Their intrinsic singularity matches
the line singularity of the analytic continuation of the solution and the
series solution converges much faster than the standard spherical harmonic
solution. Here we extend this approach to internal potentials using the Kelvin
transformation, i.e. radial inversion, of the spheroidal coordinate system.
This transform converts the standard series solution involving regular solid
spherical harmonics into a series of irregular spherical harmonics. We then
substitute the expansion of irregular spherical harmonics in terms of
transformed irregular spheroidal harmonics into the potential. The spheroidal
harmonic solution fits the image line singularity of the solution exactly and
converges much faster. We also discuss why a solution in terms of regular solid
spheroidal harmonics cannot work, even though these functions are finite
everywhere in the sphere. We also present the analogous solution for an
internal point source, and two new relationships between the solid spherical
and spheroidal harmonics.
