\[ \int x^5 \left (d+e x^2\right ) \left (a+b \sec ^{-1}(c x)\right ) \, dx \]
Optimal antiderivative \[ \frac {d \,x^{6} \left (a +b \,\mathrm {arcsec}\left (c x \right )\right )}{6}+\frac {e \,x^{8} \left (a +b \,\mathrm {arcsec}\left (c x \right )\right )}{8}-\frac {b \left (8 c^{2} d +9 e \right ) x \left (c^{2} x^{2}-1\right )^{\frac {3}{2}}}{72 c^{7} \sqrt {c^{2} x^{2}}}-\frac {b \left (4 c^{2} d +9 e \right ) x \left (c^{2} x^{2}-1\right )^{\frac {5}{2}}}{120 c^{7} \sqrt {c^{2} x^{2}}}-\frac {b e x \left (c^{2} x^{2}-1\right )^{\frac {7}{2}}}{56 c^{7} \sqrt {c^{2} x^{2}}}-\frac {b \left (4 c^{2} d +3 e \right ) x \sqrt {c^{2} x^{2}-1}}{24 c^{7} \sqrt {c^{2} x^{2}}} \]
command
integrate(x^5*(e*x^2+d)*(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} \]_______________________________________________________