I found a nice general formula for the Taylor series of $\tan x$:
$$\tan x = \sum_{n\,=\,1}^\infty \frac {(-1)^{n-1}2^{2n} (2^{2n}-1) B_{2n}} {(2n)!} x^{2n - 1} $$
where $B_n$ are the Bernoulli numbers and $|x| < \dfrac {\pi} 2$.
I've tried Googling for a proof but didn't find anything. Hints would be appreciated too.
I am using the typical definition of the Bernoulli numbers:
$$\frac x {e^x - 1} = \sum_{n\,=\,0}^\infty \frac {B_n x^n} {n!}$$
$$\frac z{e^z-1}+\frac z 2=1+\sum_{n=2}^\infty\frac{B_n}{n!}z^n$$ Replace $z$ with $2iz$ to get $$\color{red}{z\cot(z)}=\frac{iz(e^{iz}+e^{-iz})}{e^{iz}-e^{-iz}}=1+\sum_{n=2}^\infty\frac{B_n}{n!}(2iz)^n=1+\sum_{n=1}^\infty\frac{B_{2n}}{(2n)!}(-1)^n(2z)^{2n}$$ Now use following trigonometric formula $$\tan(z)=\cot(z)-2\cot(2z).$$