In *Simon's book Harmonic Analysis*, **example 3.5.12** shows:

Fix $x \in \mathbb{R}^{\nu}$ (with $\left.\nu \geq 3\right), x \neq 0$ and define $f$ on $\{y:| y|<| x |\}$ by
$$
f(y)=|x-y|^{-(\nu-2)}.
$$
Then we have$,$
$$
|x-y|^{-(\nu-2)}=\sum_{d=0}^{\infty} \frac{|y|^{d}}{|x|^{d+\nu-2}} \alpha_{\nu, d} Z_{d}\left(\frac{y}{|y|}, \frac{x}{|x|}\right) \label{1}\tag{$\ast$}
$$
where $\alpha_{\nu, d}$ is a constant which is dependent on $\nu, d$ and the $Z_{d}$ (the zonal harmonic of degree d) term at $y=0$ is interpreted as 1 if $d=0$ and is irrelevant if $d \geq 1$.

**My Question** is$,$ for fixed $x \in \mathbb{R}^{\nu}$ (with $\left.\nu \geq 3\right), x \neq 0$ and function $f(y)=|x-y|^{-\lambda}$ on $\{y|| y|<| x |\}$ with $\lambda>0,$ can we obtain a similar formula like \eqref{1} $?$

Any reference and suggestion will be greatly appreciated. Thanks!