I'm stuck trying to prove the following identity, given $\zeta \in \mathbb{D}$, where $\mathbb{D} \subset \mathbb{C}$ is the unit disk. Here, $z = x +iy$.
$$ \iint_{\mathbb{D}}\frac{1}{\left|1-\bar{z}\zeta\right|^4}\mathrm dx\mathrm dy = \frac{\pi}{\left(1-|\zeta|^2\right)^2} $$
I was hoping just to expand in $x,y$ and then convert to polar, but things got real nasty and it wasn't working out. I'm now convinced there's some simple trick that I'm missing.
Let $a = \left|{\zeta}\right|$. We have
\begin{equation}I = \int_{\mathbb{D}}^{}\frac{1}{{\left|1-\overline{z} a\right|}^{4}} d x d y = \int_{\mathbb{D}}^{}\frac{1}{{\left({\left(1-a x\right)}^{2}+{\left(a y\right)}^{2}\right)}^{2}} d x d y\end{equation}
Using polar coordinates gives
\begin{equation}I = \int_{0}^{1}r d r \left(\int_{-\pi}^{{\pi}}\frac{1}{{\left(1-2 a r \cos {\theta}+{a}^{2} {r}^{2}\right)}^{2}} d {\theta}\right) \equiv \int_{0}^{1}J \left(a r\right) r d r\end{equation}
Now using the substitution $t = \tan \left({\theta}/2\right)$, one gets
\begin{equation}\renewcommand{\arraystretch}{2.0} \begin{array}{rcl}J \left(k\right)&=&\displaystyle \int_{-\infty}^{\infty }\frac{2 \left(1+{t}^{2}\right)}{{\left(\left(1+{k}^{2}\right) \left(1+{t}^{2}\right)-2 k \left(1-{t}^{2}\right)\right)}^{2}} d t\\ &=&\displaystyle \int_{-\infty}^{\infty }\frac{2 \left(1+{t}^{2}\right)}{{\left({\left(1+k\right)}^{2} {t}^{2}+{\left(1-k\right)}^{2}\right)}^{2}} d t \end{array}\end{equation}
Now using the substitution $t = u \left(1-k\right)/\left(1+k\right)$, this leads to
\begin{equation}\renewcommand{\arraystretch}{2.0} \begin{array}{rcl}J \left(k\right)&=&\displaystyle \frac{2}{{\left(1-k\right)}^{3} {\left(1+k\right)}^{3}} \int_{-\infty}^{\infty }\frac{{\left(1+k\right)}^{2}+{u}^{2} {\left(1-k\right)}^{2}}{{\left(1+{u}^{2}\right)}^{2}} d u\\ &=&\displaystyle \frac{2}{\left(1-k\right) {\left(1+k\right)}^{3}} \int_{-\infty}^{\infty }\frac{1}{1+{u}^{2}} d u+\frac{8 k}{{\left(1-k\right)}^{3} {\left(1+k\right)}^{3}} \int_{-\infty}^{\infty }\frac{1}{{\left(1+{u}^{2}\right)}^{2}} d u\\ &=&\displaystyle \frac{1+{k}^{2}}{{\left(1-{k}^{2}\right)}^{3}} {2 \pi} \end{array}\end{equation}
Hence
\begin{equation}\renewcommand{\arraystretch}{2.0} \begin{array}{rcl}I&=&\displaystyle \frac{{2\pi}}{{a}^{2}} \int_{0}^{a}\frac{1+{k}^{2}}{{\left(1-{k}^{2}\right)}^{3}} k d k\\ &=&\displaystyle \frac{{2\pi}}{{a}^{2}} {\left[\frac{{k}^{2}}{2 {\left(1-{k}^{2}\right)}^{2}}\right]}_{0}^{a}\\ &=&\displaystyle \frac{{\pi}}{{\left(1-{a}^{2}\right)}^{2}} \end{array}\end{equation}