\(\int (c+d x)^3 (a+b \sin (e+f x))^2 \, dx\) [157]

Optimal result
Mathematica [A] (verified)
Rubi [A] (verified)
Maple [A] (verified)
Fricas [A] (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 = 20, antiderivative size = 237 \[ \int (c+d x)^3 (a+b \sin (e+f x))^2 \, dx=-\frac {3 b^2 d (c+d x)^2}{8 f^2}+\frac {a^2 (c+d x)^4}{4 d}+\frac {b^2 (c+d x)^4}{8 d}+\frac {12 a b d^2 (c+d x) \cos (e+f x)}{f^3}-\frac {2 a b (c+d x)^3 \cos (e+f x)}{f}-\frac {12 a b d^3 \sin (e+f x)}{f^4}+\frac {6 a b d (c+d x)^2 \sin (e+f x)}{f^2}+\frac {3 b^2 d^2 (c+d x) \cos (e+f x) \sin (e+f x)}{4 f^3}-\frac {b^2 (c+d x)^3 \cos (e+f x) \sin (e+f x)}{2 f}-\frac {3 b^2 d^3 \sin ^2(e+f x)}{8 f^4}+\frac {3 b^2 d (c+d x)^2 \sin ^2(e+f x)}{4 f^2} \] Output:

-3/8*b^2*d*(d*x+c)^2/f^2+1/4*a^2*(d*x+c)^4/d+1/8*b^2*(d*x+c)^4/d+12*a*b*d^ 
2*(d*x+c)*cos(f*x+e)/f^3-2*a*b*(d*x+c)^3*cos(f*x+e)/f-12*a*b*d^3*sin(f*x+e 
)/f^4+6*a*b*d*(d*x+c)^2*sin(f*x+e)/f^2+3/4*b^2*d^2*(d*x+c)*cos(f*x+e)*sin( 
f*x+e)/f^3-1/2*b^2*(d*x+c)^3*cos(f*x+e)*sin(f*x+e)/f-3/8*b^2*d^3*sin(f*x+e 
)^2/f^4+3/4*b^2*d*(d*x+c)^2*sin(f*x+e)^2/f^2
 

Mathematica [A] (verified)

Time = 1.17 (sec) , antiderivative size = 232, normalized size of antiderivative = 0.98 \[ \int (c+d x)^3 (a+b \sin (e+f x))^2 \, dx=\frac {2 \left (2 a^2+b^2\right ) f^4 x \left (4 c^3+6 c^2 d x+4 c d^2 x^2+d^3 x^3\right )-32 a b f (c+d x) \left (c^2 f^2+2 c d f^2 x+d^2 \left (-6+f^2 x^2\right )\right ) \cos (e+f x)-3 b^2 d \left (2 c^2 f^2+4 c d f^2 x+d^2 \left (-1+2 f^2 x^2\right )\right ) \cos (2 (e+f x))+96 a b d \left (c^2 f^2+2 c d f^2 x+d^2 \left (-2+f^2 x^2\right )\right ) \sin (e+f x)-2 b^2 f (c+d x) \left (2 c^2 f^2+4 c d f^2 x+d^2 \left (-3+2 f^2 x^2\right )\right ) \sin (2 (e+f x))}{16 f^4} \] Input:

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

Output:

(2*(2*a^2 + b^2)*f^4*x*(4*c^3 + 6*c^2*d*x + 4*c*d^2*x^2 + d^3*x^3) - 32*a* 
b*f*(c + d*x)*(c^2*f^2 + 2*c*d*f^2*x + d^2*(-6 + f^2*x^2))*Cos[e + f*x] - 
3*b^2*d*(2*c^2*f^2 + 4*c*d*f^2*x + d^2*(-1 + 2*f^2*x^2))*Cos[2*(e + f*x)] 
+ 96*a*b*d*(c^2*f^2 + 2*c*d*f^2*x + d^2*(-2 + f^2*x^2))*Sin[e + f*x] - 2*b 
^2*f*(c + d*x)*(2*c^2*f^2 + 4*c*d*f^2*x + d^2*(-3 + 2*f^2*x^2))*Sin[2*(e + 
 f*x)])/(16*f^4)
 

Rubi [A] (verified)

Time = 0.52 (sec) , antiderivative size = 237, normalized size of antiderivative = 1.00, number of steps used = 3, number of rules used = 3, \(\frac {\text {number of rules}}{\text {integrand size}}\) = 0.150, Rules used = {3042, 3798, 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 (a+b \sin (e+f x))^2 \, dx\)

\(\Big \downarrow \) 3042

\(\displaystyle \int (c+d x)^3 (a+b \sin (e+f x))^2dx\)

\(\Big \downarrow \) 3798

\(\displaystyle \int \left (a^2 (c+d x)^3+2 a b (c+d x)^3 \sin (e+f x)+b^2 (c+d x)^3 \sin ^2(e+f x)\right )dx\)

\(\Big \downarrow \) 2009

\(\displaystyle \frac {a^2 (c+d x)^4}{4 d}+\frac {12 a b d^2 (c+d x) \cos (e+f x)}{f^3}+\frac {6 a b d (c+d x)^2 \sin (e+f x)}{f^2}-\frac {2 a b (c+d x)^3 \cos (e+f x)}{f}-\frac {12 a b d^3 \sin (e+f x)}{f^4}+\frac {3 b^2 d^2 (c+d x) \sin (e+f x) \cos (e+f x)}{4 f^3}+\frac {3 b^2 d (c+d x)^2 \sin ^2(e+f x)}{4 f^2}-\frac {b^2 (c+d x)^3 \sin (e+f x) \cos (e+f x)}{2 f}-\frac {3 b^2 d (c+d x)^2}{8 f^2}+\frac {b^2 (c+d x)^4}{8 d}-\frac {3 b^2 d^3 \sin ^2(e+f x)}{8 f^4}\)

Input:

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

Output:

(-3*b^2*d*(c + d*x)^2)/(8*f^2) + (a^2*(c + d*x)^4)/(4*d) + (b^2*(c + d*x)^ 
4)/(8*d) + (12*a*b*d^2*(c + d*x)*Cos[e + f*x])/f^3 - (2*a*b*(c + d*x)^3*Co 
s[e + f*x])/f - (12*a*b*d^3*Sin[e + f*x])/f^4 + (6*a*b*d*(c + d*x)^2*Sin[e 
 + f*x])/f^2 + (3*b^2*d^2*(c + d*x)*Cos[e + f*x]*Sin[e + f*x])/(4*f^3) - ( 
b^2*(c + d*x)^3*Cos[e + f*x]*Sin[e + f*x])/(2*f) - (3*b^2*d^3*Sin[e + f*x] 
^2)/(8*f^4) + (3*b^2*d*(c + d*x)^2*Sin[e + f*x]^2)/(4*f^2)
 

Defintions of rubi rules used

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

rule 3042
Int[u_, x_Symbol] :> Int[DeactivateTrig[u, x], x] /; FunctionOfTrigOfLinear 
Q[u, x]
 

rule 3798
Int[((c_.) + (d_.)*(x_))^(m_.)*((a_) + (b_.)*sin[(e_.) + (f_.)*(x_)])^(n_.) 
, x_Symbol] :> Int[ExpandIntegrand[(c + d*x)^m, (a + b*Sin[e + f*x])^n, x], 
 x] /; FreeQ[{a, b, c, d, e, f, m}, x] && IGtQ[n, 0] && (EqQ[n, 1] || IGtQ[ 
m, 0] || NeQ[a^2 - b^2, 0])
 
Maple [A] (verified)

Time = 2.20 (sec) , antiderivative size = 213, normalized size of antiderivative = 0.90

method result size
parallelrisch \(\frac {-4 \left (\left (d x +c \right )^{2} f^{2}-\frac {3 d^{2}}{2}\right ) \left (d x +c \right ) f \,b^{2} \sin \left (2 f x +2 e \right )-6 d \,b^{2} \left (\left (d x +c \right )^{2} f^{2}-\frac {d^{2}}{2}\right ) \cos \left (2 f x +2 e \right )-32 \left (d x +c \right ) a f b \left (\left (d x +c \right )^{2} f^{2}-6 d^{2}\right ) \cos \left (f x +e \right )+96 \left (\left (d x +c \right )^{2} f^{2}-2 d^{2}\right ) a d b \sin \left (f x +e \right )+16 x \left (\frac {d x}{2}+c \right ) \left (a^{2}+\frac {b^{2}}{2}\right ) \left (\frac {1}{2} x^{2} d^{2}+c d x +c^{2}\right ) f^{4}-32 a b \,c^{3} f^{3}+6 b^{2} c^{2} d \,f^{2}+192 a b c \,d^{2} f -3 b^{2} d^{3}}{16 f^{4}}\) \(213\)
risch \(\frac {d^{3} a^{2} x^{4}}{4}+\frac {d^{3} b^{2} x^{4}}{8}+d^{2} a^{2} c \,x^{3}+\frac {d^{2} b^{2} c \,x^{3}}{2}+\frac {3 d \,a^{2} c^{2} x^{2}}{2}+\frac {3 d \,b^{2} c^{2} x^{2}}{4}+a^{2} c^{3} x +\frac {b^{2} c^{3} x}{2}+\frac {a^{2} c^{4}}{4 d}+\frac {b^{2} c^{4}}{8 d}-\frac {2 a b \left (d^{3} f^{2} x^{3}+3 c \,d^{2} f^{2} x^{2}+3 c^{2} d \,f^{2} x +c^{3} f^{2}-6 d^{3} x -6 c \,d^{2}\right ) \cos \left (f x +e \right )}{f^{3}}+\frac {6 a b d \left (d^{2} x^{2} f^{2}+2 c d \,f^{2} x +c^{2} f^{2}-2 d^{2}\right ) \sin \left (f x +e \right )}{f^{4}}-\frac {3 b^{2} d \left (2 d^{2} x^{2} f^{2}+4 c d \,f^{2} x +2 c^{2} f^{2}-d^{2}\right ) \cos \left (2 f x +2 e \right )}{16 f^{4}}-\frac {b^{2} \left (2 d^{3} f^{2} x^{3}+6 c \,d^{2} f^{2} x^{2}+6 c^{2} d \,f^{2} x +2 c^{3} f^{2}-3 d^{3} x -3 c \,d^{2}\right ) \sin \left (2 f x +2 e \right )}{8 f^{3}}\) \(342\)
parts \(\text {Expression too large to display}\) \(916\)
norman \(\frac {\left (\frac {1}{4} d^{3} a^{2}+\frac {1}{8} b^{2} d^{3}\right ) x^{4}+\left (\frac {1}{2} d^{3} a^{2}+\frac {1}{4} b^{2} d^{3}\right ) x^{4} \tan \left (\frac {f x}{2}+\frac {e}{2}\right )^{2}+\left (\frac {1}{4} d^{3} a^{2}+\frac {1}{8} b^{2} d^{3}\right ) x^{4} \tan \left (\frac {f x}{2}+\frac {e}{2}\right )^{4}+\frac {b^{2} d^{3} x^{3} \tan \left (\frac {f x}{2}+\frac {e}{2}\right )^{3}}{f}+d^{2} c \left (2 a^{2}+b^{2}\right ) x^{3} \tan \left (\frac {f x}{2}+\frac {e}{2}\right )^{2}-\frac {8 a b \,c^{3} f^{3}+6 b^{2} c^{2} d \,f^{2}-48 a b c \,d^{2} f -3 b^{2} d^{3}}{4 f^{4}}+\frac {\left (4 a^{2} c^{3} f^{3}+2 b^{2} c^{3} f^{3}-24 a b \,c^{2} d \,f^{2}-3 b^{2} c \,d^{2} f +48 a b \,d^{3}\right ) x}{4 f^{3}}+\frac {\left (8 a b \,c^{3} f^{3}-6 b^{2} c^{2} d \,f^{2}-48 a b c \,d^{2} f +3 b^{2} d^{3}\right ) \tan \left (\frac {f x}{2}+\frac {e}{2}\right )^{4}}{4 f^{4}}+\frac {b \left (-2 b \,c^{3} f^{3}+24 a \,c^{2} d \,f^{2}+3 b c \,d^{2} f -48 a \,d^{3}\right ) \tan \left (\frac {f x}{2}+\frac {e}{2}\right )}{2 f^{4}}+\frac {b \left (2 b \,c^{3} f^{3}+24 a \,c^{2} d \,f^{2}-3 b c \,d^{2} f -48 a \,d^{3}\right ) \tan \left (\frac {f x}{2}+\frac {e}{2}\right )^{3}}{2 f^{4}}+\frac {3 d \left (4 a^{2} c^{2} f^{2}+2 b^{2} c^{2} f^{2}-16 a b c d f -b^{2} d^{2}\right ) x^{2}}{8 f^{2}}+\frac {d^{2} \left (2 a^{2} c f +b^{2} c f -4 a b d \right ) x^{3}}{2 f}+\frac {\left (4 a^{2} c^{3} f^{3}+2 b^{2} c^{3} f^{3}+24 a b \,c^{2} d \,f^{2}-3 b^{2} c \,d^{2} f -48 a b \,d^{3}\right ) x \tan \left (\frac {f x}{2}+\frac {e}{2}\right )^{4}}{4 f^{3}}-\frac {b^{2} d^{3} x^{3} \tan \left (\frac {f x}{2}+\frac {e}{2}\right )}{f}+\frac {c \left (4 a^{2} c^{2} f^{2}+2 b^{2} c^{2} f^{2}+9 b^{2} d^{2}\right ) x \tan \left (\frac {f x}{2}+\frac {e}{2}\right )^{2}}{2 f^{2}}+\frac {3 d \left (4 a^{2} c^{2} f^{2}+2 b^{2} c^{2} f^{2}+3 b^{2} d^{2}\right ) x^{2} \tan \left (\frac {f x}{2}+\frac {e}{2}\right )^{2}}{4 f^{2}}+\frac {3 d \left (4 a^{2} c^{2} f^{2}+2 b^{2} c^{2} f^{2}+16 a b c d f -b^{2} d^{2}\right ) x^{2} \tan \left (\frac {f x}{2}+\frac {e}{2}\right )^{4}}{8 f^{2}}+\frac {d^{2} \left (2 a^{2} c f +b^{2} c f +4 a b d \right ) x^{3} \tan \left (\frac {f x}{2}+\frac {e}{2}\right )^{4}}{2 f}+\frac {3 b d \left (2 b \,c^{2} f^{2}+16 a c d f -d^{2} b \right ) x \tan \left (\frac {f x}{2}+\frac {e}{2}\right )^{3}}{2 f^{3}}+\frac {3 d b \left (-2 b \,c^{2} f^{2}+16 a c d f +d^{2} b \right ) x \tan \left (\frac {f x}{2}+\frac {e}{2}\right )}{2 f^{3}}+\frac {3 d^{2} b \left (-b c f +4 a d \right ) x^{2} \tan \left (\frac {f x}{2}+\frac {e}{2}\right )}{f^{2}}+\frac {3 d^{2} b \left (b c f +4 a d \right ) x^{2} \tan \left (\frac {f x}{2}+\frac {e}{2}\right )^{3}}{f^{2}}}{\left (1+\tan \left (\frac {f x}{2}+\frac {e}{2}\right )^{2}\right )^{2}}\) \(922\)
derivativedivides \(\text {Expression too large to display}\) \(1125\)
default \(\text {Expression too large to display}\) \(1125\)
orering \(\text {Expression too large to display}\) \(1317\)

Input:

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

Output:

1/16*(-4*((d*x+c)^2*f^2-3/2*d^2)*(d*x+c)*f*b^2*sin(2*f*x+2*e)-6*d*b^2*((d* 
x+c)^2*f^2-1/2*d^2)*cos(2*f*x+2*e)-32*(d*x+c)*a*f*b*((d*x+c)^2*f^2-6*d^2)* 
cos(f*x+e)+96*((d*x+c)^2*f^2-2*d^2)*a*d*b*sin(f*x+e)+16*x*(1/2*d*x+c)*(a^2 
+1/2*b^2)*(1/2*x^2*d^2+c*d*x+c^2)*f^4-32*a*b*c^3*f^3+6*b^2*c^2*d*f^2+192*a 
*b*c*d^2*f-3*b^2*d^3)/f^4
 

Fricas [A] (verification not implemented)

Time = 0.09 (sec) , antiderivative size = 382, normalized size of antiderivative = 1.61 \[ \int (c+d x)^3 (a+b \sin (e+f x))^2 \, dx=\frac {{\left (2 \, a^{2} + b^{2}\right )} d^{3} f^{4} x^{4} + 4 \, {\left (2 \, a^{2} + b^{2}\right )} c d^{2} f^{4} x^{3} + 3 \, {\left (2 \, {\left (2 \, a^{2} + b^{2}\right )} c^{2} d f^{4} + b^{2} d^{3} f^{2}\right )} x^{2} - 3 \, {\left (2 \, b^{2} d^{3} f^{2} x^{2} + 4 \, b^{2} c d^{2} f^{2} x + 2 \, b^{2} c^{2} d f^{2} - b^{2} d^{3}\right )} \cos \left (f x + e\right )^{2} + 2 \, {\left (2 \, {\left (2 \, a^{2} + b^{2}\right )} c^{3} f^{4} + 3 \, b^{2} c d^{2} f^{2}\right )} x - 16 \, {\left (a b d^{3} f^{3} x^{3} + 3 \, a b c d^{2} f^{3} x^{2} + a b c^{3} f^{3} - 6 \, a b c d^{2} f + 3 \, {\left (a b c^{2} d f^{3} - 2 \, a b d^{3} f\right )} x\right )} \cos \left (f x + e\right ) + 2 \, {\left (24 \, a b d^{3} f^{2} x^{2} + 48 \, a b c d^{2} f^{2} x + 24 \, a b c^{2} d f^{2} - 48 \, a b d^{3} - {\left (2 \, b^{2} d^{3} f^{3} x^{3} + 6 \, b^{2} c d^{2} f^{3} x^{2} + 2 \, b^{2} c^{3} f^{3} - 3 \, b^{2} c d^{2} f + 3 \, {\left (2 \, b^{2} c^{2} d f^{3} - b^{2} d^{3} f\right )} x\right )} \cos \left (f x + e\right )\right )} \sin \left (f x + e\right )}{8 \, f^{4}} \] Input:

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

Output:

1/8*((2*a^2 + b^2)*d^3*f^4*x^4 + 4*(2*a^2 + b^2)*c*d^2*f^4*x^3 + 3*(2*(2*a 
^2 + b^2)*c^2*d*f^4 + b^2*d^3*f^2)*x^2 - 3*(2*b^2*d^3*f^2*x^2 + 4*b^2*c*d^ 
2*f^2*x + 2*b^2*c^2*d*f^2 - b^2*d^3)*cos(f*x + e)^2 + 2*(2*(2*a^2 + b^2)*c 
^3*f^4 + 3*b^2*c*d^2*f^2)*x - 16*(a*b*d^3*f^3*x^3 + 3*a*b*c*d^2*f^3*x^2 + 
a*b*c^3*f^3 - 6*a*b*c*d^2*f + 3*(a*b*c^2*d*f^3 - 2*a*b*d^3*f)*x)*cos(f*x + 
 e) + 2*(24*a*b*d^3*f^2*x^2 + 48*a*b*c*d^2*f^2*x + 24*a*b*c^2*d*f^2 - 48*a 
*b*d^3 - (2*b^2*d^3*f^3*x^3 + 6*b^2*c*d^2*f^3*x^2 + 2*b^2*c^3*f^3 - 3*b^2* 
c*d^2*f + 3*(2*b^2*c^2*d*f^3 - b^2*d^3*f)*x)*cos(f*x + e))*sin(f*x + e))/f 
^4
 

Sympy [B] (verification not implemented)

Leaf count of result is larger than twice the leaf count of optimal. 779 vs. \(2 (240) = 480\).

Time = 0.42 (sec) , antiderivative size = 779, normalized size of antiderivative = 3.29 \[ \int (c+d x)^3 (a+b \sin (e+f x))^2 \, dx =\text {Too large to display} \] Input:

integrate((d*x+c)**3*(a+b*sin(f*x+e))**2,x)
 

Output:

Piecewise((a**2*c**3*x + 3*a**2*c**2*d*x**2/2 + a**2*c*d**2*x**3 + a**2*d* 
*3*x**4/4 - 2*a*b*c**3*cos(e + f*x)/f - 6*a*b*c**2*d*x*cos(e + f*x)/f + 6* 
a*b*c**2*d*sin(e + f*x)/f**2 - 6*a*b*c*d**2*x**2*cos(e + f*x)/f + 12*a*b*c 
*d**2*x*sin(e + f*x)/f**2 + 12*a*b*c*d**2*cos(e + f*x)/f**3 - 2*a*b*d**3*x 
**3*cos(e + f*x)/f + 6*a*b*d**3*x**2*sin(e + f*x)/f**2 + 12*a*b*d**3*x*cos 
(e + f*x)/f**3 - 12*a*b*d**3*sin(e + f*x)/f**4 + b**2*c**3*x*sin(e + f*x)* 
*2/2 + b**2*c**3*x*cos(e + f*x)**2/2 - b**2*c**3*sin(e + f*x)*cos(e + f*x) 
/(2*f) + 3*b**2*c**2*d*x**2*sin(e + f*x)**2/4 + 3*b**2*c**2*d*x**2*cos(e + 
 f*x)**2/4 - 3*b**2*c**2*d*x*sin(e + f*x)*cos(e + f*x)/(2*f) + 3*b**2*c**2 
*d*sin(e + f*x)**2/(4*f**2) + b**2*c*d**2*x**3*sin(e + f*x)**2/2 + b**2*c* 
d**2*x**3*cos(e + f*x)**2/2 - 3*b**2*c*d**2*x**2*sin(e + f*x)*cos(e + f*x) 
/(2*f) + 3*b**2*c*d**2*x*sin(e + f*x)**2/(4*f**2) - 3*b**2*c*d**2*x*cos(e 
+ f*x)**2/(4*f**2) + 3*b**2*c*d**2*sin(e + f*x)*cos(e + f*x)/(4*f**3) + b* 
*2*d**3*x**4*sin(e + f*x)**2/8 + b**2*d**3*x**4*cos(e + f*x)**2/8 - b**2*d 
**3*x**3*sin(e + f*x)*cos(e + f*x)/(2*f) + 3*b**2*d**3*x**2*sin(e + f*x)** 
2/(8*f**2) - 3*b**2*d**3*x**2*cos(e + f*x)**2/(8*f**2) + 3*b**2*d**3*x*sin 
(e + f*x)*cos(e + f*x)/(4*f**3) - 3*b**2*d**3*sin(e + f*x)**2/(8*f**4), Ne 
(f, 0)), ((a + b*sin(e))**2*(c**3*x + 3*c**2*d*x**2/2 + c*d**2*x**3 + d**3 
*x**4/4), True))
 

Maxima [B] (verification not implemented)

Leaf count of result is larger than twice the leaf count of optimal. 959 vs. \(2 (223) = 446\).

Time = 0.08 (sec) , antiderivative size = 959, normalized size of antiderivative = 4.05 \[ \int (c+d x)^3 (a+b \sin (e+f x))^2 \, dx=\text {Too large to display} \] Input:

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

Output:

1/16*(16*(f*x + e)*a^2*c^3 + 4*(2*f*x + 2*e - sin(2*f*x + 2*e))*b^2*c^3 + 
4*(f*x + e)^4*a^2*d^3/f^3 - 16*(f*x + e)^3*a^2*d^3*e/f^3 + 24*(f*x + e)^2* 
a^2*d^3*e^2/f^3 - 16*(f*x + e)*a^2*d^3*e^3/f^3 - 4*(2*f*x + 2*e - sin(2*f* 
x + 2*e))*b^2*d^3*e^3/f^3 + 16*(f*x + e)^3*a^2*c*d^2/f^2 - 48*(f*x + e)^2* 
a^2*c*d^2*e/f^2 + 48*(f*x + e)*a^2*c*d^2*e^2/f^2 + 12*(2*f*x + 2*e - sin(2 
*f*x + 2*e))*b^2*c*d^2*e^2/f^2 + 24*(f*x + e)^2*a^2*c^2*d/f - 48*(f*x + e) 
*a^2*c^2*d*e/f - 12*(2*f*x + 2*e - sin(2*f*x + 2*e))*b^2*c^2*d*e/f - 32*a* 
b*c^3*cos(f*x + e) + 32*a*b*d^3*e^3*cos(f*x + e)/f^3 - 96*a*b*c*d^2*e^2*co 
s(f*x + e)/f^2 + 96*a*b*c^2*d*e*cos(f*x + e)/f - 96*((f*x + e)*cos(f*x + e 
) - sin(f*x + e))*a*b*d^3*e^2/f^3 + 6*(2*(f*x + e)^2 - 2*(f*x + e)*sin(2*f 
*x + 2*e) - cos(2*f*x + 2*e))*b^2*d^3*e^2/f^3 + 192*((f*x + e)*cos(f*x + e 
) - sin(f*x + e))*a*b*c*d^2*e/f^2 - 12*(2*(f*x + e)^2 - 2*(f*x + e)*sin(2* 
f*x + 2*e) - cos(2*f*x + 2*e))*b^2*c*d^2*e/f^2 - 96*((f*x + e)*cos(f*x + e 
) - sin(f*x + e))*a*b*c^2*d/f + 6*(2*(f*x + e)^2 - 2*(f*x + e)*sin(2*f*x + 
 2*e) - cos(2*f*x + 2*e))*b^2*c^2*d/f + 96*(((f*x + e)^2 - 2)*cos(f*x + e) 
 - 2*(f*x + e)*sin(f*x + e))*a*b*d^3*e/f^3 - 2*(4*(f*x + e)^3 - 6*(f*x + e 
)*cos(2*f*x + 2*e) - 3*(2*(f*x + e)^2 - 1)*sin(2*f*x + 2*e))*b^2*d^3*e/f^3 
 - 96*(((f*x + e)^2 - 2)*cos(f*x + e) - 2*(f*x + e)*sin(f*x + e))*a*b*c*d^ 
2/f^2 + 2*(4*(f*x + e)^3 - 6*(f*x + e)*cos(2*f*x + 2*e) - 3*(2*(f*x + e)^2 
 - 1)*sin(2*f*x + 2*e))*b^2*c*d^2/f^2 - 32*(((f*x + e)^3 - 6*f*x - 6*e)...
 

Giac [A] (verification not implemented)

Time = 0.78 (sec) , antiderivative size = 367, normalized size of antiderivative = 1.55 \[ \int (c+d x)^3 (a+b \sin (e+f x))^2 \, dx=\frac {1}{4} \, a^{2} d^{3} x^{4} + \frac {1}{8} \, b^{2} d^{3} x^{4} + a^{2} c d^{2} x^{3} + \frac {1}{2} \, b^{2} c d^{2} x^{3} + \frac {3}{2} \, a^{2} c^{2} d x^{2} + \frac {3}{4} \, b^{2} c^{2} d x^{2} + a^{2} c^{3} x + \frac {1}{2} \, b^{2} c^{3} x - \frac {3 \, {\left (2 \, b^{2} d^{3} f^{2} x^{2} + 4 \, b^{2} c d^{2} f^{2} x + 2 \, b^{2} c^{2} d f^{2} - b^{2} d^{3}\right )} \cos \left (2 \, f x + 2 \, e\right )}{16 \, f^{4}} - \frac {2 \, {\left (a b d^{3} f^{3} x^{3} + 3 \, a b c d^{2} f^{3} x^{2} + 3 \, a b c^{2} d f^{3} x + a b c^{3} f^{3} - 6 \, a b d^{3} f x - 6 \, a b c d^{2} f\right )} \cos \left (f x + e\right )}{f^{4}} - \frac {{\left (2 \, b^{2} d^{3} f^{3} x^{3} + 6 \, b^{2} c d^{2} f^{3} x^{2} + 6 \, b^{2} c^{2} d f^{3} x + 2 \, b^{2} c^{3} f^{3} - 3 \, b^{2} d^{3} f x - 3 \, b^{2} c d^{2} f\right )} \sin \left (2 \, f x + 2 \, e\right )}{8 \, f^{4}} + \frac {6 \, {\left (a b d^{3} f^{2} x^{2} + 2 \, a b c d^{2} f^{2} x + a b c^{2} d f^{2} - 2 \, a b d^{3}\right )} \sin \left (f x + e\right )}{f^{4}} \] Input:

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

Output:

1/4*a^2*d^3*x^4 + 1/8*b^2*d^3*x^4 + a^2*c*d^2*x^3 + 1/2*b^2*c*d^2*x^3 + 3/ 
2*a^2*c^2*d*x^2 + 3/4*b^2*c^2*d*x^2 + a^2*c^3*x + 1/2*b^2*c^3*x - 3/16*(2* 
b^2*d^3*f^2*x^2 + 4*b^2*c*d^2*f^2*x + 2*b^2*c^2*d*f^2 - b^2*d^3)*cos(2*f*x 
 + 2*e)/f^4 - 2*(a*b*d^3*f^3*x^3 + 3*a*b*c*d^2*f^3*x^2 + 3*a*b*c^2*d*f^3*x 
 + a*b*c^3*f^3 - 6*a*b*d^3*f*x - 6*a*b*c*d^2*f)*cos(f*x + e)/f^4 - 1/8*(2* 
b^2*d^3*f^3*x^3 + 6*b^2*c*d^2*f^3*x^2 + 6*b^2*c^2*d*f^3*x + 2*b^2*c^3*f^3 
- 3*b^2*d^3*f*x - 3*b^2*c*d^2*f)*sin(2*f*x + 2*e)/f^4 + 6*(a*b*d^3*f^2*x^2 
 + 2*a*b*c*d^2*f^2*x + a*b*c^2*d*f^2 - 2*a*b*d^3)*sin(f*x + e)/f^4
 

Mupad [B] (verification not implemented)

Time = 38.10 (sec) , antiderivative size = 497, normalized size of antiderivative = 2.10 \[ \int (c+d x)^3 (a+b \sin (e+f x))^2 \, dx=\frac {\frac {3\,b^2\,d^3\,\cos \left (2\,e+2\,f\,x\right )}{2}+8\,a^2\,c^3\,f^4\,x+4\,b^2\,c^3\,f^4\,x-96\,a\,b\,d^3\,\sin \left (e+f\,x\right )-2\,b^2\,c^3\,f^3\,\sin \left (2\,e+2\,f\,x\right )+2\,a^2\,d^3\,f^4\,x^4+b^2\,d^3\,f^4\,x^4-16\,a\,b\,c^3\,f^3\,\cos \left (e+f\,x\right )-3\,b^2\,d^3\,f^2\,x^2\,\cos \left (2\,e+2\,f\,x\right )-2\,b^2\,d^3\,f^3\,x^3\,\sin \left (2\,e+2\,f\,x\right )+3\,b^2\,c\,d^2\,f\,\sin \left (2\,e+2\,f\,x\right )+3\,b^2\,d^3\,f\,x\,\sin \left (2\,e+2\,f\,x\right )-3\,b^2\,c^2\,d\,f^2\,\cos \left (2\,e+2\,f\,x\right )+12\,a^2\,c^2\,d\,f^4\,x^2+8\,a^2\,c\,d^2\,f^4\,x^3+6\,b^2\,c^2\,d\,f^4\,x^2+4\,b^2\,c\,d^2\,f^4\,x^3+96\,a\,b\,c\,d^2\,f\,\cos \left (e+f\,x\right )+96\,a\,b\,d^3\,f\,x\,\cos \left (e+f\,x\right )-6\,b^2\,c\,d^2\,f^2\,x\,\cos \left (2\,e+2\,f\,x\right )-6\,b^2\,c^2\,d\,f^3\,x\,\sin \left (2\,e+2\,f\,x\right )+48\,a\,b\,c^2\,d\,f^2\,\sin \left (e+f\,x\right )-6\,b^2\,c\,d^2\,f^3\,x^2\,\sin \left (2\,e+2\,f\,x\right )-16\,a\,b\,d^3\,f^3\,x^3\,\cos \left (e+f\,x\right )+48\,a\,b\,d^3\,f^2\,x^2\,\sin \left (e+f\,x\right )-48\,a\,b\,c\,d^2\,f^3\,x^2\,\cos \left (e+f\,x\right )-48\,a\,b\,c^2\,d\,f^3\,x\,\cos \left (e+f\,x\right )+96\,a\,b\,c\,d^2\,f^2\,x\,\sin \left (e+f\,x\right )}{8\,f^4} \] Input:

int((a + b*sin(e + f*x))^2*(c + d*x)^3,x)
 

Output:

((3*b^2*d^3*cos(2*e + 2*f*x))/2 + 8*a^2*c^3*f^4*x + 4*b^2*c^3*f^4*x - 96*a 
*b*d^3*sin(e + f*x) - 2*b^2*c^3*f^3*sin(2*e + 2*f*x) + 2*a^2*d^3*f^4*x^4 + 
 b^2*d^3*f^4*x^4 - 16*a*b*c^3*f^3*cos(e + f*x) - 3*b^2*d^3*f^2*x^2*cos(2*e 
 + 2*f*x) - 2*b^2*d^3*f^3*x^3*sin(2*e + 2*f*x) + 3*b^2*c*d^2*f*sin(2*e + 2 
*f*x) + 3*b^2*d^3*f*x*sin(2*e + 2*f*x) - 3*b^2*c^2*d*f^2*cos(2*e + 2*f*x) 
+ 12*a^2*c^2*d*f^4*x^2 + 8*a^2*c*d^2*f^4*x^3 + 6*b^2*c^2*d*f^4*x^2 + 4*b^2 
*c*d^2*f^4*x^3 + 96*a*b*c*d^2*f*cos(e + f*x) + 96*a*b*d^3*f*x*cos(e + f*x) 
 - 6*b^2*c*d^2*f^2*x*cos(2*e + 2*f*x) - 6*b^2*c^2*d*f^3*x*sin(2*e + 2*f*x) 
 + 48*a*b*c^2*d*f^2*sin(e + f*x) - 6*b^2*c*d^2*f^3*x^2*sin(2*e + 2*f*x) - 
16*a*b*d^3*f^3*x^3*cos(e + f*x) + 48*a*b*d^3*f^2*x^2*sin(e + f*x) - 48*a*b 
*c*d^2*f^3*x^2*cos(e + f*x) - 48*a*b*c^2*d*f^3*x*cos(e + f*x) + 96*a*b*c*d 
^2*f^2*x*sin(e + f*x))/(8*f^4)
 

Reduce [B] (verification not implemented)

Time = 0.37 (sec) , antiderivative size = 558, normalized size of antiderivative = 2.35 \[ \int (c+d x)^3 (a+b \sin (e+f x))^2 \, dx=\frac {6 \cos \left (f x +e \right ) \sin \left (f x +e \right ) b^{2} d^{3} f x +96 \cos \left (f x +e \right ) a b c \,d^{2} f -16 \cos \left (f x +e \right ) a b \,d^{3} f^{3} x^{3}+96 \cos \left (f x +e \right ) a b \,d^{3} f x +12 \sin \left (f x +e \right )^{2} b^{2} c \,d^{2} f^{2} x +48 \sin \left (f x +e \right ) a b \,c^{2} d \,f^{2}+48 \sin \left (f x +e \right ) a b \,d^{3} f^{2} x^{2}+6 \cos \left (f x +e \right ) \sin \left (f x +e \right ) b^{2} c \,d^{2} f -4 \cos \left (f x +e \right ) \sin \left (f x +e \right ) b^{2} d^{3} f^{3} x^{3}-3 \sin \left (f x +e \right )^{2} b^{2} d^{3}+4 b^{2} c \,d^{2} f^{4} x^{3}-6 b^{2} c \,d^{2} f^{2} x +6 b^{2} d^{3}+b^{2} d^{3} f^{4} x^{4}-96 \sin \left (f x +e \right ) a b \,d^{3}+8 a^{2} c^{3} f^{4} x +2 a^{2} d^{3} f^{4} x^{4}+4 b^{2} c^{3} f^{4} x -12 b^{2} c^{2} d \,f^{2}-3 b^{2} d^{3} f^{2} x^{2}-4 \cos \left (f x +e \right ) \sin \left (f x +e \right ) b^{2} c^{3} f^{3}-16 \cos \left (f x +e \right ) a b \,c^{3} f^{3}+6 \sin \left (f x +e \right )^{2} b^{2} c^{2} d \,f^{2}+6 \sin \left (f x +e \right )^{2} b^{2} d^{3} f^{2} x^{2}+12 a^{2} c^{2} d \,f^{4} x^{2}+8 a^{2} c \,d^{2} f^{4} x^{3}+6 b^{2} c^{2} d \,f^{4} x^{2}-12 \cos \left (f x +e \right ) \sin \left (f x +e \right ) b^{2} c^{2} d \,f^{3} x -12 \cos \left (f x +e \right ) \sin \left (f x +e \right ) b^{2} c \,d^{2} f^{3} x^{2}-48 \cos \left (f x +e \right ) a b \,c^{2} d \,f^{3} x -48 \cos \left (f x +e \right ) a b c \,d^{2} f^{3} x^{2}+96 \sin \left (f x +e \right ) a b c \,d^{2} f^{2} x}{8 f^{4}} \] Input:

int((d*x+c)^3*(a+b*sin(f*x+e))^2,x)
 

Output:

( - 4*cos(e + f*x)*sin(e + f*x)*b**2*c**3*f**3 - 12*cos(e + f*x)*sin(e + f 
*x)*b**2*c**2*d*f**3*x - 12*cos(e + f*x)*sin(e + f*x)*b**2*c*d**2*f**3*x** 
2 + 6*cos(e + f*x)*sin(e + f*x)*b**2*c*d**2*f - 4*cos(e + f*x)*sin(e + f*x 
)*b**2*d**3*f**3*x**3 + 6*cos(e + f*x)*sin(e + f*x)*b**2*d**3*f*x - 16*cos 
(e + f*x)*a*b*c**3*f**3 - 48*cos(e + f*x)*a*b*c**2*d*f**3*x - 48*cos(e + f 
*x)*a*b*c*d**2*f**3*x**2 + 96*cos(e + f*x)*a*b*c*d**2*f - 16*cos(e + f*x)* 
a*b*d**3*f**3*x**3 + 96*cos(e + f*x)*a*b*d**3*f*x + 6*sin(e + f*x)**2*b**2 
*c**2*d*f**2 + 12*sin(e + f*x)**2*b**2*c*d**2*f**2*x + 6*sin(e + f*x)**2*b 
**2*d**3*f**2*x**2 - 3*sin(e + f*x)**2*b**2*d**3 + 48*sin(e + f*x)*a*b*c** 
2*d*f**2 + 96*sin(e + f*x)*a*b*c*d**2*f**2*x + 48*sin(e + f*x)*a*b*d**3*f* 
*2*x**2 - 96*sin(e + f*x)*a*b*d**3 + 8*a**2*c**3*f**4*x + 12*a**2*c**2*d*f 
**4*x**2 + 8*a**2*c*d**2*f**4*x**3 + 2*a**2*d**3*f**4*x**4 + 4*b**2*c**3*f 
**4*x + 6*b**2*c**2*d*f**4*x**2 - 12*b**2*c**2*d*f**2 + 4*b**2*c*d**2*f**4 
*x**3 - 6*b**2*c*d**2*f**2*x + b**2*d**3*f**4*x**4 - 3*b**2*d**3*f**2*x**2 
 + 6*b**2*d**3)/(8*f**4)