\(\int \frac {1}{(a+b \sin (e+f x))^{3/2} \sqrt {c+d \sin (e+f x)}} \, dx\) [799]

Optimal result
Mathematica [B] (warning: unable to verify)
Rubi [A] (verified)
Maple [B] (warning: unable to verify)
Fricas [F]
Sympy [F]
Maxima [F]
Giac [F]
Mupad [F(-1)]
Reduce [F]

Optimal result

Integrand size = 29, antiderivative size = 405 \[ \int \frac {1}{(a+b \sin (e+f x))^{3/2} \sqrt {c+d \sin (e+f x)}} \, dx=\frac {2 b (c-d) \sqrt {c+d} E\left (\arcsin \left (\frac {\sqrt {a+b} \sqrt {c+d \sin (e+f x)}}{\sqrt {c+d} \sqrt {a+b \sin (e+f x)}}\right )|\frac {(a-b) (c+d)}{(a+b) (c-d)}\right ) \sec (e+f x) \sqrt {-\frac {(b c-a d) (1-\sin (e+f x))}{(c+d) (a+b \sin (e+f x))}} \sqrt {\frac {(b c-a d) (1+\sin (e+f x))}{(c-d) (a+b \sin (e+f x))}} (a+b \sin (e+f x))}{(a-b) \sqrt {a+b} (b c-a d)^2 f}+\frac {2 \sqrt {a+b} \operatorname {EllipticF}\left (\arcsin \left (\frac {\sqrt {c+d} \sqrt {a+b \sin (e+f x)}}{\sqrt {a+b} \sqrt {c+d \sin (e+f x)}}\right ),\frac {(a+b) (c-d)}{(a-b) (c+d)}\right ) \sec (e+f x) \sqrt {\frac {(b c-a d) (1-\sin (e+f x))}{(a+b) (c+d \sin (e+f x))}} \sqrt {-\frac {(b c-a d) (1+\sin (e+f x))}{(a-b) (c+d \sin (e+f x))}} (c+d \sin (e+f x))}{(a-b) \sqrt {c+d} (b c-a d) f} \] Output:

2*b*(c-d)*(c+d)^(1/2)*EllipticE((a+b)^(1/2)*(c+d*sin(f*x+e))^(1/2)/(c+d)^( 
1/2)/(a+b*sin(f*x+e))^(1/2),((a-b)*(c+d)/(a+b)/(c-d))^(1/2))*sec(f*x+e)*(- 
(-a*d+b*c)*(1-sin(f*x+e))/(c+d)/(a+b*sin(f*x+e)))^(1/2)*((-a*d+b*c)*(1+sin 
(f*x+e))/(c-d)/(a+b*sin(f*x+e)))^(1/2)*(a+b*sin(f*x+e))/(a-b)/(a+b)^(1/2)/ 
(-a*d+b*c)^2/f+2*(a+b)^(1/2)*EllipticF((c+d)^(1/2)*(a+b*sin(f*x+e))^(1/2)/ 
(a+b)^(1/2)/(c+d*sin(f*x+e))^(1/2),((a+b)*(c-d)/(a-b)/(c+d))^(1/2))*sec(f* 
x+e)*((-a*d+b*c)*(1-sin(f*x+e))/(a+b)/(c+d*sin(f*x+e)))^(1/2)*(-(-a*d+b*c) 
*(1+sin(f*x+e))/(a-b)/(c+d*sin(f*x+e)))^(1/2)*(c+d*sin(f*x+e))/(a-b)/(c+d) 
^(1/2)/(-a*d+b*c)/f
 

Mathematica [B] (warning: unable to verify)

Leaf count is larger than twice the leaf count of optimal. \(88729\) vs. \(2(405)=810\).

Time = 39.08 (sec) , antiderivative size = 88729, normalized size of antiderivative = 219.08 \[ \int \frac {1}{(a+b \sin (e+f x))^{3/2} \sqrt {c+d \sin (e+f x)}} \, dx=\text {Result too large to show} \] Input:

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

Output:

Result too large to show
 

Rubi [A] (verified)

Time = 0.96 (sec) , antiderivative size = 405, normalized size of antiderivative = 1.00, number of steps used = 5, number of rules used = 5, \(\frac {\text {number of rules}}{\text {integrand size}}\) = 0.172, Rules used = {3042, 3280, 3042, 3297, 3475}

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

\(\Big \downarrow \) 3042

\(\displaystyle \int \frac {1}{(a+b \sin (e+f x))^{3/2} \sqrt {c+d \sin (e+f x)}}dx\)

\(\Big \downarrow \) 3280

\(\displaystyle \frac {\int \frac {1}{\sqrt {a+b \sin (e+f x)} \sqrt {c+d \sin (e+f x)}}dx}{a-b}-\frac {b \int \frac {\sin (e+f x)+1}{(a+b \sin (e+f x))^{3/2} \sqrt {c+d \sin (e+f x)}}dx}{a-b}\)

\(\Big \downarrow \) 3042

\(\displaystyle \frac {\int \frac {1}{\sqrt {a+b \sin (e+f x)} \sqrt {c+d \sin (e+f x)}}dx}{a-b}-\frac {b \int \frac {\sin (e+f x)+1}{(a+b \sin (e+f x))^{3/2} \sqrt {c+d \sin (e+f x)}}dx}{a-b}\)

\(\Big \downarrow \) 3297

\(\displaystyle \frac {2 \sqrt {a+b} \sec (e+f x) (c+d \sin (e+f x)) \sqrt {\frac {(b c-a d) (1-\sin (e+f x))}{(a+b) (c+d \sin (e+f x))}} \sqrt {-\frac {(b c-a d) (\sin (e+f x)+1)}{(a-b) (c+d \sin (e+f x))}} \operatorname {EllipticF}\left (\arcsin \left (\frac {\sqrt {c+d} \sqrt {a+b \sin (e+f x)}}{\sqrt {a+b} \sqrt {c+d \sin (e+f x)}}\right ),\frac {(a+b) (c-d)}{(a-b) (c+d)}\right )}{f (a-b) \sqrt {c+d} (b c-a d)}-\frac {b \int \frac {\sin (e+f x)+1}{(a+b \sin (e+f x))^{3/2} \sqrt {c+d \sin (e+f x)}}dx}{a-b}\)

\(\Big \downarrow \) 3475

\(\displaystyle \frac {2 \sqrt {a+b} \sec (e+f x) (c+d \sin (e+f x)) \sqrt {\frac {(b c-a d) (1-\sin (e+f x))}{(a+b) (c+d \sin (e+f x))}} \sqrt {-\frac {(b c-a d) (\sin (e+f x)+1)}{(a-b) (c+d \sin (e+f x))}} \operatorname {EllipticF}\left (\arcsin \left (\frac {\sqrt {c+d} \sqrt {a+b \sin (e+f x)}}{\sqrt {a+b} \sqrt {c+d \sin (e+f x)}}\right ),\frac {(a+b) (c-d)}{(a-b) (c+d)}\right )}{f (a-b) \sqrt {c+d} (b c-a d)}+\frac {2 b (c-d) \sqrt {c+d} \sec (e+f x) (a+b \sin (e+f x)) \sqrt {-\frac {(b c-a d) (1-\sin (e+f x))}{(c+d) (a+b \sin (e+f x))}} \sqrt {\frac {(b c-a d) (\sin (e+f x)+1)}{(c-d) (a+b \sin (e+f x))}} E\left (\arcsin \left (\frac {\sqrt {a+b} \sqrt {c+d \sin (e+f x)}}{\sqrt {c+d} \sqrt {a+b \sin (e+f x)}}\right )|\frac {(a-b) (c+d)}{(a+b) (c-d)}\right )}{f (a-b) \sqrt {a+b} (b c-a d)^2}\)

Input:

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

Output:

(2*b*(c - d)*Sqrt[c + d]*EllipticE[ArcSin[(Sqrt[a + b]*Sqrt[c + d*Sin[e + 
f*x]])/(Sqrt[c + d]*Sqrt[a + b*Sin[e + f*x]])], ((a - b)*(c + d))/((a + b) 
*(c - d))]*Sec[e + f*x]*Sqrt[-(((b*c - a*d)*(1 - Sin[e + f*x]))/((c + d)*( 
a + b*Sin[e + f*x])))]*Sqrt[((b*c - a*d)*(1 + Sin[e + f*x]))/((c - d)*(a + 
 b*Sin[e + f*x]))]*(a + b*Sin[e + f*x]))/((a - b)*Sqrt[a + b]*(b*c - a*d)^ 
2*f) + (2*Sqrt[a + b]*EllipticF[ArcSin[(Sqrt[c + d]*Sqrt[a + b*Sin[e + f*x 
]])/(Sqrt[a + b]*Sqrt[c + d*Sin[e + f*x]])], ((a + b)*(c - d))/((a - b)*(c 
 + d))]*Sec[e + f*x]*Sqrt[((b*c - a*d)*(1 - Sin[e + f*x]))/((a + b)*(c + d 
*Sin[e + f*x]))]*Sqrt[-(((b*c - a*d)*(1 + Sin[e + f*x]))/((a - b)*(c + d*S 
in[e + f*x])))]*(c + d*Sin[e + f*x]))/((a - b)*Sqrt[c + d]*(b*c - a*d)*f)
 

Defintions of rubi rules used

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

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

rule 3297
Int[1/(Sqrt[(a_) + (b_.)*sin[(e_.) + (f_.)*(x_)]]*Sqrt[(c_) + (d_.)*sin[(e_ 
.) + (f_.)*(x_)]]), x_Symbol] :> Simp[2*((c + d*Sin[e + f*x])/(f*(b*c - a*d 
)*Rt[(c + d)/(a + b), 2]*Cos[e + f*x]))*Sqrt[(b*c - a*d)*((1 - Sin[e + f*x] 
)/((a + b)*(c + d*Sin[e + f*x])))]*Sqrt[(-(b*c - a*d))*((1 + Sin[e + f*x])/ 
((a - b)*(c + d*Sin[e + f*x])))]*EllipticF[ArcSin[Rt[(c + d)/(a + b), 2]*(S 
qrt[a + b*Sin[e + f*x]]/Sqrt[c + d*Sin[e + f*x]])], (a + b)*((c - d)/((a - 
b)*(c + d)))], x] /; FreeQ[{a, b, c, d, e, f}, x] && NeQ[b*c - a*d, 0] && N 
eQ[a^2 - b^2, 0] && NeQ[c^2 - d^2, 0] && PosQ[(c + d)/(a + b)]
 

rule 3475
Int[((A_) + (B_.)*sin[(e_.) + (f_.)*(x_)])/(((a_) + (b_.)*sin[(e_.) + (f_.) 
*(x_)])^(3/2)*Sqrt[(c_) + (d_.)*sin[(e_.) + (f_.)*(x_)]]), x_Symbol] :> Sim 
p[-2*A*(c - d)*((a + b*Sin[e + f*x])/(f*(b*c - a*d)^2*Rt[(a + b)/(c + d), 2 
]*Cos[e + f*x]))*Sqrt[(b*c - a*d)*((1 + Sin[e + f*x])/((c - d)*(a + b*Sin[e 
 + f*x])))]*Sqrt[(-(b*c - a*d))*((1 - Sin[e + f*x])/((c + d)*(a + b*Sin[e + 
 f*x])))]*EllipticE[ArcSin[Rt[(a + b)/(c + d), 2]*(Sqrt[c + d*Sin[e + f*x]] 
/Sqrt[a + b*Sin[e + f*x]])], (a - b)*((c + d)/((a + b)*(c - d)))], x] /; Fr 
eeQ[{a, b, c, d, e, f, A, B}, x] && NeQ[b*c - a*d, 0] && NeQ[a^2 - b^2, 0] 
&& NeQ[c^2 - d^2, 0] && EqQ[A, B] && PosQ[(a + b)/(c + d)]
 
Maple [B] (warning: unable to verify)

Leaf count of result is larger than twice the leaf count of optimal. \(36511\) vs. \(2(375)=750\).

Time = 11.82 (sec) , antiderivative size = 36512, normalized size of antiderivative = 90.15

method result size
default \(\text {Expression too large to display}\) \(36512\)

Input:

int(1/(a+b*sin(f*x+e))^(3/2)/(c+d*sin(f*x+e))^(1/2),x,method=_RETURNVERBOS 
E)
 

Output:

result too large to display
 

Fricas [F]

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

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

Output:

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

Sympy [F]

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

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

Output:

Integral(1/((a + b*sin(e + f*x))**(3/2)*sqrt(c + d*sin(e + f*x))), x)
 

Maxima [F]

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

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

Output:

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

Giac [F]

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

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

Output:

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

Mupad [F(-1)]

Timed out. \[ \int \frac {1}{(a+b \sin (e+f x))^{3/2} \sqrt {c+d \sin (e+f x)}} \, dx=\int \frac {1}{{\left (a+b\,\sin \left (e+f\,x\right )\right )}^{3/2}\,\sqrt {c+d\,\sin \left (e+f\,x\right )}} \,d x \] Input:

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

Output:

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

Reduce [F]

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

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

Output:

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