<s>[INST] Integral $\int_0^{1}\frac{(1+x^2)\ln(1+x^4)}{1+x^4}dx$ 
Find: $$I=\int_0^{1}\frac{(1+x^2)\ln(1+x^4)}{1+x^4}\,dx$$

I need to use the definition of do Trigamma function but I don't see how.
My try was as following : 
Let $y=x^4$ then $dy=4x^{3}\,dx$ so $\;dx=\frac{dy}{4y^{\frac{3}{4}}}$.
$$I=\int_0^{1}\frac{(1+x^{\frac{3}{4}})\ln(1+x)}{x^{\frac{3}{4}}+x^{\frac{7}{4}}}dx$$
Use : $\ln(1+x)=-\sum_{n=1}^{\infty}\frac{(-1)^{n}x^n}{n}$ 
We obtain 
$$J=\int_0^{1}\frac{x^n+x^{(3+4n)/4}}{x^{3/4}+x^{7/4}}dx$$
Then I use definition of digamma function
$$\int_0^{1}\frac{x^{a-1}\ln x}{1+x}\,dx=\beta'(a)$$
Where:  $$2\beta(a)=\psi\left(\frac{1+a}{2}\right)-\psi\left(\frac{a}{2}\right)$$
But I find sum related with digamma function.
So I don't have other ideas to approach it! 
 [/INST] $$\small \boxed{\int_0^1 \frac{(1+x^2)\ln(1+x^4)}{1+x^4}dx=\frac{3\pi}{2\sqrt 2}\ln 2-\frac1{32}\left(\psi_1\left(\frac{1}{8}\right)+\psi_1\left(\frac{3}{8}\right)-\psi_1\left(\frac{5}{8}\right)-\psi_1\left(\frac{7}{8}\right)\right)}$$
To see that write the integral as:
$$I=\int_0^1 \frac{(1+x^2)\ln(1+x^4)}{1+x^4}dx\overset{x\to \frac{1}{x}}=\int_1^\infty \frac{(1+x^2) (\ln(1+x^4)-\ln(x^4))}{1+x^4}dx$$
$$\Rightarrow 2I=\int_0^\infty \frac{(1+x^2)\ln(1+x^4)}{1+x^4}dx-4\int_1^\infty \frac{(1+x^2)\ln x}{1+x^4}dx$$

For the first integral consider:
$$I(a)=\int_0^\infty \frac{(1+x^2)\ln(a+x^4)}{1+x^4}dx\Rightarrow I'(a)=\int_0^\infty \frac{1+x^2}{(a+x^4)(1+x^4)}dx$$
$$=\frac{1}{a-1}\int_0^\infty \frac{1+x^2}{1+x^4}dx-\frac{1}{a-1}\int_0^\infty \frac{1+x^2}{a+x^4}dx$$
$$\int_0^\infty \frac{1+x^2}{1+x^4}dx=\int_0^\infty \frac{1+\frac{1}{x^2}}{\left(x-\frac{1}{x}\right)^2+2}=\frac{1}{\sqrt 2}\arctan\left(\frac{x^2-1}{\sqrt 2 x}\right)\bigg|_0^\infty =\frac{\pi}{\sqrt 2}$$
$$\int_0^\infty \frac{1+x^2}{a+x^4}dx\overset{x=\sqrt[4]at}=\sqrt[4]a\int_0^\infty \frac{1+\sqrt at^2}{a(1+t^4)}dt=a^{-3/4}\int_0^\infty \frac{dt}{1+t^4}+a^{-1/4}\int_0^\infty \frac{t^2dt}{1+t^4}$$
$$\int_0^\infty \frac{1}{1+t^4}dt\overset{t\to \frac{1}{t}}=\int_0^\infty \frac{t^2}{1+t^4}dt=\frac12 \int_0^\infty \frac{1+t^2}{1+t^4}dt=\frac{\pi}{2\sqrt 2}$$
$$\Rightarrow I'(a)= \frac{\pi}{2\sqrt{2}}\left(\frac{2-a^{-1/4}-a^{-3/4}}{a-1}\right) $$
We are looking here for $I(1)$, but $I(0)=0$, (put $x=\frac{1}{x}$ to see that), so:
$$\Rightarrow I(1)= \frac{\pi}{2\sqrt{2}}\int_0^1 \left(\frac{1-a^{-1/4}+1-a^{-3/4}}{a-1}\right)da=\frac{3\pi}{\sqrt 2}\ln 2 $$
Above follows by splitting into two parts and using this.

For the second integral we will use the trigamma function as you mentioned in your question.
$$\int_1^\infty \frac{(1+x^2)\ln x}{1+x^4}dx\overset{x\to \frac{1}{x}}=-\int_0^1 \frac{\ln x}{1+x^4}dx-\int_0^1 \frac{x^2\ln x}{1+x^4}dx$$
$$=\frac1{64} \left(\psi_1\left(\frac{1}{8}\right)+\psi_1\left(\frac{3}{8}\right)-\psi_1\left(\frac{5}{8}\right)-\psi_1\left(\frac{7}{8}\right)\right)$$
One can use trigamma's reflection formula in order to simplify two more terms, but I think it looks nice this way.
</s>