Ясно, что относительно оси Оz имеет место аксиальная симметрия. Можно выписать общее решение уравнения Лапласа

в таком случае. Учитывая, что в начале координат потенциал должен быть конечен, а на бесконечности равен нулю, для двух областей пространства выпишем решение:
=\sum\limits_{n=0}^{n=\infty}A_{n}r^{n}P_{n}(\cos\theta))
при

=\sum\limits_{n=0}^{n=\infty}\frac{B_n}{r^{n+1}}P_{n}(\cos\theta))
при

Осталось только найти коэффициенты

. Для этого сначала воспользуемся тем, что потенциал непрерывен:
=\varphi_{2}(R,\theta),)
откуда получим

. Далее, в силу симметрии задачи потенциал в плоскости z=0 внутри сферы следует положить равным 0. Но тогда согласно свойствам полиномов Лежандра в разложении останутся только полиномы с нечетными номерами

. Воспользуемся еще тем, что при переходе через сферу нормальная компонента поля терпит разрыв:

при
)
. Получится вроде бы так:
A_{2l+1}R^{2l}+(2l+2)B_{2l+1}R^{-2l-3})P_{2l+1}(\cos\theta))
.
Теперь домножим это соотношение на
\sin\theta)
и проинтегрируем от 0 до

, используя тождества
P_{m}(x)\,dx=\frac{1}{2n+1}\delta_{mn})
при одновременно нечетных m и n
\,dx=\frac{(-1)^{k}(2k)!}{2^{2k+1}k!(k+1)!})
Написанного здесь достаточно, чтобы получить коэффициенты и записать ответ в виде ряда.