I'm looking for a closed form for the coefficients of the series expansion of $\arctan(x)$ at some real $x_0$. If $x_0 = 0$, the expansion is of course the well-known
$$\sum_{n\geq 0} \frac{(-1)^n z^{2n+1}}{2n+1}\ .$$
In general, due to analyticity, the expansion is $$\sum_{n\geq 0} \arctan^{(n)}(x_0) \cdot (x - x_0)^n = \\\arctan(x_0) + \sum_{n\geq 0}\left(\frac{d^n}{dx^n}\frac{1}{1+x^2}\right)_{\!|\, x=x_0}\! \cdot (x - x_0)^{n+1}$$
but I'm struggling to find a closed form for these coefficients $$\left(\frac{d^n}{dx^n}\frac{1}{1+x^2}\right)_{\!|\, x=x_0}$$
In fact, I don't necessarily need a closed form as such; a recursive formula would also be all right.
UPDATE: Mathematica says the $n$-th coefficient in the series expansion of $\arctan(x)$ at $x = x_0$ is: $$\frac{-i}{2}\left((-i-x_0)^{-(n+1)} - (i-x_0)^{-(n+1)}\right)$$ One can probably derive this using a partial fraction expansion like Daniel Fischer suggested in his comment and then the geometric series. But I'm struggling to simplify this result to a ‘purely real’ formula.
UPDATE: After some rewriting, I found the following formula: $$\arctan(x) = \arctan(0) + \sum_{n=1}^\infty \frac{\sin(n\arg(i-x_0))}{n(x_0^2+1)^{\frac{n}{2}}} (x - x_0)^n$$
I guess one can simplify the $\arg(i-x_0)$ using $\arctan$. Not sure if the result is the easiest form one can get though.