\(\int (a+a \sin (e+f x))^{7/2} (c-c \sin (e+f x))^{7/2} \, dx\) [378]

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

Optimal result

Integrand size = 30, antiderivative size = 179 \[ \int (a+a \sin (e+f x))^{7/2} (c-c \sin (e+f x))^{7/2} \, dx=-\frac {2 a^4 \cos (e+f x) (c-c \sin (e+f x))^{7/2}}{35 f \sqrt {a+a \sin (e+f x)}}-\frac {4 a^3 \cos (e+f x) \sqrt {a+a \sin (e+f x)} (c-c \sin (e+f x))^{7/2}}{35 f}-\frac {a^2 \cos (e+f x) (a+a \sin (e+f x))^{3/2} (c-c \sin (e+f x))^{7/2}}{7 f}-\frac {a \cos (e+f x) (a+a \sin (e+f x))^{5/2} (c-c \sin (e+f x))^{7/2}}{7 f} \] Output:

-2/35*a^4*cos(f*x+e)*(c-c*sin(f*x+e))^(7/2)/f/(a+a*sin(f*x+e))^(1/2)-4/35* 
a^3*cos(f*x+e)*(a+a*sin(f*x+e))^(1/2)*(c-c*sin(f*x+e))^(7/2)/f-1/7*a^2*cos 
(f*x+e)*(a+a*sin(f*x+e))^(3/2)*(c-c*sin(f*x+e))^(7/2)/f-1/7*a*cos(f*x+e)*( 
a+a*sin(f*x+e))^(5/2)*(c-c*sin(f*x+e))^(7/2)/f
 

Mathematica [A] (verified)

Time = 2.11 (sec) , antiderivative size = 80, normalized size of antiderivative = 0.45 \[ \int (a+a \sin (e+f x))^{7/2} (c-c \sin (e+f x))^{7/2} \, dx=-\frac {a^3 c^3 \sqrt {a (1+\sin (e+f x))} \sqrt {c-c \sin (e+f x)} \left (-35+35 \sin ^2(e+f x)-21 \sin ^4(e+f x)+5 \sin ^6(e+f x)\right ) \tan (e+f x)}{35 f} \] Input:

Integrate[(a + a*Sin[e + f*x])^(7/2)*(c - c*Sin[e + f*x])^(7/2),x]
 

Output:

-1/35*(a^3*c^3*Sqrt[a*(1 + Sin[e + f*x])]*Sqrt[c - c*Sin[e + f*x]]*(-35 + 
35*Sin[e + f*x]^2 - 21*Sin[e + f*x]^4 + 5*Sin[e + f*x]^6)*Tan[e + f*x])/f
 

Rubi [A] (verified)

Time = 0.85 (sec) , antiderivative size = 187, normalized size of antiderivative = 1.04, number of steps used = 8, number of rules used = 8, \(\frac {\text {number of rules}}{\text {integrand size}}\) = 0.267, Rules used = {3042, 3219, 3042, 3219, 3042, 3219, 3042, 3217}

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 (a \sin (e+f x)+a)^{7/2} (c-c \sin (e+f x))^{7/2} \, dx\)

\(\Big \downarrow \) 3042

\(\displaystyle \int (a \sin (e+f x)+a)^{7/2} (c-c \sin (e+f x))^{7/2}dx\)

\(\Big \downarrow \) 3219

\(\displaystyle \frac {6}{7} a \int (\sin (e+f x) a+a)^{5/2} (c-c \sin (e+f x))^{7/2}dx-\frac {a \cos (e+f x) (a \sin (e+f x)+a)^{5/2} (c-c \sin (e+f x))^{7/2}}{7 f}\)

\(\Big \downarrow \) 3042

\(\displaystyle \frac {6}{7} a \int (\sin (e+f x) a+a)^{5/2} (c-c \sin (e+f x))^{7/2}dx-\frac {a \cos (e+f x) (a \sin (e+f x)+a)^{5/2} (c-c \sin (e+f x))^{7/2}}{7 f}\)

\(\Big \downarrow \) 3219

\(\displaystyle \frac {6}{7} a \left (\frac {2}{3} a \int (\sin (e+f x) a+a)^{3/2} (c-c \sin (e+f x))^{7/2}dx-\frac {a \cos (e+f x) (a \sin (e+f x)+a)^{3/2} (c-c \sin (e+f x))^{7/2}}{6 f}\right )-\frac {a \cos (e+f x) (a \sin (e+f x)+a)^{5/2} (c-c \sin (e+f x))^{7/2}}{7 f}\)

\(\Big \downarrow \) 3042

\(\displaystyle \frac {6}{7} a \left (\frac {2}{3} a \int (\sin (e+f x) a+a)^{3/2} (c-c \sin (e+f x))^{7/2}dx-\frac {a \cos (e+f x) (a \sin (e+f x)+a)^{3/2} (c-c \sin (e+f x))^{7/2}}{6 f}\right )-\frac {a \cos (e+f x) (a \sin (e+f x)+a)^{5/2} (c-c \sin (e+f x))^{7/2}}{7 f}\)

\(\Big \downarrow \) 3219

\(\displaystyle \frac {6}{7} a \left (\frac {2}{3} a \left (\frac {2}{5} a \int \sqrt {\sin (e+f x) a+a} (c-c \sin (e+f x))^{7/2}dx-\frac {a \cos (e+f x) \sqrt {a \sin (e+f x)+a} (c-c \sin (e+f x))^{7/2}}{5 f}\right )-\frac {a \cos (e+f x) (a \sin (e+f x)+a)^{3/2} (c-c \sin (e+f x))^{7/2}}{6 f}\right )-\frac {a \cos (e+f x) (a \sin (e+f x)+a)^{5/2} (c-c \sin (e+f x))^{7/2}}{7 f}\)

\(\Big \downarrow \) 3042

\(\displaystyle \frac {6}{7} a \left (\frac {2}{3} a \left (\frac {2}{5} a \int \sqrt {\sin (e+f x) a+a} (c-c \sin (e+f x))^{7/2}dx-\frac {a \cos (e+f x) \sqrt {a \sin (e+f x)+a} (c-c \sin (e+f x))^{7/2}}{5 f}\right )-\frac {a \cos (e+f x) (a \sin (e+f x)+a)^{3/2} (c-c \sin (e+f x))^{7/2}}{6 f}\right )-\frac {a \cos (e+f x) (a \sin (e+f x)+a)^{5/2} (c-c \sin (e+f x))^{7/2}}{7 f}\)

\(\Big \downarrow \) 3217

\(\displaystyle \frac {6}{7} a \left (\frac {2}{3} a \left (-\frac {a^2 \cos (e+f x) (c-c \sin (e+f x))^{7/2}}{10 f \sqrt {a \sin (e+f x)+a}}-\frac {a \cos (e+f x) \sqrt {a \sin (e+f x)+a} (c-c \sin (e+f x))^{7/2}}{5 f}\right )-\frac {a \cos (e+f x) (a \sin (e+f x)+a)^{3/2} (c-c \sin (e+f x))^{7/2}}{6 f}\right )-\frac {a \cos (e+f x) (a \sin (e+f x)+a)^{5/2} (c-c \sin (e+f x))^{7/2}}{7 f}\)

Input:

Int[(a + a*Sin[e + f*x])^(7/2)*(c - c*Sin[e + f*x])^(7/2),x]
 

Output:

-1/7*(a*Cos[e + f*x]*(a + a*Sin[e + f*x])^(5/2)*(c - c*Sin[e + f*x])^(7/2) 
)/f + (6*a*(-1/6*(a*Cos[e + f*x]*(a + a*Sin[e + f*x])^(3/2)*(c - c*Sin[e + 
 f*x])^(7/2))/f + (2*a*(-1/10*(a^2*Cos[e + f*x]*(c - c*Sin[e + f*x])^(7/2) 
)/(f*Sqrt[a + a*Sin[e + f*x]]) - (a*Cos[e + f*x]*Sqrt[a + a*Sin[e + f*x]]* 
(c - c*Sin[e + f*x])^(7/2))/(5*f)))/3))/7
 

Defintions of rubi rules used

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

rule 3217
Int[Sqrt[(a_) + (b_.)*sin[(e_.) + (f_.)*(x_)]]*((c_) + (d_.)*sin[(e_.) + (f 
_.)*(x_)])^(n_), x_Symbol] :> Simp[-2*b*Cos[e + f*x]*((c + d*Sin[e + f*x])^ 
n/(f*(2*n + 1)*Sqrt[a + b*Sin[e + f*x]])), x] /; FreeQ[{a, b, c, d, e, f, n 
}, x] && EqQ[b*c + a*d, 0] && EqQ[a^2 - b^2, 0] && NeQ[n, -2^(-1)]
 

rule 3219
Int[((a_) + (b_.)*sin[(e_.) + (f_.)*(x_)])^(m_)*((c_) + (d_.)*sin[(e_.) + ( 
f_.)*(x_)])^(n_), x_Symbol] :> Simp[(-b)*Cos[e + f*x]*(a + b*Sin[e + f*x])^ 
(m - 1)*((c + d*Sin[e + f*x])^n/(f*(m + n))), x] + Simp[a*((2*m - 1)/(m + n 
))   Int[(a + b*Sin[e + f*x])^(m - 1)*(c + d*Sin[e + f*x])^n, x], x] /; Fre 
eQ[{a, b, c, d, e, f, n}, x] && EqQ[b*c + a*d, 0] && EqQ[a^2 - b^2, 0] && I 
GtQ[m - 1/2, 0] &&  !LtQ[n, -1] &&  !(IGtQ[n - 1/2, 0] && LtQ[n, m]) &&  !( 
ILtQ[m + n, 0] && GtQ[2*m + n + 1, 0])
 
Maple [A] (verified)

Time = 3.32 (sec) , antiderivative size = 127, normalized size of antiderivative = 0.71

method result size
default \(\frac {16 \sin \left (\frac {\pi }{4}+\frac {f x}{2}+\frac {e}{2}\right )^{6} \tan \left (\frac {\pi }{4}+\frac {f x}{2}+\frac {e}{2}\right ) \sqrt {a \sin \left (\frac {\pi }{4}+\frac {f x}{2}+\frac {e}{2}\right )^{2}}\, \sqrt {c \cos \left (\frac {\pi }{4}+\frac {f x}{2}+\frac {e}{2}\right )^{2}}\, \left (20 \cos \left (\frac {\pi }{4}+\frac {f x}{2}+\frac {e}{2}\right )^{6}+10 \cos \left (\frac {\pi }{4}+\frac {f x}{2}+\frac {e}{2}\right )^{4}+4 \cos \left (\frac {\pi }{4}+\frac {f x}{2}+\frac {e}{2}\right )^{2}+1\right ) a^{3} c^{3} \sqrt {4}}{35 f}\) \(127\)

Input:

int((a+sin(f*x+e)*a)^(7/2)*(c-c*sin(f*x+e))^(7/2),x,method=_RETURNVERBOSE)
 

Output:

16/35/f*sin(1/4*Pi+1/2*f*x+1/2*e)^6*tan(1/4*Pi+1/2*f*x+1/2*e)*(a*sin(1/4*P 
i+1/2*f*x+1/2*e)^2)^(1/2)*(c*cos(1/4*Pi+1/2*f*x+1/2*e)^2)^(1/2)*(20*cos(1/ 
4*Pi+1/2*f*x+1/2*e)^6+10*cos(1/4*Pi+1/2*f*x+1/2*e)^4+4*cos(1/4*Pi+1/2*f*x+ 
1/2*e)^2+1)*a^3*c^3*4^(1/2)
 

Fricas [A] (verification not implemented)

Time = 0.10 (sec) , antiderivative size = 101, normalized size of antiderivative = 0.56 \[ \int (a+a \sin (e+f x))^{7/2} (c-c \sin (e+f x))^{7/2} \, dx=\frac {{\left (5 \, a^{3} c^{3} \cos \left (f x + e\right )^{6} + 6 \, a^{3} c^{3} \cos \left (f x + e\right )^{4} + 8 \, a^{3} c^{3} \cos \left (f x + e\right )^{2} + 16 \, a^{3} c^{3}\right )} \sqrt {a \sin \left (f x + e\right ) + a} \sqrt {-c \sin \left (f x + e\right ) + c} \sin \left (f x + e\right )}{35 \, f \cos \left (f x + e\right )} \] Input:

integrate((a+a*sin(f*x+e))^(7/2)*(c-c*sin(f*x+e))^(7/2),x, algorithm="fric 
as")
 

Output:

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

Sympy [F(-1)]

Timed out. \[ \int (a+a \sin (e+f x))^{7/2} (c-c \sin (e+f x))^{7/2} \, dx=\text {Timed out} \] Input:

integrate((a+a*sin(f*x+e))**(7/2)*(c-c*sin(f*x+e))**(7/2),x)
 

Output:

Timed out
 

Maxima [F]

\[ \int (a+a \sin (e+f x))^{7/2} (c-c \sin (e+f x))^{7/2} \, dx=\int { {\left (a \sin \left (f x + e\right ) + a\right )}^{\frac {7}{2}} {\left (-c \sin \left (f x + e\right ) + c\right )}^{\frac {7}{2}} \,d x } \] Input:

integrate((a+a*sin(f*x+e))^(7/2)*(c-c*sin(f*x+e))^(7/2),x, algorithm="maxi 
ma")
 

Output:

integrate((a*sin(f*x + e) + a)^(7/2)*(-c*sin(f*x + e) + c)^(7/2), x)
 

Giac [A] (verification not implemented)

Time = 0.16 (sec) , antiderivative size = 204, normalized size of antiderivative = 1.14 \[ \int (a+a \sin (e+f x))^{7/2} (c-c \sin (e+f x))^{7/2} \, dx=-\frac {32 \, {\left (20 \, a^{3} c^{3} \mathrm {sgn}\left (\cos \left (-\frac {1}{4} \, \pi + \frac {1}{2} \, f x + \frac {1}{2} \, e\right )\right ) \mathrm {sgn}\left (\sin \left (-\frac {1}{4} \, \pi + \frac {1}{2} \, f x + \frac {1}{2} \, e\right )\right ) \sin \left (-\frac {1}{4} \, \pi + \frac {1}{2} \, f x + \frac {1}{2} \, e\right )^{14} - 70 \, a^{3} c^{3} \mathrm {sgn}\left (\cos \left (-\frac {1}{4} \, \pi + \frac {1}{2} \, f x + \frac {1}{2} \, e\right )\right ) \mathrm {sgn}\left (\sin \left (-\frac {1}{4} \, \pi + \frac {1}{2} \, f x + \frac {1}{2} \, e\right )\right ) \sin \left (-\frac {1}{4} \, \pi + \frac {1}{2} \, f x + \frac {1}{2} \, e\right )^{12} + 84 \, a^{3} c^{3} \mathrm {sgn}\left (\cos \left (-\frac {1}{4} \, \pi + \frac {1}{2} \, f x + \frac {1}{2} \, e\right )\right ) \mathrm {sgn}\left (\sin \left (-\frac {1}{4} \, \pi + \frac {1}{2} \, f x + \frac {1}{2} \, e\right )\right ) \sin \left (-\frac {1}{4} \, \pi + \frac {1}{2} \, f x + \frac {1}{2} \, e\right )^{10} - 35 \, a^{3} c^{3} \mathrm {sgn}\left (\cos \left (-\frac {1}{4} \, \pi + \frac {1}{2} \, f x + \frac {1}{2} \, e\right )\right ) \mathrm {sgn}\left (\sin \left (-\frac {1}{4} \, \pi + \frac {1}{2} \, f x + \frac {1}{2} \, e\right )\right ) \sin \left (-\frac {1}{4} \, \pi + \frac {1}{2} \, f x + \frac {1}{2} \, e\right )^{8}\right )} \sqrt {a} \sqrt {c}}{35 \, f} \] Input:

integrate((a+a*sin(f*x+e))^(7/2)*(c-c*sin(f*x+e))^(7/2),x, algorithm="giac 
")
 

Output:

-32/35*(20*a^3*c^3*sgn(cos(-1/4*pi + 1/2*f*x + 1/2*e))*sgn(sin(-1/4*pi + 1 
/2*f*x + 1/2*e))*sin(-1/4*pi + 1/2*f*x + 1/2*e)^14 - 70*a^3*c^3*sgn(cos(-1 
/4*pi + 1/2*f*x + 1/2*e))*sgn(sin(-1/4*pi + 1/2*f*x + 1/2*e))*sin(-1/4*pi 
+ 1/2*f*x + 1/2*e)^12 + 84*a^3*c^3*sgn(cos(-1/4*pi + 1/2*f*x + 1/2*e))*sgn 
(sin(-1/4*pi + 1/2*f*x + 1/2*e))*sin(-1/4*pi + 1/2*f*x + 1/2*e)^10 - 35*a^ 
3*c^3*sgn(cos(-1/4*pi + 1/2*f*x + 1/2*e))*sgn(sin(-1/4*pi + 1/2*f*x + 1/2* 
e))*sin(-1/4*pi + 1/2*f*x + 1/2*e)^8)*sqrt(a)*sqrt(c)/f
 

Mupad [B] (verification not implemented)

Time = 20.90 (sec) , antiderivative size = 179, normalized size of antiderivative = 1.00 \[ \int (a+a \sin (e+f x))^{7/2} (c-c \sin (e+f x))^{7/2} \, dx=\frac {\frac {1225\,a^3\,c^3\,\sin \left (e+f\,x\right )\,\sqrt {a+a\,\sin \left (e+f\,x\right )}\,\sqrt {c-c\,\sin \left (e+f\,x\right )}}{32}+\frac {245\,a^3\,c^3\,\sin \left (3\,e+3\,f\,x\right )\,\sqrt {a+a\,\sin \left (e+f\,x\right )}\,\sqrt {c-c\,\sin \left (e+f\,x\right )}}{32}+\frac {49\,a^3\,c^3\,\sin \left (5\,e+5\,f\,x\right )\,\sqrt {a+a\,\sin \left (e+f\,x\right )}\,\sqrt {c-c\,\sin \left (e+f\,x\right )}}{32}+\frac {5\,a^3\,c^3\,\sin \left (7\,e+7\,f\,x\right )\,\sqrt {a+a\,\sin \left (e+f\,x\right )}\,\sqrt {c-c\,\sin \left (e+f\,x\right )}}{32}}{70\,f\,\cos \left (e+f\,x\right )} \] Input:

int((a + a*sin(e + f*x))^(7/2)*(c - c*sin(e + f*x))^(7/2),x)
 

Output:

((1225*a^3*c^3*sin(e + f*x)*(a + a*sin(e + f*x))^(1/2)*(c - c*sin(e + f*x) 
)^(1/2))/32 + (245*a^3*c^3*sin(3*e + 3*f*x)*(a + a*sin(e + f*x))^(1/2)*(c 
- c*sin(e + f*x))^(1/2))/32 + (49*a^3*c^3*sin(5*e + 5*f*x)*(a + a*sin(e + 
f*x))^(1/2)*(c - c*sin(e + f*x))^(1/2))/32 + (5*a^3*c^3*sin(7*e + 7*f*x)*( 
a + a*sin(e + f*x))^(1/2)*(c - c*sin(e + f*x))^(1/2))/32)/(70*f*cos(e + f* 
x))
 

Reduce [F]

\[ \int (a+a \sin (e+f x))^{7/2} (c-c \sin (e+f x))^{7/2} \, dx=\sqrt {c}\, \sqrt {a}\, a^{3} c^{3} \left (-\left (\int \sqrt {\sin \left (f x +e \right )+1}\, \sqrt {-\sin \left (f x +e \right )+1}\, \sin \left (f x +e \right )^{6}d x \right )+3 \left (\int \sqrt {\sin \left (f x +e \right )+1}\, \sqrt {-\sin \left (f x +e \right )+1}\, \sin \left (f x +e \right )^{4}d x \right )-3 \left (\int \sqrt {\sin \left (f x +e \right )+1}\, \sqrt {-\sin \left (f x +e \right )+1}\, \sin \left (f x +e \right )^{2}d x \right )+\int \sqrt {\sin \left (f x +e \right )+1}\, \sqrt {-\sin \left (f x +e \right )+1}d x \right ) \] Input:

int((a+a*sin(f*x+e))^(7/2)*(c-c*sin(f*x+e))^(7/2),x)
 

Output:

sqrt(c)*sqrt(a)*a**3*c**3*( - int(sqrt(sin(e + f*x) + 1)*sqrt( - sin(e + f 
*x) + 1)*sin(e + f*x)**6,x) + 3*int(sqrt(sin(e + f*x) + 1)*sqrt( - sin(e + 
 f*x) + 1)*sin(e + f*x)**4,x) - 3*int(sqrt(sin(e + f*x) + 1)*sqrt( - sin(e 
 + f*x) + 1)*sin(e + f*x)**2,x) + int(sqrt(sin(e + f*x) + 1)*sqrt( - sin(e 
 + f*x) + 1),x))