Closed-form of $\int_0^1\left(\frac{\arctan x}{x}\right)^n\,dx$

Using substitution and integration by parts, \begin{align}&\int_0^1\left(\frac{\tan^{-1}x}{x}\right)^n\,dx=\int_0^{\pi/4}\frac{x^n}{\tan^n x}\sec^2 x\,dx\\&=\left[\frac{x^n}{\tan^n x}\tan x\right]_0^{\pi/4}-\int_0^{\pi/4}\frac{nx^{n-1}}{\tan^{n-1}x}-\frac{nx^n}{\tan^n x}\sec^2x\,dx\end{align} So we have $$\int_0^1\left(\frac{\tan^{-1}x}{x}\right)^n\,dx=\frac{n}{n-1}\int_0^{\pi/4}\frac{x^{n-1}}{\tan^{n-1}x}\,dx-\frac{1}{n-1}\left(\frac{\pi}{4}\right)^n$$ Now define $\displaystyle I(m,n)=\int_0^{\pi/4}x^m\cot^nx\,dx$ for $m\ge n$. Using integration by parts, we get a recurrence formula. \begin{align}I(m,n)&=\int_0^{\pi/4}x^{m}\cot^{n}x\,dx\\&=\left[\frac{1}{m+1}x^{m+1}\cot^n x\right]_0^{\pi/4}+\frac{n}{m+1}\int_0^{\pi/4}x^{m+1}\cot^{n-1}x(\cot^2x+1)\,dx\end{align} So we have $\displaystyle I(m,n)=-I(m,n-2)+\frac{m}{n-1}I(m-1,n-1)-\frac{1}{n-1}\left(\frac\pi4\right)^m$. We now need to determine $I(m,0)$ and $I(m,1)$ to find $I(n-1,n-1)$. It is easily seen $$I(m,0)=\frac{1}{m+1}\left(\frac\pi4\right)^{m+1}$$ $I(m,1)$ is a bit of a problem. This is where things start to get messy.(All zetas, catalans, and polygammas appear here) \begin{align}I(m,1)&=\int_0^{\pi/4}x^m\cot x\,dx\\&=\left[x^m\log\sin x\right]_0^{\pi/4}-m\int_0^{\pi/4}x^{m-1}\log\sin x\,dx\\&=\frac{\log 2}{2}\left(\frac{\pi}{4}\right)^m+\sum_{n=1}^\infty\frac mn\int_0^{\pi/4}x^{m-1}\cos 2nx\,dx\end{align} Now using integration by parts,$$\int_0^{\pi/4}x^{m-1}\cos 2nx\,dx=\frac{1}{2n}\left(\frac\pi4\right)^{m-1}\sin\frac{n\pi}2+\frac{m-1}{(2n)^2}\left(\frac\pi4\right)^{m-2}\cos\frac{n\pi}2-\frac{(m-1)(m-2)}{(2n)^2}\int_0^{\pi/4}x^{m-3}\cos 2nx\,dx$$ This yields \begin{align}\int_0^{\pi/4}x^{m-1}\cos 2nx\,dx=&\sin\frac{n\pi}{2}\sum_{k=1}^{[\frac{m+1}2]}\frac{(-1)^{k-1}}{(2n)^{2k-1}}\left(\frac\pi4\right)^{m-2k+1}\frac{(m-1)!}{(m-2k+1)!}\\&+\cos\frac{n\pi}{2}\sum_{k=1}^{[\frac m2]}\frac{(-1)^{k-1}}{(2n)^{2k}}\left(\frac\pi4\right)^{m-2k}\frac{(m-1)!}{(m-2k)!}+\left(\frac{(m-1)!}{(2n)^m}(-1)^{m/2}\right)_{2|m}\end{align} where the last term $\displaystyle\left(\frac{(m-1)!}{(2n)^m}(-1)^{m/2}\right)$ is present iff $m$ is even. Now plug this up on the above summation and we get the $I(m,1)$. \begin{align}I(m,1)=&\sum_{k=1}^{[\frac{m+1}2]}\frac{(-1)^{k-1}\pi^{m+1-2k}}{2^{2m+1+4k}}\binom{m}{2k-1}\left(\psi_{2k-1}\left(\frac14\right)-\psi_{2k-1}\left(\frac34\right)\right)\\&+\sum_{k=1}^{[\frac m2]}\frac{(-1)^{k-1}\pi^{m-2k}(2^{2k}-1)}{2^{2m+2k+1}}\frac{m!}{(m-2k)!}\zeta(2k+1)+\left(\frac{(-1)^{m/2}m!}{2^m}\zeta(m+1)\right)_{2|m}\end{align}