The Laplace expansion is in fact the expansion of the inverse distance between two points. Let the points have position vectors \textbf{r} and \textbf{r}' , then the Laplace expansion is \frac{1}{\|\mathbf{r}-\mathbf{r}'\|} = \sum_{\ell=0}^\infty \frac{4\pi}{2\ell+1} \sum_{m=-\ell}^{\ell} (-1)^m \frac{r_^\ell }{r_{\scriptscriptstyle>}^{\ell+1} } Y^{-m}_\ell(\theta, \varphi) Y^m_\ell(\theta', \varphi'). Here \textbf{r} has the spherical polar coordinates (r, \theta, \varphi) and \textbf{r}' has (r', \theta', \varphi') with homogeneous polynomials of degree \ell . Further
r< is min(
r,
r′) and
r> is max(
r,
r′). The function Y^m_\ell is a normalized
spherical harmonic function. The expansion takes a simpler form when written in terms of
solid harmonics, \frac{1}{\|\mathbf{r}-\mathbf{r}'\|} = \sum_{\ell=0}^\infty \sum_{m=-\ell}^\ell (-1)^m I^{-m}_\ell(\mathbf{r}) R^{m}_\ell(\mathbf{r}')\quad\text{with}\quad \|\mathbf{r}\| > \|\mathbf{r}'\|. ==Derivation==