\[ \int (d+e x)^m \left (a^2+2 a b x+b^2 x^2\right )^3 \, dx \]
Optimal antiderivative \[ \frac {\left (-a e +b d \right )^{6} \left (e x +d \right )^{1+m}}{e^{7} \left (1+m \right )}-\frac {6 b \left (-a e +b d \right )^{5} \left (e x +d \right )^{2+m}}{e^{7} \left (2+m \right )}+\frac {15 b^{2} \left (-a e +b d \right )^{4} \left (e x +d \right )^{3+m}}{e^{7} \left (3+m \right )}-\frac {20 b^{3} \left (-a e +b d \right )^{3} \left (e x +d \right )^{4+m}}{e^{7} \left (4+m \right )}+\frac {15 b^{4} \left (-a e +b d \right )^{2} \left (e x +d \right )^{5+m}}{e^{7} \left (5+m \right )}-\frac {6 b^{5} \left (-a e +b d \right ) \left (e x +d \right )^{6+m}}{e^{7} \left (6+m \right )}+\frac {b^{6} \left (e x +d \right )^{7+m}}{e^{7} \left (7+m \right )} \]
command
integrate((e*x+d)**m*(b**2*x**2+2*a*b*x+a**2)**3,x)
Sympy 1.10.1 under Python 3.10.4 output
\[ \text {output too large to display} \]
Sympy 1.8 under Python 3.8.8 output \[ \text {Timed out} \]_____________________________________________________