Integral $\int_a^\infty \frac{\arctan(x+b)}{x^2+c}dx$ I was playing around with some integrals and noticed that some integrals of the form:  $$I(a,b,c)=\int_a^\infty \frac{\arctan(x+b)}{x^2+c}dx$$ Does have a closed form.
I am trying to find for what constant $c$ will this work. In case you wonder why only $c$ is problematic, I will try to show by an example. 
$$I=I(1,3,16)=\int_1^\infty \frac{\arctan(x+3)}{x^2+16}dx$$
Let's start by letting $x-1=t\,$ thus:
$$I=\int_0^\infty \frac{{\arctan(\color{blue}{t+4})}}{t^2+2t+17}dt$$
With $\displaystyle{t=\frac{17}{y}\rightarrow dt=-\frac{17}{y^2}dy}$
$$I=\int_0^\infty \frac{{\arctan\left(\color{red}{\frac{17}{y}+4}\right)}}{\left(\frac{17}{y}\right)^2 +\frac{34}{y}+17 }\frac{17}{y^2}dy\overset{y=t}=\int_0^\infty \frac{{\arctan\left(\color{red}{\frac{17}{t}+4}\right)}}{t^2+2t+17}dt$$
$$2I=\int_0^\infty \frac{{\arctan(\color{blue}{t+4})+{\arctan\left(\color{red}{\frac{17}{t}+4}\right)}}}{t^2+2t+17}dt$$
$${\arctan(\color{blue}{t+4})+{\arctan\left(\color{red}{\frac{17}{t}+4}\right)}}=\arctan\left(\frac{\color{blue}{t+4}+\color{red}{\frac{17}{t}+4}}{1-(\color{blue}{t+4})\left(\color{red}{\frac{17}{t}+4}\right)}\right)$$
$$=\arctan\left(\frac{x^2+8x+17}{x}\frac{x}{-4(x^2+8x+17}\right)=\pi-\arctan\left(\frac14\right)$$
Above follows since the original integral is positive so we take $\arctan(-x)$ as $\pi-\arctan x $ and therefore getting a negative answer will not be an issue.
$$I=\frac12 \left(\pi -\arctan\left(\frac14\right)\right)\int_0^\infty \frac{1}{t^2+2t+17}dt$$
Well, now the inner integral is not hard to compute and the final answer happens to be: $$I=\frac12 \left(\pi -\arctan\left(\frac14\right)\right)\frac14\arctan\left(\frac{t+1}{4}\right)\bigg|_0^\infty =\frac{\pi^2}{16}-\frac{3\pi}{16}\arctan\left(\frac14\right)+\frac18\arctan^2\left(\frac14\right) $$
There are more examples that I found by checking and try such as: 
$$I(1,2,9)=\int_1^\infty \frac{\arctan(x+2)}{x^2+9}dx$$
$$I(2,1,6)=\int_2^\infty \frac{\arctan(x+1)}{x^2+6}dx$$
$$I(2,2,13)=\int_2^\infty \frac{\arctan(x+2)}{x^2+13}dx$$
And so on... All those can be solved by the same method: First substitute $x-a=t$, then let $t=\frac{\alpha}{y}$, where $\alpha$ is the "free of x" coefficient from the denominator. 
The problem is that I tried more than $100$ combinations to get those integrals which is not that nice. How can we "smartly" find $c$ so that $I(a,b,c)$ is evaluable by symmetry? Or put in other word what should be $c$ if one wants to compute by symmetry $I(7,13,c)$? 
 A: The simplest approaches are sometimes the right ones. Therefore choose any arbitrarily numbers $a,b$ and $c$ and apply your algorithm. So by firstly  setting $t=x-a$ we get
$$\begin{align}
I(a,b,c)=\int_a^{\infty}\frac{\arctan(x+b)}{x^2+c}dx=\int_0^{\infty}\frac{\arctan(t+a+b)}{t^2+2at+(c+a^2)}dt
\end{align}$$
Now set $\alpha=c+a^2$ and then $\displaystyle t=\frac{c+a^2}{y}$ to further get
$$\begin{align}
I(a,b,c)=\int_0^{\infty}\frac{\arctan(t+a+b)}{t^2+2at+(c+a^2)}dt&=\int_0^{\infty}\frac{\arctan\left(\frac{c+a^2}{y}+a+b\right)}{\left(\frac{c+a^2}{y}\right)^2+2a\left(\frac{c+a^2}{y}\right)+(c+a^2)}\frac{c+a^2}{y^2}dy\\
&\stackrel{y=t}{=}\int_0^{\infty}\frac{\arctan\left(\frac{c+a^2}{y}+a+b\right)}{t^2+2at+(c+a^2)}dt
\end{align}$$
Adding the first and the second form up results in using the addition theorem of the inverse tangent function. This addition looks like the following
$$\small\begin{align}
\arctan(\color{blue}{t+a+b})+\arctan\left(\color{red}{\frac{c+a^2}{y}+a+b}\right)&=\arctan\left(\frac{\color{blue}{t+a+b}+\color{red}{\frac{c+a^2}{y}+a+b}}{1-(\color{blue}{t+a+b})\left(\color{red}{\frac{c+a^2}{y}+a+b}\right)}\right)\\
&=\arctan\left(\frac{t^2+2(a+b)t+a^2+c}{-(a+b)\left(t^2+\frac{2a^2+2ab+b^2+c-1}{a+b}t+a^2+c\right)}\right)
\end{align}$$
In order to make the polynomial $t^2+2(a+b)t+a^2+c$ vanish the following condition has to be fulfilled
$$\frac{2a^2+2ab+b^2+c-1}{a+b}=2(a+b)$$
From hereon we can deduce a relation between $a,b$ and $c$ that has to be satisfied. To be precise

$$\begin{align}
\frac{2a^2+2ab+b^2+c-1}{a+b}=2(a+b)&\Leftrightarrow 2a^2+2ab+b^2+c-1=2(a+b)^2\\
&\Leftrightarrow 2a^2+2ab+b^2+c-1=2a^2+2b^2+4ab\\
&\Leftrightarrow c=b^2+2ab+1
\end{align}$$

And indeed your given case $I(\color{red}{1},\color{blue}{3},\color{green}{16})$ fulfills this relation as $\color{green}{16}=\color{blue}{3}^2+2\cdot\color{blue}{3}\cdot\color{red}{1}+1$. Another example would be $I(\color{red}{1},\color{blue}{1},\color{green}{4})$ which again turns out to work since $\color{green}{4}=\color{blue}{1}^2+2\cdot\color{blue}{1}\cdot\color{red}{1}+1$. 

Therefore for your given case $I(7,13,c)$ you have to choose $c=352$ in order to let the integral be solvable via symmetry.

To finish the evaluation of $I(a,b,c)$, now under the restriction of $c=b^2+2ab+1$, we arrive at
$$\small\arctan\left(\frac{t^2+2(a+b)t+a^2+c}{-(a+b)\left(t^2+\frac{2a^2+2ab+b^2+c-1}{a+b}t+a^2+c\right)}\right)=\arctan\left(\frac{-1}{a+b}\right)=\pi-\arctan\left(\frac1{a+b}\right)$$
following the same argumentation as you did. So for $I(a,b,c)$ as a whole we get
$$\begin{align}
2I(a,b,c)&=\left(\pi-\arctan\left(\frac1{a+b}\right)\right)\int_0^{\infty}\frac{dt}{t^2+2at+(c+a^2)}\\
&=\left(\pi-\arctan\left(\frac1{a+b}\right)\right)\int_0^{\infty}\frac{dt}{(t+a)^2+c}\\
&=\left(\pi-\arctan\left(\frac1{a+b}\right)\right)\left[\frac1{\sqrt{c}}\arctan\left(\frac{t+a}{\sqrt{c}}\right)\right]_0^{\infty}\\
\Leftrightarrow I(a,b,c)&=\frac1{2\sqrt{c}}\left(\pi-\arctan\left(\frac1{a+b}\right)\right)\left[\frac{\pi}2-\arctan\left(\frac a{\sqrt{c}}\right)\right]
\end{align}$$
Where the final formula produces the right value for your example integral $I(1,3,16)$.

I do not claim that this deduced relation between $a,b$ and $c$ is the only one for which the integral can be evaluated via symmetry but in fact it is one possibilty.
