\[ \int x^3 \left (a+b x^2\right )^3 \left (A+B x+C x^2+D x^3\right ) \, dx \]
Optimal antiderivative \[ \frac {a^{3} A \,x^{4}}{4}+\frac {a^{3} B \,x^{5}}{5}+\frac {a^{2} \left (3 A b +a C \right ) x^{6}}{6}+\frac {a^{2} \left (3 b B +a D\right ) x^{7}}{7}+\frac {3 a b \left (A b +a C \right ) x^{8}}{8}+\frac {a b \left (b B +a D\right ) x^{9}}{3}+\frac {b^{2} \left (A b +3 a C \right ) x^{10}}{10}+\frac {b^{2} \left (b B +3 a D\right ) x^{11}}{11}+\frac {b^{3} C \,x^{12}}{12}+\frac {b^{3} D x^{13}}{13} \]
command
integrate(x^3*(b*x^2+a)^3*(D*x^3+C*x^2+B*x+A),x, algorithm="fricas")
Fricas 1.3.8 (sbcl 2.2.11.debian) via sagemath 9.6 output
\[ \frac {1}{13} \, D b^{3} x^{13} + \frac {1}{12} \, C b^{3} x^{12} + \frac {1}{11} \, {\left (3 \, D a b^{2} + B b^{3}\right )} x^{11} + \frac {1}{10} \, {\left (3 \, C a b^{2} + A b^{3}\right )} x^{10} + \frac {1}{3} \, {\left (D a^{2} b + B a b^{2}\right )} x^{9} + \frac {1}{5} \, B a^{3} x^{5} + \frac {3}{8} \, {\left (C a^{2} b + A a b^{2}\right )} x^{8} + \frac {1}{4} \, A a^{3} x^{4} + \frac {1}{7} \, {\left (D a^{3} + 3 \, B a^{2} b\right )} x^{7} + \frac {1}{6} \, {\left (C a^{3} + 3 \, A a^{2} b\right )} x^{6} \]
Fricas 1.3.7 via sagemath 9.3 output
\[ \text {Exception raised: TypeError} \]