\[ \int \left (d+e x^2\right )^2 \left (a+b \sec ^{-1}(c x)\right ) \, dx \]
Optimal antiderivative \[ d^{2} x \left (a +b \,\mathrm {arcsec}\left (c x \right )\right )+\frac {2 d e \,x^{3} \left (a +b \,\mathrm {arcsec}\left (c x \right )\right )}{3}+\frac {e^{2} x^{5} \left (a +b \,\mathrm {arcsec}\left (c x \right )\right )}{5}-\frac {b \left (120 c^{4} d^{2}+40 c^{2} d e +9 e^{2}\right ) x \arctanh \left (\frac {c x}{\sqrt {c^{2} x^{2}-1}}\right )}{120 c^{4} \sqrt {c^{2} x^{2}}}-\frac {b e \left (40 c^{2} d +9 e \right ) x^{2} \sqrt {c^{2} x^{2}-1}}{120 c^{3} \sqrt {c^{2} x^{2}}}-\frac {b \,e^{2} x^{4} \sqrt {c^{2} x^{2}-1}}{20 c \sqrt {c^{2} x^{2}}} \]
command
integrate((e*x^2+d)^2*(a+b*arcsec(c*x)),x, algorithm="giac")
Giac 1.9.0-11 via sagemath 9.6 output
\[ \text {output too large to display} \]
Giac 1.7.0 via sagemath 9.3 output \[ \text {Timed out} \]_______________________________________________________