\(\int (c+d x)^3 \cos ^3(a+b x) \sin ^3(a+b x) \, dx\) [156]

Optimal result
Mathematica [A] (verified)
Rubi [A] (verified)
Maple [A] (verified)
Fricas [B] (verification not implemented)
Sympy [B] (verification not implemented)
Maxima [B] (verification not implemented)
Giac [A] (verification not implemented)
Mupad [B] (verification not implemented)
Reduce [B] (verification not implemented)

Optimal result

Integrand size = 24, antiderivative size = 181 \[ \int (c+d x)^3 \cos ^3(a+b x) \sin ^3(a+b x) \, dx=\frac {9 d^2 (c+d x) \cos (2 a+2 b x)}{128 b^3}-\frac {3 (c+d x)^3 \cos (2 a+2 b x)}{64 b}-\frac {d^2 (c+d x) \cos (6 a+6 b x)}{1152 b^3}+\frac {(c+d x)^3 \cos (6 a+6 b x)}{192 b}-\frac {9 d^3 \sin (2 a+2 b x)}{256 b^4}+\frac {9 d (c+d x)^2 \sin (2 a+2 b x)}{128 b^2}+\frac {d^3 \sin (6 a+6 b x)}{6912 b^4}-\frac {d (c+d x)^2 \sin (6 a+6 b x)}{384 b^2} \] Output:

9/128*d^2*(d*x+c)*cos(2*b*x+2*a)/b^3-3/64*(d*x+c)^3*cos(2*b*x+2*a)/b-1/115 
2*d^2*(d*x+c)*cos(6*b*x+6*a)/b^3+1/192*(d*x+c)^3*cos(6*b*x+6*a)/b-9/256*d^ 
3*sin(2*b*x+2*a)/b^4+9/128*d*(d*x+c)^2*sin(2*b*x+2*a)/b^2+1/6912*d^3*sin(6 
*b*x+6*a)/b^4-1/384*d*(d*x+c)^2*sin(6*b*x+6*a)/b^2
 

Mathematica [A] (verified)

Time = 1.58 (sec) , antiderivative size = 132, normalized size of antiderivative = 0.73 \[ \int (c+d x)^3 \cos ^3(a+b x) \sin ^3(a+b x) \, dx=\frac {-324 b (c+d x) \left (-3 d^2+2 b^2 (c+d x)^2\right ) \cos (2 (a+b x))+12 b (c+d x) \left (-d^2+6 b^2 (c+d x)^2\right ) \cos (6 (a+b x))-4 d \left (121 d^2-234 b^2 (c+d x)^2+\left (-d^2+18 b^2 (c+d x)^2\right ) \cos (4 (a+b x))\right ) \sin (2 (a+b x))}{13824 b^4} \] Input:

Integrate[(c + d*x)^3*Cos[a + b*x]^3*Sin[a + b*x]^3,x]
 

Output:

(-324*b*(c + d*x)*(-3*d^2 + 2*b^2*(c + d*x)^2)*Cos[2*(a + b*x)] + 12*b*(c 
+ d*x)*(-d^2 + 6*b^2*(c + d*x)^2)*Cos[6*(a + b*x)] - 4*d*(121*d^2 - 234*b^ 
2*(c + d*x)^2 + (-d^2 + 18*b^2*(c + d*x)^2)*Cos[4*(a + b*x)])*Sin[2*(a + b 
*x)])/(13824*b^4)
 

Rubi [A] (verified)

Time = 0.43 (sec) , antiderivative size = 181, normalized size of antiderivative = 1.00, number of steps used = 2, number of rules used = 2, \(\frac {\text {number of rules}}{\text {integrand size}}\) = 0.083, Rules used = {4906, 2009}

Below are the steps used by Rubi to obtain the solution. The rule number used for the transformation is given above next to the arrow. The rules definitions used are listed below.

\(\displaystyle \int (c+d x)^3 \sin ^3(a+b x) \cos ^3(a+b x) \, dx\)

\(\Big \downarrow \) 4906

\(\displaystyle \int \left (\frac {3}{32} (c+d x)^3 \sin (2 a+2 b x)-\frac {1}{32} (c+d x)^3 \sin (6 a+6 b x)\right )dx\)

\(\Big \downarrow \) 2009

\(\displaystyle -\frac {9 d^3 \sin (2 a+2 b x)}{256 b^4}+\frac {d^3 \sin (6 a+6 b x)}{6912 b^4}+\frac {9 d^2 (c+d x) \cos (2 a+2 b x)}{128 b^3}-\frac {d^2 (c+d x) \cos (6 a+6 b x)}{1152 b^3}+\frac {9 d (c+d x)^2 \sin (2 a+2 b x)}{128 b^2}-\frac {d (c+d x)^2 \sin (6 a+6 b x)}{384 b^2}-\frac {3 (c+d x)^3 \cos (2 a+2 b x)}{64 b}+\frac {(c+d x)^3 \cos (6 a+6 b x)}{192 b}\)

Input:

Int[(c + d*x)^3*Cos[a + b*x]^3*Sin[a + b*x]^3,x]
 

Output:

(9*d^2*(c + d*x)*Cos[2*a + 2*b*x])/(128*b^3) - (3*(c + d*x)^3*Cos[2*a + 2* 
b*x])/(64*b) - (d^2*(c + d*x)*Cos[6*a + 6*b*x])/(1152*b^3) + ((c + d*x)^3* 
Cos[6*a + 6*b*x])/(192*b) - (9*d^3*Sin[2*a + 2*b*x])/(256*b^4) + (9*d*(c + 
 d*x)^2*Sin[2*a + 2*b*x])/(128*b^2) + (d^3*Sin[6*a + 6*b*x])/(6912*b^4) - 
(d*(c + d*x)^2*Sin[6*a + 6*b*x])/(384*b^2)
 

Defintions of rubi rules used

rule 2009
Int[u_, x_Symbol] :> Simp[IntSum[u, x], x] /; SumQ[u]
 

rule 4906
Int[Cos[(a_.) + (b_.)*(x_)]^(p_.)*((c_.) + (d_.)*(x_))^(m_.)*Sin[(a_.) + (b 
_.)*(x_)]^(n_.), x_Symbol] :> Int[ExpandTrigReduce[(c + d*x)^m, Sin[a + b*x 
]^n*Cos[a + b*x]^p, x], x] /; FreeQ[{a, b, c, d, m}, x] && IGtQ[n, 0] && IG 
tQ[p, 0]
 
Maple [A] (verified)

Time = 2.82 (sec) , antiderivative size = 148, normalized size of antiderivative = 0.82

method result size
parallelrisch \(\frac {-54 \left (\left (d x +c \right )^{2} b^{2}-\frac {3 d^{2}}{2}\right ) \left (d x +c \right ) b \cos \left (2 b x +2 a \right )+6 \left (\left (d x +c \right )^{2} b^{2}-\frac {d^{2}}{6}\right ) \left (d x +c \right ) b \cos \left (6 b x +6 a \right )+81 d \left (\left (d x +c \right )^{2} b^{2}-\frac {d^{2}}{2}\right ) \sin \left (2 b x +2 a \right )-3 d \left (\left (d x +c \right )^{2} b^{2}-\frac {d^{2}}{18}\right ) \sin \left (6 b x +6 a \right )+48 b^{3} c^{3}-80 c \,d^{2} b}{1152 b^{4}}\) \(148\)
risch \(\frac {\left (6 b^{2} d^{3} x^{3}+18 b^{2} c \,d^{2} x^{2}+18 b^{2} c^{2} d x +6 b^{2} c^{3}-d^{3} x -c \,d^{2}\right ) \cos \left (6 b x +6 a \right )}{1152 b^{3}}-\frac {d \left (18 x^{2} d^{2} b^{2}+36 b^{2} c d x +18 b^{2} c^{2}-d^{2}\right ) \sin \left (6 b x +6 a \right )}{6912 b^{4}}-\frac {3 \left (2 b^{2} d^{3} x^{3}+6 b^{2} c \,d^{2} x^{2}+6 b^{2} c^{2} d x +2 b^{2} c^{3}-3 d^{3} x -3 c \,d^{2}\right ) \cos \left (2 b x +2 a \right )}{128 b^{3}}+\frac {9 d \left (2 x^{2} d^{2} b^{2}+4 b^{2} c d x +2 b^{2} c^{2}-d^{2}\right ) \sin \left (2 b x +2 a \right )}{256 b^{4}}\) \(234\)
orering \(\frac {5 d \left (18 d^{4} x^{4} b^{4}+72 b^{4} c \,d^{3} x^{3}+108 b^{4} c^{2} d^{2} x^{2}+72 b^{4} c^{3} d x +18 b^{4} c^{4}-11 b^{2} d^{4} x^{2}-22 b^{2} c \,d^{3} x -11 b^{2} c^{2} d^{2}-9 d^{4}\right ) \cos \left (b x +a \right )^{3} \sin \left (b x +a \right )^{3}}{54 b^{6} \left (d x +c \right )^{2}}-\frac {5 \left (12 d^{4} x^{4} b^{4}+48 b^{4} c \,d^{3} x^{3}+72 b^{4} c^{2} d^{2} x^{2}+48 b^{4} c^{3} d x +12 b^{4} c^{4}-2 b^{2} d^{4} x^{2}-4 b^{2} c \,d^{3} x -2 b^{2} c^{2} d^{2}-21 d^{4}\right ) \left (3 \left (d x +c \right )^{2} \cos \left (b x +a \right )^{3} \sin \left (b x +a \right )^{3} d -3 \left (d x +c \right )^{3} \cos \left (b x +a \right )^{2} \sin \left (b x +a \right )^{4} b +3 \left (d x +c \right )^{3} \cos \left (b x +a \right )^{4} \sin \left (b x +a \right )^{2} b \right )}{216 \left (d x +c \right )^{4} b^{6}}+\frac {d \left (18 x^{2} d^{2} b^{2}+36 b^{2} c d x +18 b^{2} c^{2}-25 d^{2}\right ) \left (6 \left (d x +c \right ) \cos \left (b x +a \right )^{3} \sin \left (b x +a \right )^{3} d^{2}-18 \left (d x +c \right )^{2} \cos \left (b x +a \right )^{2} \sin \left (b x +a \right )^{4} d b +18 \left (d x +c \right )^{2} \cos \left (b x +a \right )^{4} \sin \left (b x +a \right )^{2} d b +6 \left (d x +c \right )^{3} \cos \left (b x +a \right ) \sin \left (b x +a \right )^{5} b^{2}-24 \left (d x +c \right )^{3} \cos \left (b x +a \right )^{3} \sin \left (b x +a \right )^{3} b^{2}+6 \left (d x +c \right )^{3} \cos \left (b x +a \right )^{5} \sin \left (b x +a \right ) b^{2}\right )}{216 b^{6} \left (d x +c \right )^{3}}-\frac {\left (3 x^{2} d^{2} b^{2}+6 b^{2} c d x +3 b^{2} c^{2}-5 d^{2}\right ) \left (6 d^{3} \cos \left (b x +a \right )^{3} \sin \left (b x +a \right )^{3}-54 \left (d x +c \right ) \cos \left (b x +a \right )^{2} \sin \left (b x +a \right )^{4} d^{2} b +54 \left (d x +c \right ) \cos \left (b x +a \right )^{4} \sin \left (b x +a \right )^{2} d^{2} b +54 \left (d x +c \right )^{2} \cos \left (b x +a \right ) \sin \left (b x +a \right )^{5} d \,b^{2}-216 \left (d x +c \right )^{2} \cos \left (b x +a \right )^{3} \sin \left (b x +a \right )^{3} d \,b^{2}+54 \left (d x +c \right )^{2} \cos \left (b x +a \right )^{5} \sin \left (b x +a \right ) d \,b^{2}-6 \left (d x +c \right )^{3} b^{3} \sin \left (b x +a \right )^{6}+102 \left (d x +c \right )^{3} \cos \left (b x +a \right )^{2} \sin \left (b x +a \right )^{4} b^{3}-102 \left (d x +c \right )^{3} \cos \left (b x +a \right )^{4} \sin \left (b x +a \right )^{2} b^{3}+6 \left (d x +c \right )^{3} \cos \left (b x +a \right )^{6} b^{3}\right )}{432 \left (d x +c \right )^{2} b^{6}}\) \(815\)
derivativedivides \(\text {Expression too large to display}\) \(1132\)
default \(\text {Expression too large to display}\) \(1132\)

Input:

int((d*x+c)^3*cos(b*x+a)^3*sin(b*x+a)^3,x,method=_RETURNVERBOSE)
 

Output:

1/1152*(-54*((d*x+c)^2*b^2-3/2*d^2)*(d*x+c)*b*cos(2*b*x+2*a)+6*((d*x+c)^2* 
b^2-1/6*d^2)*(d*x+c)*b*cos(6*b*x+6*a)+81*d*((d*x+c)^2*b^2-1/2*d^2)*sin(2*b 
*x+2*a)-3*d*((d*x+c)^2*b^2-1/18*d^2)*sin(6*b*x+6*a)+48*b^3*c^3-80*c*d^2*b) 
/b^4
 

Fricas [B] (verification not implemented)

Leaf count of result is larger than twice the leaf count of optimal. 349 vs. \(2 (165) = 330\).

Time = 0.08 (sec) , antiderivative size = 349, normalized size of antiderivative = 1.93 \[ \int (c+d x)^3 \cos ^3(a+b x) \sin ^3(a+b x) \, dx=\frac {9 \, b^{3} d^{3} x^{3} + 27 \, b^{3} c d^{2} x^{2} + 6 \, {\left (6 \, b^{3} d^{3} x^{3} + 18 \, b^{3} c d^{2} x^{2} + 6 \, b^{3} c^{3} - b c d^{2} + {\left (18 \, b^{3} c^{2} d - b d^{3}\right )} x\right )} \cos \left (b x + a\right )^{6} - 9 \, {\left (6 \, b^{3} d^{3} x^{3} + 18 \, b^{3} c d^{2} x^{2} + 6 \, b^{3} c^{3} - b c d^{2} + {\left (18 \, b^{3} c^{2} d - b d^{3}\right )} x\right )} \cos \left (b x + a\right )^{4} + 27 \, {\left (b d^{3} x + b c d^{2}\right )} \cos \left (b x + a\right )^{2} + 3 \, {\left (9 \, b^{3} c^{2} d - 5 \, b d^{3}\right )} x - {\left ({\left (18 \, b^{2} d^{3} x^{2} + 36 \, b^{2} c d^{2} x + 18 \, b^{2} c^{2} d - d^{3}\right )} \cos \left (b x + a\right )^{5} - {\left (18 \, b^{2} d^{3} x^{2} + 36 \, b^{2} c d^{2} x + 18 \, b^{2} c^{2} d - d^{3}\right )} \cos \left (b x + a\right )^{3} - 3 \, {\left (9 \, b^{2} d^{3} x^{2} + 18 \, b^{2} c d^{2} x + 9 \, b^{2} c^{2} d - 5 \, d^{3}\right )} \cos \left (b x + a\right )\right )} \sin \left (b x + a\right )}{216 \, b^{4}} \] Input:

integrate((d*x+c)^3*cos(b*x+a)^3*sin(b*x+a)^3,x, algorithm="fricas")
 

Output:

1/216*(9*b^3*d^3*x^3 + 27*b^3*c*d^2*x^2 + 6*(6*b^3*d^3*x^3 + 18*b^3*c*d^2* 
x^2 + 6*b^3*c^3 - b*c*d^2 + (18*b^3*c^2*d - b*d^3)*x)*cos(b*x + a)^6 - 9*( 
6*b^3*d^3*x^3 + 18*b^3*c*d^2*x^2 + 6*b^3*c^3 - b*c*d^2 + (18*b^3*c^2*d - b 
*d^3)*x)*cos(b*x + a)^4 + 27*(b*d^3*x + b*c*d^2)*cos(b*x + a)^2 + 3*(9*b^3 
*c^2*d - 5*b*d^3)*x - ((18*b^2*d^3*x^2 + 36*b^2*c*d^2*x + 18*b^2*c^2*d - d 
^3)*cos(b*x + a)^5 - (18*b^2*d^3*x^2 + 36*b^2*c*d^2*x + 18*b^2*c^2*d - d^3 
)*cos(b*x + a)^3 - 3*(9*b^2*d^3*x^2 + 18*b^2*c*d^2*x + 9*b^2*c^2*d - 5*d^3 
)*cos(b*x + a))*sin(b*x + a))/b^4
 

Sympy [B] (verification not implemented)

Leaf count of result is larger than twice the leaf count of optimal. 857 vs. \(2 (178) = 356\).

Time = 1.12 (sec) , antiderivative size = 857, normalized size of antiderivative = 4.73 \[ \int (c+d x)^3 \cos ^3(a+b x) \sin ^3(a+b x) \, dx =\text {Too large to display} \] Input:

integrate((d*x+c)**3*cos(b*x+a)**3*sin(b*x+a)**3,x)
 

Output:

Piecewise((-c**3*sin(a + b*x)**2*cos(a + b*x)**4/(4*b) - c**3*cos(a + b*x) 
**6/(12*b) + c**2*d*x*sin(a + b*x)**6/(8*b) + 3*c**2*d*x*sin(a + b*x)**4*c 
os(a + b*x)**2/(8*b) - 3*c**2*d*x*sin(a + b*x)**2*cos(a + b*x)**4/(8*b) - 
c**2*d*x*cos(a + b*x)**6/(8*b) + c*d**2*x**2*sin(a + b*x)**6/(8*b) + 3*c*d 
**2*x**2*sin(a + b*x)**4*cos(a + b*x)**2/(8*b) - 3*c*d**2*x**2*sin(a + b*x 
)**2*cos(a + b*x)**4/(8*b) - c*d**2*x**2*cos(a + b*x)**6/(8*b) + d**3*x**3 
*sin(a + b*x)**6/(24*b) + d**3*x**3*sin(a + b*x)**4*cos(a + b*x)**2/(8*b) 
- d**3*x**3*sin(a + b*x)**2*cos(a + b*x)**4/(8*b) - d**3*x**3*cos(a + b*x) 
**6/(24*b) + c**2*d*sin(a + b*x)**5*cos(a + b*x)/(8*b**2) + c**2*d*sin(a + 
 b*x)**3*cos(a + b*x)**3/(3*b**2) + c**2*d*sin(a + b*x)*cos(a + b*x)**5/(8 
*b**2) + c*d**2*x*sin(a + b*x)**5*cos(a + b*x)/(4*b**2) + 2*c*d**2*x*sin(a 
 + b*x)**3*cos(a + b*x)**3/(3*b**2) + c*d**2*x*sin(a + b*x)*cos(a + b*x)** 
5/(4*b**2) + d**3*x**2*sin(a + b*x)**5*cos(a + b*x)/(8*b**2) + d**3*x**2*s 
in(a + b*x)**3*cos(a + b*x)**3/(3*b**2) + d**3*x**2*sin(a + b*x)*cos(a + b 
*x)**5/(8*b**2) - c*d**2*sin(a + b*x)**6/(24*b**3) + c*d**2*sin(a + b*x)** 
2*cos(a + b*x)**4/(6*b**3) + 7*c*d**2*cos(a + b*x)**6/(72*b**3) - 5*d**3*x 
*sin(a + b*x)**6/(72*b**3) - d**3*x*sin(a + b*x)**4*cos(a + b*x)**2/(12*b* 
*3) + d**3*x*sin(a + b*x)**2*cos(a + b*x)**4/(12*b**3) + 5*d**3*x*cos(a + 
b*x)**6/(72*b**3) - 5*d**3*sin(a + b*x)**5*cos(a + b*x)/(72*b**4) - 31*d** 
3*sin(a + b*x)**3*cos(a + b*x)**3/(216*b**4) - 5*d**3*sin(a + b*x)*cos(...
 

Maxima [B] (verification not implemented)

Leaf count of result is larger than twice the leaf count of optimal. 602 vs. \(2 (165) = 330\).

Time = 0.06 (sec) , antiderivative size = 602, normalized size of antiderivative = 3.33 \[ \int (c+d x)^3 \cos ^3(a+b x) \sin ^3(a+b x) \, dx =\text {Too large to display} \] Input:

integrate((d*x+c)^3*cos(b*x+a)^3*sin(b*x+a)^3,x, algorithm="maxima")
 

Output:

-1/6912*(576*(2*sin(b*x + a)^6 - 3*sin(b*x + a)^4)*c^3 - 1728*(2*sin(b*x + 
 a)^6 - 3*sin(b*x + a)^4)*a*c^2*d/b + 1728*(2*sin(b*x + a)^6 - 3*sin(b*x + 
 a)^4)*a^2*c*d^2/b^2 - 576*(2*sin(b*x + a)^6 - 3*sin(b*x + a)^4)*a^3*d^3/b 
^3 - 18*(6*(b*x + a)*cos(6*b*x + 6*a) - 54*(b*x + a)*cos(2*b*x + 2*a) - si 
n(6*b*x + 6*a) + 27*sin(2*b*x + 2*a))*c^2*d/b + 36*(6*(b*x + a)*cos(6*b*x 
+ 6*a) - 54*(b*x + a)*cos(2*b*x + 2*a) - sin(6*b*x + 6*a) + 27*sin(2*b*x + 
 2*a))*a*c*d^2/b^2 - 18*(6*(b*x + a)*cos(6*b*x + 6*a) - 54*(b*x + a)*cos(2 
*b*x + 2*a) - sin(6*b*x + 6*a) + 27*sin(2*b*x + 2*a))*a^2*d^3/b^3 - 6*((18 
*(b*x + a)^2 - 1)*cos(6*b*x + 6*a) - 81*(2*(b*x + a)^2 - 1)*cos(2*b*x + 2* 
a) - 6*(b*x + a)*sin(6*b*x + 6*a) + 162*(b*x + a)*sin(2*b*x + 2*a))*c*d^2/ 
b^2 + 6*((18*(b*x + a)^2 - 1)*cos(6*b*x + 6*a) - 81*(2*(b*x + a)^2 - 1)*co 
s(2*b*x + 2*a) - 6*(b*x + a)*sin(6*b*x + 6*a) + 162*(b*x + a)*sin(2*b*x + 
2*a))*a*d^3/b^3 - (6*(6*(b*x + a)^3 - b*x - a)*cos(6*b*x + 6*a) - 162*(2*( 
b*x + a)^3 - 3*b*x - 3*a)*cos(2*b*x + 2*a) - (18*(b*x + a)^2 - 1)*sin(6*b* 
x + 6*a) + 243*(2*(b*x + a)^2 - 1)*sin(2*b*x + 2*a))*d^3/b^3)/b
 

Giac [A] (verification not implemented)

Time = 0.18 (sec) , antiderivative size = 241, normalized size of antiderivative = 1.33 \[ \int (c+d x)^3 \cos ^3(a+b x) \sin ^3(a+b x) \, dx=\frac {{\left (6 \, b^{3} d^{3} x^{3} + 18 \, b^{3} c d^{2} x^{2} + 18 \, b^{3} c^{2} d x + 6 \, b^{3} c^{3} - b d^{3} x - b c d^{2}\right )} \cos \left (6 \, b x + 6 \, a\right )}{1152 \, b^{4}} - \frac {3 \, {\left (2 \, b^{3} d^{3} x^{3} + 6 \, b^{3} c d^{2} x^{2} + 6 \, b^{3} c^{2} d x + 2 \, b^{3} c^{3} - 3 \, b d^{3} x - 3 \, b c d^{2}\right )} \cos \left (2 \, b x + 2 \, a\right )}{128 \, b^{4}} - \frac {{\left (18 \, b^{2} d^{3} x^{2} + 36 \, b^{2} c d^{2} x + 18 \, b^{2} c^{2} d - d^{3}\right )} \sin \left (6 \, b x + 6 \, a\right )}{6912 \, b^{4}} + \frac {9 \, {\left (2 \, b^{2} d^{3} x^{2} + 4 \, b^{2} c d^{2} x + 2 \, b^{2} c^{2} d - d^{3}\right )} \sin \left (2 \, b x + 2 \, a\right )}{256 \, b^{4}} \] Input:

integrate((d*x+c)^3*cos(b*x+a)^3*sin(b*x+a)^3,x, algorithm="giac")
 

Output:

1/1152*(6*b^3*d^3*x^3 + 18*b^3*c*d^2*x^2 + 18*b^3*c^2*d*x + 6*b^3*c^3 - b* 
d^3*x - b*c*d^2)*cos(6*b*x + 6*a)/b^4 - 3/128*(2*b^3*d^3*x^3 + 6*b^3*c*d^2 
*x^2 + 6*b^3*c^2*d*x + 2*b^3*c^3 - 3*b*d^3*x - 3*b*c*d^2)*cos(2*b*x + 2*a) 
/b^4 - 1/6912*(18*b^2*d^3*x^2 + 36*b^2*c*d^2*x + 18*b^2*c^2*d - d^3)*sin(6 
*b*x + 6*a)/b^4 + 9/256*(2*b^2*d^3*x^2 + 4*b^2*c*d^2*x + 2*b^2*c^2*d - d^3 
)*sin(2*b*x + 2*a)/b^4
 

Mupad [B] (verification not implemented)

Time = 1.34 (sec) , antiderivative size = 366, normalized size of antiderivative = 2.02 \[ \int (c+d x)^3 \cos ^3(a+b x) \sin ^3(a+b x) \, dx=-\frac {243\,d^3\,\sin \left (2\,a+2\,b\,x\right )-d^3\,\sin \left (6\,a+6\,b\,x\right )+324\,b^3\,c^3\,\cos \left (2\,a+2\,b\,x\right )-36\,b^3\,c^3\,\cos \left (6\,a+6\,b\,x\right )-486\,b^2\,c^2\,d\,\sin \left (2\,a+2\,b\,x\right )+18\,b^2\,c^2\,d\,\sin \left (6\,a+6\,b\,x\right )+324\,b^3\,d^3\,x^3\,\cos \left (2\,a+2\,b\,x\right )-36\,b^3\,d^3\,x^3\,\cos \left (6\,a+6\,b\,x\right )-486\,b^2\,d^3\,x^2\,\sin \left (2\,a+2\,b\,x\right )+18\,b^2\,d^3\,x^2\,\sin \left (6\,a+6\,b\,x\right )-486\,b\,c\,d^2\,\cos \left (2\,a+2\,b\,x\right )+6\,b\,c\,d^2\,\cos \left (6\,a+6\,b\,x\right )-486\,b\,d^3\,x\,\cos \left (2\,a+2\,b\,x\right )+6\,b\,d^3\,x\,\cos \left (6\,a+6\,b\,x\right )+972\,b^3\,c^2\,d\,x\,\cos \left (2\,a+2\,b\,x\right )-108\,b^3\,c^2\,d\,x\,\cos \left (6\,a+6\,b\,x\right )-972\,b^2\,c\,d^2\,x\,\sin \left (2\,a+2\,b\,x\right )+36\,b^2\,c\,d^2\,x\,\sin \left (6\,a+6\,b\,x\right )+972\,b^3\,c\,d^2\,x^2\,\cos \left (2\,a+2\,b\,x\right )-108\,b^3\,c\,d^2\,x^2\,\cos \left (6\,a+6\,b\,x\right )}{6912\,b^4} \] Input:

int(cos(a + b*x)^3*sin(a + b*x)^3*(c + d*x)^3,x)
 

Output:

-(243*d^3*sin(2*a + 2*b*x) - d^3*sin(6*a + 6*b*x) + 324*b^3*c^3*cos(2*a + 
2*b*x) - 36*b^3*c^3*cos(6*a + 6*b*x) - 486*b^2*c^2*d*sin(2*a + 2*b*x) + 18 
*b^2*c^2*d*sin(6*a + 6*b*x) + 324*b^3*d^3*x^3*cos(2*a + 2*b*x) - 36*b^3*d^ 
3*x^3*cos(6*a + 6*b*x) - 486*b^2*d^3*x^2*sin(2*a + 2*b*x) + 18*b^2*d^3*x^2 
*sin(6*a + 6*b*x) - 486*b*c*d^2*cos(2*a + 2*b*x) + 6*b*c*d^2*cos(6*a + 6*b 
*x) - 486*b*d^3*x*cos(2*a + 2*b*x) + 6*b*d^3*x*cos(6*a + 6*b*x) + 972*b^3* 
c^2*d*x*cos(2*a + 2*b*x) - 108*b^3*c^2*d*x*cos(6*a + 6*b*x) - 972*b^2*c*d^ 
2*x*sin(2*a + 2*b*x) + 36*b^2*c*d^2*x*sin(6*a + 6*b*x) + 972*b^3*c*d^2*x^2 
*cos(2*a + 2*b*x) - 108*b^3*c*d^2*x^2*cos(6*a + 6*b*x))/(6912*b^4)
 

Reduce [B] (verification not implemented)

Time = 0.15 (sec) , antiderivative size = 569, normalized size of antiderivative = 3.14 \[ \int (c+d x)^3 \cos ^3(a+b x) \sin ^3(a+b x) \, dx=\frac {27 \cos \left (b x +a \right ) \sin \left (b x +a \right ) b^{2} c^{2} d +27 \cos \left (b x +a \right ) \sin \left (b x +a \right ) b^{2} d^{3} x^{2}-36 \sin \left (b x +a \right )^{6} b^{3} c^{3}-27 \sin \left (b x +a \right )^{2} b \,d^{3} x +36 \cos \left (b x +a \right ) \sin \left (b x +a \right )^{3} b^{2} c \,d^{2} x -27 b^{3} c \,d^{2} x^{2}-18 \cos \left (b x +a \right ) \sin \left (b x +a \right )^{5} b^{2} c^{2} d -18 \cos \left (b x +a \right ) \sin \left (b x +a \right )^{5} b^{2} d^{3} x^{2}-108 \sin \left (b x +a \right )^{6} b^{3} c^{2} d x -108 \sin \left (b x +a \right )^{6} b^{3} c \,d^{2} x^{2}+18 \cos \left (b x +a \right ) \sin \left (b x +a \right )^{3} b^{2} c^{2} d +18 \cos \left (b x +a \right ) \sin \left (b x +a \right )^{3} b^{2} d^{3} x^{2}+162 \sin \left (b x +a \right )^{4} b^{3} c^{2} d x +162 \sin \left (b x +a \right )^{4} b^{3} c \,d^{2} x^{2}-\cos \left (b x +a \right ) \sin \left (b x +a \right )^{3} d^{3}+54 \sin \left (b x +a \right )^{4} b^{3} c^{3}-9 b^{3} d^{3} x^{3}+18 b c \,d^{2}+54 \cos \left (b x +a \right ) \sin \left (b x +a \right ) b^{2} c \,d^{2} x +\cos \left (b x +a \right ) \sin \left (b x +a \right )^{5} d^{3}-36 \sin \left (b x +a \right )^{6} b^{3} d^{3} x^{3}+6 \sin \left (b x +a \right )^{6} b c \,d^{2}+6 \sin \left (b x +a \right )^{6} b \,d^{3} x -15 \cos \left (b x +a \right ) \sin \left (b x +a \right ) d^{3}+15 b \,d^{3} x -36 \cos \left (b x +a \right ) \sin \left (b x +a \right )^{5} b^{2} c \,d^{2} x -3 a \,d^{3}+54 \sin \left (b x +a \right )^{4} b^{3} d^{3} x^{3}-9 \sin \left (b x +a \right )^{4} b c \,d^{2}-9 \sin \left (b x +a \right )^{4} b \,d^{3} x -27 \sin \left (b x +a \right )^{2} b c \,d^{2}-27 a \,b^{2} c^{2} d -27 b^{3} c^{2} d x}{216 b^{4}} \] Input:

int((d*x+c)^3*cos(b*x+a)^3*sin(b*x+a)^3,x)
 

Output:

( - 18*cos(a + b*x)*sin(a + b*x)**5*b**2*c**2*d - 36*cos(a + b*x)*sin(a + 
b*x)**5*b**2*c*d**2*x - 18*cos(a + b*x)*sin(a + b*x)**5*b**2*d**3*x**2 + c 
os(a + b*x)*sin(a + b*x)**5*d**3 + 18*cos(a + b*x)*sin(a + b*x)**3*b**2*c* 
*2*d + 36*cos(a + b*x)*sin(a + b*x)**3*b**2*c*d**2*x + 18*cos(a + b*x)*sin 
(a + b*x)**3*b**2*d**3*x**2 - cos(a + b*x)*sin(a + b*x)**3*d**3 + 27*cos(a 
 + b*x)*sin(a + b*x)*b**2*c**2*d + 54*cos(a + b*x)*sin(a + b*x)*b**2*c*d** 
2*x + 27*cos(a + b*x)*sin(a + b*x)*b**2*d**3*x**2 - 15*cos(a + b*x)*sin(a 
+ b*x)*d**3 - 36*sin(a + b*x)**6*b**3*c**3 - 108*sin(a + b*x)**6*b**3*c**2 
*d*x - 108*sin(a + b*x)**6*b**3*c*d**2*x**2 - 36*sin(a + b*x)**6*b**3*d**3 
*x**3 + 6*sin(a + b*x)**6*b*c*d**2 + 6*sin(a + b*x)**6*b*d**3*x + 54*sin(a 
 + b*x)**4*b**3*c**3 + 162*sin(a + b*x)**4*b**3*c**2*d*x + 162*sin(a + b*x 
)**4*b**3*c*d**2*x**2 + 54*sin(a + b*x)**4*b**3*d**3*x**3 - 9*sin(a + b*x) 
**4*b*c*d**2 - 9*sin(a + b*x)**4*b*d**3*x - 27*sin(a + b*x)**2*b*c*d**2 - 
27*sin(a + b*x)**2*b*d**3*x - 27*a*b**2*c**2*d - 3*a*d**3 - 27*b**3*c**2*d 
*x - 27*b**3*c*d**2*x**2 - 9*b**3*d**3*x**3 + 18*b*c*d**2 + 15*b*d**3*x)/( 
216*b**4)