Derivation
The derivation of this expansion is simple. One writes
We find here the generating function of the Legendre polynomials :
Use of the spherical harmonic addition theorem
gives the desired result.
Read more about this topic: Laplace Expansion (potential)