\[ \int \sin ^3(a+b x) \sin ^3(c+d x) \, dx \]
Optimal antiderivative \[ -\frac {3 \sin \left (a -3 c +\left (b -3 d \right ) x \right )}{32 \left (b -3 d \right )}+\frac {9 \sin \left (a -c +\left (b -d \right ) x \right )}{32 \left (b -d \right )}+\frac {\sin \left (3 a -3 c +3 \left (b -d \right ) x \right )}{96 b -96 d}-\frac {3 \sin \left (3 a -c +\left (3 b -d \right ) x \right )}{32 \left (3 b -d \right )}-\frac {9 \sin \left (a +c +\left (b +d \right ) x \right )}{32 \left (b +d \right )}-\frac {\sin \left (3 a +3 c +3 \left (b +d \right ) x \right )}{96 \left (b +d \right )}+\frac {3 \sin \left (3 a +c +\left (3 b +d \right ) x \right )}{32 \left (3 b +d \right )}+\frac {3 \sin \left (a +3 c +\left (b +3 d \right ) x \right )}{32 \left (b +3 d \right )} \]
command
integrate(sin(b*x+a)**3*sin(d*x+c)**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} \]_____________________________________________________