Computing convolution using only the definition.

I am stuck computing a convolution because I need to be able to solve it using only the definition of the concept. I need to compute $(f*f)(x)$, where $$f(x) = \frac{1}{x^2+1}$$ So far I've only been able to do this:

$$(f*f)(x) = \int_{-\infty}^\infty f(y)\cdot f(x-y)dy = \int_{-\infty}^\infty \frac{1}{y^2-1}\frac{1}{(x-y)^2+1}dy$$ Expanding the parentheses in the denominator and solving the quadrating equation on $y$ returns the roots $y = x\pm i$, and therefore: $$(f*f)(x) = \int_{-\infty}^\infty \frac{1}{y+1}\frac{1}{y-1}\frac{1}{y-x-i}\frac{1}{y-x+i}dy$$ But I don't know how to keep going to use the residue theorem. How could I transform this into an intgral over a closed path? Any help would be appreciated.

What you should compute is$$\int_{-\infty}^\infty\frac1{y^2+1}\frac1{(x-y)^2+1}\,\mathrm dy.\tag1$$If$$f(z)=\frac1{z^2+1}\frac1{(x-z)^2+1},$$then$$\operatorname{res}_{z=i}f(z)=\frac i{2x(x-2i)}\quad\text{and}\quad\operatorname{res}_{z=x+i}f(z)=-\frac i{2x(x+2i)}.$$So, $(1)$ is equal to$$2\pi i\left(\frac i{2x(x-2i)}-\frac i{2x(x+2i)}\right)=\frac{2\pi}{x^2+4}.$$