Working on this question, I faced the problem of computing:
$$I_n=\int_1^a \frac{ \tan ^{-1}(t)}{t^{2n}\,\sqrt{t^2-1}}\,dt$$
For a given value of $n \geq 1$, Mathematica does not face any problem and the results are quite simple.
For example, if $a=\sqrt 3$ as in the linked question $$I_{5}=\frac{8346321-4992192 \sqrt{2}+3246848 \sqrt{6}}{58786560}\pi+$$ $$\frac{7867 \sqrt{2}-309123 \cot ^{-1}\left(\sqrt{2}\right)}{1088640}$$
My question is : Is it possible to find the antiderivative even in terms of special functions ?
Although I believe there is no closed form for $I_n$ for any values of $n\ge 0$, in the following post, we will provide a recurrence relation for $I_n$ in terms of well-known, closed-form functions.
We proceed by applying IBP to $I_n$: $$ I_n { = \int_1^a\frac{\tan^{-1}t}{t^{2n}\sqrt{t^2-1}}dt = \int_1^a\frac{\tan^{-1}t}{t^{2n+1}}\frac{tdt}{\sqrt{t^2-1}} = \int_{t=1}^{t=a}\frac{\tan^{-1}t}{t^{2n+1}}d(\sqrt{t^2-1}) \\= \frac{\sqrt{t^2-1}\tan^{-1}t}{t^{2n+1}}\Bigg|_1^a +(2n+1) \int_1^a\frac{\sqrt{t^2-1}\tan^{-1}t}{t^{2n+2}}dt - \int_1^a\frac{\sqrt{t^2-1}}{t^{2n+1}(1+t^2)}dt \\= \frac{\sqrt{a^2-1}\tan^{-1}a}{a^{2n+1}} +(2n+1) \int_1^a\frac{(t^2-1)\tan^{-1}t}{t^{2n+2}\sqrt{t^2-1}}dt \\- \int_0^{\cos^{-1}\frac{1}{a}}\frac{\cos^{2n}\theta-\cos^{2n+2}\theta }{1+\cos^2\theta} d\theta \\= \frac{\sqrt{a^2-1}\tan^{-1}a}{a^{2n+1}} +(2n+1)\left( \int_1^a\frac{\tan^{-1}t}{t^{2n}\sqrt{t^2-1}}dt - \int_1^a\frac{\tan^{-1}t}{t^{2n+2}\sqrt{t^2-1}}dt \right) \\- \int_0^{\cos^{-1}\frac{1}{a}} \frac{2(-1)^n}{1+\cos^2\theta} -\cos^{2n}\theta+\sum_{k=0}^{n-1}2(-1)^{n-k-1}\cos^{2k}\theta d\theta \\= \frac{\sqrt{a^2-1}\tan^{-1}a}{a^{2n+1}} -\sqrt2(-1)^n\tan^{-1}\sqrt{\frac{a^2-1}{2}} \\+(2n+1)\left( I_n - I_{n+1} \right)+J_n+\sum_{k=0}^{n-1}2(-1)^{n-k}J_k, } $$ where $J_n=\int_0^{\cos^{-1}\frac{1}{a}}\cos^{2n}\theta d\theta$. The following recurrence is then deduced: $$ I_{n+1}{= \frac{\sqrt{a^2-1}\tan^{-1}a}{(2n+1)a^{2n+1}} -\frac{\sqrt2(-1)^n\tan^{-1}\sqrt{\frac{a^2-1}{2}}}{2n+1} \\+\frac{2nI_n}{2n+1}+\frac{J_n}{2n+1}+\sum_{k=0}^{n-1}\frac{2(-1)^{n-k}}{2n+1}J_k }. $$ For $J_n$, we can write $$ J_n {= \int_0^{\cos^{-1}\frac{1}{a}}\cos^{2n}\theta d\theta = \int_0^{\cos^{-1}\frac{1}{a}}\cos^{2n-1}\theta \cos \theta d\theta \\=\cos^{2n-1}\theta \sin \theta|_0^{\cos^{-1}\frac{1}{a}}+(2n-1)\int_0^{\cos^{-1}\frac{1}{a}} \sin^2\theta\cos^{2n-2}\theta \\=\frac{\sqrt{a^2-1}}{a^{2n}}+(2n-1)\int_0^{\cos^{-1}\frac{1}{a}} \cos^{2n-2}-\cos^{2n}\theta d\theta \\=\frac{\sqrt{a^2-1}}{a^{2n}}+(2n-1)J_{n-1}-(2n-1)J_n, } $$ which yields $$ J_n=\frac{\sqrt{a^2-1}}{2na^{2n}}+\frac{2n-1}{2n}J_{n-1}\quad,\quad J_0=\cos^{-1}\frac{1}{a}. $$