I wish to evaluate
$$I=\int_{-\infty}^{\infty}\frac{dx}{(x^4+a^4)^2}$$
for a positive real number $a$. I believe the answer is $\frac{3\sqrt{2}\pi}{8a^7}$.
But I am not sure how to arrive at this. I know that $I=2\pi i\sum_{Im(s)>0}\text{res}_f(s)$ where $f=\frac{1}{(x^4+a^4)^2}$, but I am having a hard time directly computing these residues. I know that I only need to consider the residues at the points $ai(1+i)/\sqrt{2}$ and $a(1+i)/\sqrt{2}$ since these are the two poles where $Im(s)>0$, but I cannot easily compute the residues at these two points. I know that these are poles of order $2$, so I can take the holomorphic continuation $g$ of $(z-s)^2f(z)$ and compute $\text{res}_f(s)=g'(a)$, but this becomes so computationally intensive that I am certain there must a better method.
Ideally, I would like to see a residue-theoretic approach to the solution.