\(\int \frac {(g \sec (e+f x))^p}{(a+b \sin (e+f x)) (c+d \sin (e+f x))} \, dx\) [1559]

Optimal result

Integrand size = 35, antiderivative size = 308 \[ \int \frac {(g \sec (e+f x))^p}{(a+b \sin (e+f x)) (c+d \sin (e+f x))} \, dx=-\frac {\operatorname {AppellF1}\left (1+p,\frac {1+p}{2},\frac {1+p}{2},2+p,\frac {a+b}{a+b \sin (e+f x)},\frac {a-b}{a+b \sin (e+f x)}\right ) \sec (e+f x) (g \sec (e+f x))^p \left (-\frac {b (1-\sin (e+f x))}{a+b \sin (e+f x)}\right )^{\frac {1+p}{2}} \left (\frac {b (1+\sin (e+f x))}{a+b \sin (e+f x)}\right )^{\frac {1+p}{2}}}{(b c-a d) f (1+p)}+\frac {\operatorname {AppellF1}\left (1+p,\frac {1+p}{2},\frac {1+p}{2},2+p,\frac {c+d}{c+d \sin (e+f x)},\frac {c-d}{c+d \sin (e+f x)}\right ) \sec (e+f x) (g \sec (e+f x))^p \left (-\frac {d (1-\sin (e+f x))}{c+d \sin (e+f x)}\right )^{\frac {1+p}{2}} \left (\frac {d (1+\sin (e+f x))}{c+d \sin (e+f x)}\right )^{\frac {1+p}{2}}}{(b c-a d) f (1+p)} \] Output:

-AppellF1(p+1,1/2*p+1/2,1/2*p+1/2,2+p,(a-b)/(a+b*sin(f*x+e)),(a+b)/(a+b*si 
n(f*x+e)))*sec(f*x+e)*(g*sec(f*x+e))^p*(-b*(1-sin(f*x+e))/(a+b*sin(f*x+e)) 
)^(1/2*p+1/2)*(b*(1+sin(f*x+e))/(a+b*sin(f*x+e)))^(1/2*p+1/2)/(-a*d+b*c)/f 
/(p+1)+AppellF1(p+1,1/2*p+1/2,1/2*p+1/2,2+p,(c-d)/(c+d*sin(f*x+e)),(c+d)/( 
c+d*sin(f*x+e)))*sec(f*x+e)*(g*sec(f*x+e))^p*(-d*(1-sin(f*x+e))/(c+d*sin(f 
*x+e)))^(1/2*p+1/2)*(d*(1+sin(f*x+e))/(c+d*sin(f*x+e)))^(1/2*p+1/2)/(-a*d+ 
b*c)/f/(p+1)
 

Mathematica [B] (warning: unable to verify)

Leaf count is larger than twice the leaf count of optimal. \(5113\) vs. \(2(308)=616\).

Time = 34.50 (sec) , antiderivative size = 5113, normalized size of antiderivative = 16.60 \[ \int \frac {(g \sec (e+f x))^p}{(a+b \sin (e+f x)) (c+d \sin (e+f x))} \, dx=\text {Result too large to show} \] Input:

Integrate[(g*Sec[e + f*x])^p/((a + b*Sin[e + f*x])*(c + d*Sin[e + f*x])),x 
]
 

Output:

Result too large to show
 

Rubi [A] (verified)

Time = 0.89 (sec) , antiderivative size = 327, normalized size of antiderivative = 1.06, number of steps used = 5, number of rules used = 5, \(\frac {\text {number of rules}}{\text {integrand size}}\) = 0.143, Rules used = {3042, 3405, 3042, 3403, 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 \frac {(g \sec (e+f x))^p}{(a+b \sin (e+f x)) (c+d \sin (e+f x))} \, dx\)

\(\Big \downarrow \) 3042

\(\displaystyle \int \frac {(g \sec (e+f x))^p}{(a+b \sin (e+f x)) (c+d \sin (e+f x))}dx\)

\(\Big \downarrow \) 3405

\(\displaystyle (g \cos (e+f x))^p (g \sec (e+f x))^p \int \frac {(g \cos (e+f x))^{-p}}{(a+b \sin (e+f x)) (c+d \sin (e+f x))}dx\)

\(\Big \downarrow \) 3042

\(\displaystyle (g \cos (e+f x))^p (g \sec (e+f x))^p \int \frac {(g \cos (e+f x))^{-p}}{(a+b \sin (e+f x)) (c+d \sin (e+f x))}dx\)

\(\Big \downarrow \) 3403

\(\displaystyle (g \cos (e+f x))^p (g \sec (e+f x))^p \int \left (\frac {b (g \cos (e+f x))^{-p}}{(b c-a d) (a+b \sin (e+f x))}-\frac {d (g \cos (e+f x))^{-p}}{(b c-a d) (c+d \sin (e+f x))}\right )dx\)

\(\Big \downarrow \) 2009

\(\displaystyle (g \cos (e+f x))^p (g \sec (e+f x))^p \left (\frac {g (g \cos (e+f x))^{-p-1} \left (-\frac {d (1-\sin (e+f x))}{c+d \sin (e+f x)}\right )^{\frac {p+1}{2}} \left (\frac {d (\sin (e+f x)+1)}{c+d \sin (e+f x)}\right )^{\frac {p+1}{2}} \operatorname {AppellF1}\left (p+1,\frac {p+1}{2},\frac {p+1}{2},p+2,\frac {c+d}{c+d \sin (e+f x)},\frac {c-d}{c+d \sin (e+f x)}\right )}{f (p+1) (b c-a d)}-\frac {g (g \cos (e+f x))^{-p-1} \left (-\frac {b (1-\sin (e+f x))}{a+b \sin (e+f x)}\right )^{\frac {p+1}{2}} \left (\frac {b (\sin (e+f x)+1)}{a+b \sin (e+f x)}\right )^{\frac {p+1}{2}} \operatorname {AppellF1}\left (p+1,\frac {p+1}{2},\frac {p+1}{2},p+2,\frac {a+b}{a+b \sin (e+f x)},\frac {a-b}{a+b \sin (e+f x)}\right )}{f (p+1) (b c-a d)}\right )\)

Input:

Int[(g*Sec[e + f*x])^p/((a + b*Sin[e + f*x])*(c + d*Sin[e + f*x])),x]
 

Output:

(g*Cos[e + f*x])^p*(g*Sec[e + f*x])^p*(-((g*AppellF1[1 + p, (1 + p)/2, (1 
+ p)/2, 2 + p, (a + b)/(a + b*Sin[e + f*x]), (a - b)/(a + b*Sin[e + f*x])] 
*(g*Cos[e + f*x])^(-1 - p)*(-((b*(1 - Sin[e + f*x]))/(a + b*Sin[e + f*x])) 
)^((1 + p)/2)*((b*(1 + Sin[e + f*x]))/(a + b*Sin[e + f*x]))^((1 + p)/2))/( 
(b*c - a*d)*f*(1 + p))) + (g*AppellF1[1 + p, (1 + p)/2, (1 + p)/2, 2 + p, 
(c + d)/(c + d*Sin[e + f*x]), (c - d)/(c + d*Sin[e + f*x])]*(g*Cos[e + f*x 
])^(-1 - p)*(-((d*(1 - Sin[e + f*x]))/(c + d*Sin[e + f*x])))^((1 + p)/2)*( 
(d*(1 + Sin[e + f*x]))/(c + d*Sin[e + f*x]))^((1 + p)/2))/((b*c - a*d)*f*( 
1 + p)))
 

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 3403

Int[(cos[(e_.) + (f_.)*(x_)]*(g_.))^(p_)*((a_) + (b_.)*sin[(e_.) + (f_.)*(x 
_)])^(m_)*((c_) + (d_.)*sin[(e_.) + (f_.)*(x_)])^(n_), x_Symbol] :> Int[Exp 
andTrig[(g*cos[e + f*x])^p*(a + b*sin[e + f*x])^m*(c + d*sin[e + f*x])^n, x 
], x] /; FreeQ[{a, b, c, d, e, f, g, p}, x] && NeQ[a^2 - b^2, 0] && Integer 
sQ[2*m, 2*n]
 

rule 3405

Int[((g_.)*sec[(e_.) + (f_.)*(x_)])^(p_)*((a_.) + (b_.)*sin[(e_.) + (f_.)*( 
x_)])^(m_.)*((c_.) + (d_.)*sin[(e_.) + (f_.)*(x_)])^(n_.), x_Symbol] :> Sim 
p[g^(2*IntPart[p])*(g*Cos[e + f*x])^FracPart[p]*(g*Sec[e + f*x])^FracPart[p 
]   Int[(a + b*Sin[e + f*x])^m*((c + d*Sin[e + f*x])^n/(g*Cos[e + f*x])^p), 
 x], x] /; FreeQ[{a, b, c, d, e, f, g, m, n, p}, x] &&  !IntegerQ[p]
 

Maple [F]

\[\int \frac {\left (g \sec \left (f x +e \right )\right )^{p}}{\left (a +b \sin \left (f x +e \right )\right ) \left (c +d \sin \left (f x +e \right )\right )}d x\]

Input:

int((g*sec(f*x+e))^p/(a+b*sin(f*x+e))/(c+d*sin(f*x+e)),x)
 

Output:

int((g*sec(f*x+e))^p/(a+b*sin(f*x+e))/(c+d*sin(f*x+e)),x)
 

Fricas [F]

\[ \int \frac {(g \sec (e+f x))^p}{(a+b \sin (e+f x)) (c+d \sin (e+f x))} \, dx=\int { \frac {\left (g \sec \left (f x + e\right )\right )^{p}}{{\left (b \sin \left (f x + e\right ) + a\right )} {\left (d \sin \left (f x + e\right ) + c\right )}} \,d x } \] Input:

integrate((g*sec(f*x+e))^p/(a+b*sin(f*x+e))/(c+d*sin(f*x+e)),x, algorithm= 
"fricas")
 

Output:

integral(-(g*sec(f*x + e))^p/(b*d*cos(f*x + e)^2 - a*c - b*d - (b*c + a*d) 
*sin(f*x + e)), x)
 

Sympy [F]

\[ \int \frac {(g \sec (e+f x))^p}{(a+b \sin (e+f x)) (c+d \sin (e+f x))} \, dx=\int \frac {\left (g \sec {\left (e + f x \right )}\right )^{p}}{\left (a + b \sin {\left (e + f x \right )}\right ) \left (c + d \sin {\left (e + f x \right )}\right )}\, dx \] Input:

integrate((g*sec(f*x+e))**p/(a+b*sin(f*x+e))/(c+d*sin(f*x+e)),x)
 

Output:

Integral((g*sec(e + f*x))**p/((a + b*sin(e + f*x))*(c + d*sin(e + f*x))), 
x)
 

Maxima [F]

\[ \int \frac {(g \sec (e+f x))^p}{(a+b \sin (e+f x)) (c+d \sin (e+f x))} \, dx=\int { \frac {\left (g \sec \left (f x + e\right )\right )^{p}}{{\left (b \sin \left (f x + e\right ) + a\right )} {\left (d \sin \left (f x + e\right ) + c\right )}} \,d x } \] Input:

integrate((g*sec(f*x+e))^p/(a+b*sin(f*x+e))/(c+d*sin(f*x+e)),x, algorithm= 
"maxima")
 

Output:

integrate((g*sec(f*x + e))^p/((b*sin(f*x + e) + a)*(d*sin(f*x + e) + c)), 
x)
 

Giac [F]

\[ \int \frac {(g \sec (e+f x))^p}{(a+b \sin (e+f x)) (c+d \sin (e+f x))} \, dx=\int { \frac {\left (g \sec \left (f x + e\right )\right )^{p}}{{\left (b \sin \left (f x + e\right ) + a\right )} {\left (d \sin \left (f x + e\right ) + c\right )}} \,d x } \] Input:

integrate((g*sec(f*x+e))^p/(a+b*sin(f*x+e))/(c+d*sin(f*x+e)),x, algorithm= 
"giac")
 

Output:

integrate((g*sec(f*x + e))^p/((b*sin(f*x + e) + a)*(d*sin(f*x + e) + c)), 
x)
 

Mupad [F(-1)]

Timed out. \[ \int \frac {(g \sec (e+f x))^p}{(a+b \sin (e+f x)) (c+d \sin (e+f x))} \, dx=\int \frac {{\left (\frac {g}{\cos \left (e+f\,x\right )}\right )}^p}{\left (a+b\,\sin \left (e+f\,x\right )\right )\,\left (c+d\,\sin \left (e+f\,x\right )\right )} \,d x \] Input:

int((g/cos(e + f*x))^p/((a + b*sin(e + f*x))*(c + d*sin(e + f*x))),x)
 

Output:

int((g/cos(e + f*x))^p/((a + b*sin(e + f*x))*(c + d*sin(e + f*x))), x)
 

Reduce [F]

\[ \int \frac {(g \sec (e+f x))^p}{(a+b \sin (e+f x)) (c+d \sin (e+f x))} \, dx=g^{p} \left (\int \frac {\sec \left (f x +e \right )^{p}}{\sin \left (f x +e \right )^{2} b d +\sin \left (f x +e \right ) a d +\sin \left (f x +e \right ) b c +a c}d x \right ) \] Input:

int((g*sec(f*x+e))^p/(a+b*sin(f*x+e))/(c+d*sin(f*x+e)),x)
 

Output:

g**p*int(sec(e + f*x)**p/(sin(e + f*x)**2*b*d + sin(e + f*x)*a*d + sin(e + 
 f*x)*b*c + a*c),x)