\[ \int \cos ^5(c+d x) \sin ^n(c+d x) (a+b \sin (c+d x))^2 \, dx \]
Optimal antiderivative \[ \frac {a^{2} \left (\sin ^{1+n}\left (d x +c \right )\right )}{d \left (1+n \right )}+\frac {2 a b \left (\sin ^{2+n}\left (d x +c \right )\right )}{d \left (2+n \right )}-\frac {\left (2 a^{2}-b^{2}\right ) \left (\sin ^{3+n}\left (d x +c \right )\right )}{d \left (3+n \right )}-\frac {4 a b \left (\sin ^{4+n}\left (d x +c \right )\right )}{d \left (4+n \right )}+\frac {\left (a^{2}-2 b^{2}\right ) \left (\sin ^{5+n}\left (d x +c \right )\right )}{d \left (5+n \right )}+\frac {2 a b \left (\sin ^{6+n}\left (d x +c \right )\right )}{d \left (6+n \right )}+\frac {b^{2} \left (\sin ^{7+n}\left (d x +c \right )\right )}{d \left (7+n \right )} \]
command
integrate(cos(d*x+c)**5*sin(d*x+c)**n*(a+b*sin(d*x+c))**2,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} \]_____________________________________________________