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

Optimal result

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

c^(1/2)*(-a*d+b*c)^(1/2)*(a*(d*x^2+c)/c/(b*x^2+a))^(1/2)*(f*x^2+e)^(1/2)*E 
llipticE((-a*d+b*c)^(1/2)*x/c^(1/2)/(b*x^2+a)^(1/2),(c*(-a*f+b*e)/(-a*d+b* 
c)/e)^(1/2))/a/b/(d*x^2+c)^(1/2)/(a*(f*x^2+e)/e/(b*x^2+a))^(1/2)+c^(1/2)*d 
*(-a*f+b*e)*(a*(d*x^2+c)/c/(b*x^2+a))^(1/2)*(f*x^2+e)^(1/2)*EllipticF((-a* 
d+b*c)^(1/2)*x/c^(1/2)/(b*x^2+a)^(1/2),(c*(-a*f+b*e)/(-a*d+b*c)/e)^(1/2))/ 
b^2/(-a*d+b*c)^(1/2)/e/(d*x^2+c)^(1/2)/(a*(f*x^2+e)/e/(b*x^2+a))^(1/2)+a*c 
^(1/2)*d*f*(a*(d*x^2+c)/c/(b*x^2+a))^(1/2)*(f*x^2+e)^(1/2)*EllipticPi((-a* 
d+b*c)^(1/2)*x/c^(1/2)/(b*x^2+a)^(1/2),b*c/(-a*d+b*c),(c*(-a*f+b*e)/(-a*d+ 
b*c)/e)^(1/2))/b^2/(-a*d+b*c)^(1/2)/e/(d*x^2+c)^(1/2)/(a*(f*x^2+e)/e/(b*x^ 
2+a))^(1/2)
 

Mathematica [F]

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

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

Output:

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

Rubi [A] (verified)

Time = 0.48 (sec) , antiderivative size = 316, normalized size of antiderivative = 0.66, number of steps used = 6, number of rules used = 5, \(\frac {\text {number of rules}}{\text {integrand size}}\) = 0.147, Rules used = {432, 428, 412, 429, 327}

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

\(\Big \downarrow \) 432

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

\(\Big \downarrow \) 428

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

\(\Big \downarrow \) 412

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

\(\Big \downarrow \) 429

\(\displaystyle \frac {\sqrt {e+f x^2} (b c-a d) \sqrt {\frac {a \left (c+d x^2\right )}{c \left (a+b x^2\right )}} \int \frac {\sqrt {1-\frac {(b e-a f) x^2}{e \left (b x^2+a\right )}}}{\sqrt {1-\frac {(b c-a d) x^2}{c \left (b x^2+a\right )}}}d\frac {x}{\sqrt {b x^2+a}}}{a b \sqrt {c+d x^2} \sqrt {\frac {a \left (e+f x^2\right )}{e \left (a+b x^2\right )}}}+\frac {d e \sqrt {a+b x^2} \sqrt {\frac {e \left (c+d x^2\right )}{c \left (e+f x^2\right )}} \operatorname {EllipticPi}\left (-\frac {a f}{b e-a f},\arcsin \left (\frac {\sqrt {a f-b e} x}{\sqrt {a} \sqrt {f x^2+e}}\right ),\frac {a (d e-c f)}{c (b e-a f)}\right )}{\sqrt {a} b \sqrt {c+d x^2} \sqrt {a f-b e} \sqrt {\frac {e \left (a+b x^2\right )}{a \left (e+f x^2\right )}}}\)

\(\Big \downarrow \) 327

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

Input:

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

Output:

(Sqrt[c]*Sqrt[b*c - a*d]*Sqrt[(a*(c + d*x^2))/(c*(a + b*x^2))]*Sqrt[e + f* 
x^2]*EllipticE[ArcSin[(Sqrt[b*c - a*d]*x)/(Sqrt[c]*Sqrt[a + b*x^2])], (c*( 
b*e - a*f))/((b*c - a*d)*e)])/(a*b*Sqrt[c + d*x^2]*Sqrt[(a*(e + f*x^2))/(e 
*(a + b*x^2))]) + (d*e*Sqrt[a + b*x^2]*Sqrt[(e*(c + d*x^2))/(c*(e + f*x^2) 
)]*EllipticPi[-((a*f)/(b*e - a*f)), ArcSin[(Sqrt[-(b*e) + a*f]*x)/(Sqrt[a] 
*Sqrt[e + f*x^2])], (a*(d*e - c*f))/(c*(b*e - a*f))])/(Sqrt[a]*b*Sqrt[-(b* 
e) + a*f]*Sqrt[c + d*x^2]*Sqrt[(e*(a + b*x^2))/(a*(e + f*x^2))])
 

Defintions of rubi rules used

rule 327

Int[Sqrt[(a_) + (b_.)*(x_)^2]/Sqrt[(c_) + (d_.)*(x_)^2], x_Symbol] :> Simp[ 
(Sqrt[a]/(Sqrt[c]*Rt[-d/c, 2]))*EllipticE[ArcSin[Rt[-d/c, 2]*x], b*(c/(a*d) 
)], x] /; FreeQ[{a, b, c, d}, x] && NegQ[d/c] && GtQ[c, 0] && GtQ[a, 0]
 

rule 412

Int[1/(((a_) + (b_.)*(x_)^2)*Sqrt[(c_) + (d_.)*(x_)^2]*Sqrt[(e_) + (f_.)*(x 
_)^2]), x_Symbol] :> Simp[(1/(a*Sqrt[c]*Sqrt[e]*Rt[-d/c, 2]))*EllipticPi[b* 
(c/(a*d)), ArcSin[Rt[-d/c, 2]*x], c*(f/(d*e))], x] /; FreeQ[{a, b, c, d, e, 
 f}, x] &&  !GtQ[d/c, 0] && GtQ[c, 0] && GtQ[e, 0] &&  !( !GtQ[f/e, 0] && S 
implerSqrtQ[-f/e, -d/c])
 

rule 428

Int[Sqrt[(a_) + (b_.)*(x_)^2]/(Sqrt[(c_) + (d_.)*(x_)^2]*Sqrt[(e_) + (f_.)* 
(x_)^2]), x_Symbol] :> Simp[a*Sqrt[c + d*x^2]*(Sqrt[a*((e + f*x^2)/(e*(a + 
b*x^2)))]/(c*Sqrt[e + f*x^2]*Sqrt[a*((c + d*x^2)/(c*(a + b*x^2)))]))   Subs 
t[Int[1/((1 - b*x^2)*Sqrt[1 - (b*c - a*d)*(x^2/c)]*Sqrt[1 - (b*e - a*f)*(x^ 
2/e)]), x], x, x/Sqrt[a + b*x^2]], x] /; FreeQ[{a, b, c, d, e, f}, x]
 

rule 429

Int[Sqrt[(c_) + (d_.)*(x_)^2]/(((a_) + (b_.)*(x_)^2)^(3/2)*Sqrt[(e_) + (f_. 
)*(x_)^2]), x_Symbol] :> Simp[Sqrt[c + d*x^2]*(Sqrt[a*((e + f*x^2)/(e*(a + 
b*x^2)))]/(a*Sqrt[e + f*x^2]*Sqrt[a*((c + d*x^2)/(c*(a + b*x^2)))]))   Subs 
t[Int[Sqrt[1 - (b*c - a*d)*(x^2/c)]/Sqrt[1 - (b*e - a*f)*(x^2/e)], x], x, x 
/Sqrt[a + b*x^2]], x] /; FreeQ[{a, b, c, d, e, f}, x]
 

rule 432

Int[(Sqrt[(a_) + (b_.)*(x_)^2]*Sqrt[(c_) + (d_.)*(x_)^2])/((e_) + (f_.)*(x_ 
)^2)^(3/2), x_Symbol] :> Simp[b/f   Int[Sqrt[c + d*x^2]/(Sqrt[a + b*x^2]*Sq 
rt[e + f*x^2]), x], x] - Simp[(b*e - a*f)/f   Int[Sqrt[c + d*x^2]/(Sqrt[a + 
 b*x^2]*(e + f*x^2)^(3/2)), x], x] /; FreeQ[{a, b, c, d, e, f}, x]
 

Maple [F]

\[\int \frac {\sqrt {x^{2} d +c}\, \sqrt {f \,x^{2}+e}}{\left (b \,x^{2}+a \right )^{\frac {3}{2}}}d x\]

Input:

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

Output:

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

Fricas [F(-1)]

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

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

Output:

Timed out
 

Sympy [F]

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

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

Output:

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

Maxima [F]

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

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

Output:

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

Giac [F]

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

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

Output:

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

Mupad [F(-1)]

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

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

Output:

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

Reduce [F]

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

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

Output:

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