\(\int \sec ^4(e+f x) (a+a \sin (e+f x))^m (A+B \sin (e+f x)) \, dx\) [1027]

Optimal result

Integrand size = 31, antiderivative size = 130 \[ \int \sec ^4(e+f x) (a+a \sin (e+f x))^m (A+B \sin (e+f x)) \, dx=\frac {B \sec ^3(e+f x) (a+a \sin (e+f x))^m}{f (3-m)}+\frac {2^{-\frac {3}{2}+m} (A (3-m)-B m) \operatorname {Hypergeometric2F1}\left (-\frac {3}{2},\frac {5}{2}-m,-\frac {1}{2},\frac {1}{2} (1-\sin (e+f x))\right ) \sec ^3(e+f x) (1+\sin (e+f x))^{\frac {3}{2}-m} (a+a \sin (e+f x))^m}{3 f (3-m)} \] Output:

B*sec(f*x+e)^3*(a+a*sin(f*x+e))^m/f/(3-m)+1/3*2^(-3/2+m)*(A*(3-m)-B*m)*hyp 
ergeom([-3/2, 5/2-m],[-1/2],1/2-1/2*sin(f*x+e))*sec(f*x+e)^3*(1+sin(f*x+e) 
)^(3/2-m)*(a+a*sin(f*x+e))^m/f/(3-m)
 

Mathematica [A] (verified)

Time = 1.16 (sec) , antiderivative size = 97, normalized size of antiderivative = 0.75 \[ \int \sec ^4(e+f x) (a+a \sin (e+f x))^m (A+B \sin (e+f x)) \, dx=\frac {\sec ^3(e+f x) (a (1+\sin (e+f x)))^m \left (-6 B+2^{-\frac {1}{2}+m} (A (-3+m)+B m) \operatorname {Hypergeometric2F1}\left (-\frac {3}{2},\frac {5}{2}-m,-\frac {1}{2},\frac {1}{2} (1-\sin (e+f x))\right ) (1+\sin (e+f x))^{\frac {3}{2}-m}\right )}{6 f (-3+m)} \] Input:

Integrate[Sec[e + f*x]^4*(a + a*Sin[e + f*x])^m*(A + B*Sin[e + f*x]),x]
 

Output:

(Sec[e + f*x]^3*(a*(1 + Sin[e + f*x]))^m*(-6*B + 2^(-1/2 + m)*(A*(-3 + m) 
+ B*m)*Hypergeometric2F1[-3/2, 5/2 - m, -1/2, (1 - Sin[e + f*x])/2]*(1 + S 
in[e + f*x])^(3/2 - m)))/(6*f*(-3 + m))
 

Rubi [A] (verified)

Time = 0.45 (sec) , antiderivative size = 129, normalized size of antiderivative = 0.99, number of steps used = 7, number of rules used = 6, \(\frac {\text {number of rules}}{\text {integrand size}}\) = 0.194, Rules used = {3042, 3339, 3042, 3168, 80, 79}

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 \sec ^4(e+f x) (a \sin (e+f x)+a)^m (A+B \sin (e+f x)) \, dx\)

\(\Big \downarrow \) 3042

\(\displaystyle \int \frac {(a \sin (e+f x)+a)^m (A+B \sin (e+f x))}{\cos (e+f x)^4}dx\)

\(\Big \downarrow \) 3339

\(\displaystyle \left (A-\frac {B m}{3-m}\right ) \int \sec ^4(e+f x) (\sin (e+f x) a+a)^mdx+\frac {B \sec ^3(e+f x) (a \sin (e+f x)+a)^m}{f (3-m)}\)

\(\Big \downarrow \) 3042

\(\displaystyle \left (A-\frac {B m}{3-m}\right ) \int \frac {(\sin (e+f x) a+a)^m}{\cos (e+f x)^4}dx+\frac {B \sec ^3(e+f x) (a \sin (e+f x)+a)^m}{f (3-m)}\)

\(\Big \downarrow \) 3168

\(\displaystyle \frac {a^2 \left (A-\frac {B m}{3-m}\right ) \sec ^3(e+f x) (a-a \sin (e+f x))^{3/2} (a \sin (e+f x)+a)^{3/2} \int \frac {(\sin (e+f x) a+a)^{m-\frac {5}{2}}}{(a-a \sin (e+f x))^{5/2}}d\sin (e+f x)}{f}+\frac {B \sec ^3(e+f x) (a \sin (e+f x)+a)^m}{f (3-m)}\)

\(\Big \downarrow \) 80

\(\displaystyle \frac {2^{m-\frac {5}{2}} \left (A-\frac {B m}{3-m}\right ) \sec ^3(e+f x) (a-a \sin (e+f x))^{3/2} (\sin (e+f x)+1)^{\frac {1}{2}-m} (a \sin (e+f x)+a)^{m+1} \int \frac {\left (\frac {1}{2} \sin (e+f x)+\frac {1}{2}\right )^{m-\frac {5}{2}}}{(a-a \sin (e+f x))^{5/2}}d\sin (e+f x)}{f}+\frac {B \sec ^3(e+f x) (a \sin (e+f x)+a)^m}{f (3-m)}\)

\(\Big \downarrow \) 79

\(\displaystyle \frac {2^{m-\frac {3}{2}} \left (A-\frac {B m}{3-m}\right ) \sec ^3(e+f x) (\sin (e+f x)+1)^{\frac {1}{2}-m} (a \sin (e+f x)+a)^{m+1} \operatorname {Hypergeometric2F1}\left (-\frac {3}{2},\frac {5}{2}-m,-\frac {1}{2},\frac {1}{2} (1-\sin (e+f x))\right )}{3 a f}+\frac {B \sec ^3(e+f x) (a \sin (e+f x)+a)^m}{f (3-m)}\)

Input:

Int[Sec[e + f*x]^4*(a + a*Sin[e + f*x])^m*(A + B*Sin[e + f*x]),x]
 

Output:

(B*Sec[e + f*x]^3*(a + a*Sin[e + f*x])^m)/(f*(3 - m)) + (2^(-3/2 + m)*(A - 
 (B*m)/(3 - m))*Hypergeometric2F1[-3/2, 5/2 - m, -1/2, (1 - Sin[e + f*x])/ 
2]*Sec[e + f*x]^3*(1 + Sin[e + f*x])^(1/2 - m)*(a + a*Sin[e + f*x])^(1 + m 
))/(3*a*f)
 

Defintions of rubi rules used

rule 79

Int[((a_) + (b_.)*(x_))^(m_)*((c_) + (d_.)*(x_))^(n_), x_Symbol] :> Simp[(( 
a + b*x)^(m + 1)/(b*(m + 1)*(b/(b*c - a*d))^n))*Hypergeometric2F1[-n, m + 1 
, m + 2, (-d)*((a + b*x)/(b*c - a*d))], x] /; FreeQ[{a, b, c, d, m, n}, x] 
&&  !IntegerQ[m] &&  !IntegerQ[n] && GtQ[b/(b*c - a*d), 0] && (RationalQ[m] 
 ||  !(RationalQ[n] && GtQ[-d/(b*c - a*d), 0]))
 

rule 80

Int[((a_) + (b_.)*(x_))^(m_)*((c_) + (d_.)*(x_))^(n_), x_Symbol] :> Simp[(c 
 + d*x)^FracPart[n]/((b/(b*c - a*d))^IntPart[n]*(b*((c + d*x)/(b*c - a*d))) 
^FracPart[n])   Int[(a + b*x)^m*Simp[b*(c/(b*c - a*d)) + b*d*(x/(b*c - a*d) 
), x]^n, x], x] /; FreeQ[{a, b, c, d, m, n}, x] &&  !IntegerQ[m] &&  !Integ 
erQ[n] && (RationalQ[m] ||  !SimplerQ[n + 1, m + 1])
 

rule 3042

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

rule 3168

Int[(cos[(e_.) + (f_.)*(x_)]*(g_.))^(p_)*((a_) + (b_.)*sin[(e_.) + (f_.)*(x 
_)])^(m_.), x_Symbol] :> Simp[a^2*((g*Cos[e + f*x])^(p + 1)/(f*g*(a + b*Sin 
[e + f*x])^((p + 1)/2)*(a - b*Sin[e + f*x])^((p + 1)/2)))   Subst[Int[(a + 
b*x)^(m + (p - 1)/2)*(a - b*x)^((p - 1)/2), x], x, Sin[e + f*x]], x] /; Fre 
eQ[{a, b, e, f, g, m, p}, x] && EqQ[a^2 - b^2, 0] &&  !IntegerQ[m]
 

rule 3339

Int[(cos[(e_.) + (f_.)*(x_)]*(g_.))^(p_)*((a_) + (b_.)*sin[(e_.) + (f_.)*(x 
_)])^(m_.)*((c_.) + (d_.)*sin[(e_.) + (f_.)*(x_)]), x_Symbol] :> Simp[(-d)* 
(g*Cos[e + f*x])^(p + 1)*((a + b*Sin[e + f*x])^m/(f*g*(m + p + 1))), x] + S 
imp[(a*d*m + b*c*(m + p + 1))/(b*(m + p + 1))   Int[(g*Cos[e + f*x])^p*(a + 
 b*Sin[e + f*x])^m, x], x] /; FreeQ[{a, b, c, d, e, f, g, m, p}, x] && EqQ[ 
a^2 - b^2, 0] && NeQ[m + p + 1, 0]
 

Maple [F]

\[\int \sec \left (f x +e \right )^{4} \left (a +a \sin \left (f x +e \right )\right )^{m} \left (A +B \sin \left (f x +e \right )\right )d x\]

Input:

int(sec(f*x+e)^4*(a+a*sin(f*x+e))^m*(A+B*sin(f*x+e)),x)
 

Output:

int(sec(f*x+e)^4*(a+a*sin(f*x+e))^m*(A+B*sin(f*x+e)),x)
 

Fricas [F]

\[ \int \sec ^4(e+f x) (a+a \sin (e+f x))^m (A+B \sin (e+f x)) \, dx=\int { {\left (B \sin \left (f x + e\right ) + A\right )} {\left (a \sin \left (f x + e\right ) + a\right )}^{m} \sec \left (f x + e\right )^{4} \,d x } \] Input:

integrate(sec(f*x+e)^4*(a+a*sin(f*x+e))^m*(A+B*sin(f*x+e)),x, algorithm="f 
ricas")
 

Output:

integral((B*sec(f*x + e)^4*sin(f*x + e) + A*sec(f*x + e)^4)*(a*sin(f*x + e 
) + a)^m, x)
 

Sympy [F(-2)]

Exception generated. \[ \int \sec ^4(e+f x) (a+a \sin (e+f x))^m (A+B \sin (e+f x)) \, dx=\text {Exception raised: HeuristicGCDFailed} \] Input:

integrate(sec(f*x+e)**4*(a+a*sin(f*x+e))**m*(A+B*sin(f*x+e)),x)
 

Output:

Exception raised: HeuristicGCDFailed >> no luck
 

Maxima [F]

\[ \int \sec ^4(e+f x) (a+a \sin (e+f x))^m (A+B \sin (e+f x)) \, dx=\int { {\left (B \sin \left (f x + e\right ) + A\right )} {\left (a \sin \left (f x + e\right ) + a\right )}^{m} \sec \left (f x + e\right )^{4} \,d x } \] Input:

integrate(sec(f*x+e)^4*(a+a*sin(f*x+e))^m*(A+B*sin(f*x+e)),x, algorithm="m 
axima")
 

Output:

integrate((B*sin(f*x + e) + A)*(a*sin(f*x + e) + a)^m*sec(f*x + e)^4, x)
 

Giac [F]

\[ \int \sec ^4(e+f x) (a+a \sin (e+f x))^m (A+B \sin (e+f x)) \, dx=\int { {\left (B \sin \left (f x + e\right ) + A\right )} {\left (a \sin \left (f x + e\right ) + a\right )}^{m} \sec \left (f x + e\right )^{4} \,d x } \] Input:

integrate(sec(f*x+e)^4*(a+a*sin(f*x+e))^m*(A+B*sin(f*x+e)),x, algorithm="g 
iac")
 

Output:

integrate((B*sin(f*x + e) + A)*(a*sin(f*x + e) + a)^m*sec(f*x + e)^4, x)
 

Mupad [F(-1)]

Timed out. \[ \int \sec ^4(e+f x) (a+a \sin (e+f x))^m (A+B \sin (e+f x)) \, dx=\int \frac {\left (A+B\,\sin \left (e+f\,x\right )\right )\,{\left (a+a\,\sin \left (e+f\,x\right )\right )}^m}{{\cos \left (e+f\,x\right )}^4} \,d x \] Input:

int(((A + B*sin(e + f*x))*(a + a*sin(e + f*x))^m)/cos(e + f*x)^4,x)
 

Output:

int(((A + B*sin(e + f*x))*(a + a*sin(e + f*x))^m)/cos(e + f*x)^4, x)
 

Reduce [F]

\[ \int \sec ^4(e+f x) (a+a \sin (e+f x))^m (A+B \sin (e+f x)) \, dx=\left (\int \left (a +a \sin \left (f x +e \right )\right )^{m} \sec \left (f x +e \right )^{4} \sin \left (f x +e \right )d x \right ) b +\left (\int \left (a +a \sin \left (f x +e \right )\right )^{m} \sec \left (f x +e \right )^{4}d x \right ) a \] Input:

int(sec(f*x+e)^4*(a+a*sin(f*x+e))^m*(A+B*sin(f*x+e)),x)
 

Output:

int((sin(e + f*x)*a + a)**m*sec(e + f*x)**4*sin(e + f*x),x)*b + int((sin(e 
 + f*x)*a + a)**m*sec(e + f*x)**4,x)*a