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