\(\int \frac {A+B x^2+C x^4+D x^6}{x^4 \sqrt {d+e x^2} \sqrt {a+b x^2+c x^4}} \, dx\) [131]

Optimal result

Integrand size = 48, antiderivative size = 719 \[ \int \frac {A+B x^2+C x^4+D x^6}{x^4 \sqrt {d+e x^2} \sqrt {a+b x^2+c x^4}} \, dx=-\frac {A \sqrt {d+e x^2} \sqrt {a+b x^2+c x^4}}{3 a d x^3}-\frac {\sqrt {b^2-4 a c} (3 a B d-2 A (b d+a e)) \sqrt {-\frac {a \left (c+\frac {a}{x^4}+\frac {b}{x^2}\right )}{b^2-4 a c}} x \sqrt {d+e x^2} E\left (\arcsin \left (\frac {\sqrt {1+\frac {b+\frac {2 a}{x^2}}{\sqrt {b^2-4 a c}}}}{\sqrt {2}}\right )|\frac {2 \sqrt {b^2-4 a c} d}{b d+\sqrt {b^2-4 a c} d-2 a e}\right )}{3 \sqrt {2} a^2 d^2 \sqrt {-\frac {a \left (e+\frac {d}{x^2}\right )}{\left (b+\sqrt {b^2-4 a c}\right ) d-2 a e}} \sqrt {a+b x^2+c x^4}}-\frac {\sqrt {2} \sqrt {b^2-4 a c} \left (3 a d (C d-B e)-A \left (c d^2-e (b d+2 a e)\right )\right ) \sqrt {-\frac {a \left (c+\frac {a}{x^4}+\frac {b}{x^2}\right )}{b^2-4 a c}} \sqrt {-\frac {a \left (e+\frac {d}{x^2}\right )}{\left (b+\sqrt {b^2-4 a c}\right ) d-2 a e}} x^3 \operatorname {EllipticF}\left (\arcsin \left (\frac {\sqrt {1+\frac {b+\frac {2 a}{x^2}}{\sqrt {b^2-4 a c}}}}{\sqrt {2}}\right ),\frac {2 \sqrt {b^2-4 a c} d}{b d+\sqrt {b^2-4 a c} d-2 a e}\right )}{3 a^2 d^2 \sqrt {d+e x^2} \sqrt {a+b x^2+c x^4}}+\frac {2 \sqrt {2} \sqrt {b^2-4 a c} D \sqrt {-\frac {a \left (c+\frac {a}{x^4}+\frac {b}{x^2}\right )}{b^2-4 a c}} \sqrt {-\frac {a \left (e+\frac {d}{x^2}\right )}{\left (b+\sqrt {b^2-4 a c}\right ) d-2 a e}} x^3 \operatorname {EllipticPi}\left (\frac {2 \sqrt {b^2-4 a c}}{b+\sqrt {b^2-4 a c}},\arcsin \left (\frac {\sqrt {1+\frac {b+\frac {2 a}{x^2}}{\sqrt {b^2-4 a c}}}}{\sqrt {2}}\right ),\frac {2 \sqrt {b^2-4 a c} d}{b d+\sqrt {b^2-4 a c} d-2 a e}\right )}{\left (b+\sqrt {b^2-4 a c}\right ) \sqrt {d+e x^2} \sqrt {a+b x^2+c x^4}} \] Output:

-1/3*A*(e*x^2+d)^(1/2)*(c*x^4+b*x^2+a)^(1/2)/a/d/x^3-1/6*(-4*a*c+b^2)^(1/2 
)*(3*B*a*d-2*A*(a*e+b*d))*(-a*(c+a/x^4+b/x^2)/(-4*a*c+b^2))^(1/2)*x*(e*x^2 
+d)^(1/2)*EllipticE(1/2*(1+(b+2*a/x^2)/(-4*a*c+b^2)^(1/2))^(1/2)*2^(1/2),2 
^(1/2)*((-4*a*c+b^2)^(1/2)*d/(b*d+(-4*a*c+b^2)^(1/2)*d-2*a*e))^(1/2))*2^(1 
/2)/a^2/d^2/(-a*(e+d/x^2)/((b+(-4*a*c+b^2)^(1/2))*d-2*a*e))^(1/2)/(c*x^4+b 
*x^2+a)^(1/2)-1/3*2^(1/2)*(-4*a*c+b^2)^(1/2)*(3*a*d*(-B*e+C*d)-A*(c*d^2-e* 
(2*a*e+b*d)))*(-a*(c+a/x^4+b/x^2)/(-4*a*c+b^2))^(1/2)*(-a*(e+d/x^2)/((b+(- 
4*a*c+b^2)^(1/2))*d-2*a*e))^(1/2)*x^3*EllipticF(1/2*(1+(b+2*a/x^2)/(-4*a*c 
+b^2)^(1/2))^(1/2)*2^(1/2),2^(1/2)*((-4*a*c+b^2)^(1/2)*d/(b*d+(-4*a*c+b^2) 
^(1/2)*d-2*a*e))^(1/2))/a^2/d^2/(e*x^2+d)^(1/2)/(c*x^4+b*x^2+a)^(1/2)+2*2^ 
(1/2)*(-4*a*c+b^2)^(1/2)*D*(-a*(c+a/x^4+b/x^2)/(-4*a*c+b^2))^(1/2)*(-a*(e+ 
d/x^2)/((b+(-4*a*c+b^2)^(1/2))*d-2*a*e))^(1/2)*x^3*EllipticPi(1/2*(1+(b+2* 
a/x^2)/(-4*a*c+b^2)^(1/2))^(1/2)*2^(1/2),2*(-4*a*c+b^2)^(1/2)/(b+(-4*a*c+b 
^2)^(1/2)),2^(1/2)*((-4*a*c+b^2)^(1/2)*d/(b*d+(-4*a*c+b^2)^(1/2)*d-2*a*e)) 
^(1/2))/(b+(-4*a*c+b^2)^(1/2))/(e*x^2+d)^(1/2)/(c*x^4+b*x^2+a)^(1/2)
 

Mathematica [F]

\[ \int \frac {A+B x^2+C x^4+D x^6}{x^4 \sqrt {d+e x^2} \sqrt {a+b x^2+c x^4}} \, dx=\int \frac {A+B x^2+C x^4+D x^6}{x^4 \sqrt {d+e x^2} \sqrt {a+b x^2+c x^4}} \, dx \] Input:

Integrate[(A + B*x^2 + C*x^4 + D*x^6)/(x^4*Sqrt[d + e*x^2]*Sqrt[a + b*x^2 
+ c*x^4]),x]
 

Output:

Integrate[(A + B*x^2 + C*x^4 + D*x^6)/(x^4*Sqrt[d + e*x^2]*Sqrt[a + b*x^2 
+ c*x^4]), x]
 

Rubi [F]

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 {A+B x^2+C x^4+D x^6}{x^4 \sqrt {d+e x^2} \sqrt {a+b x^2+c x^4}} \, dx\)

\(\Big \downarrow \) 2250

\(\displaystyle \int \frac {A+B x^2+C x^4+D x^6}{x^4 \sqrt {d+e x^2} \sqrt {a+b x^2+c x^4}}dx\)

Input:

Int[(A + B*x^2 + C*x^4 + D*x^6)/(x^4*Sqrt[d + e*x^2]*Sqrt[a + b*x^2 + c*x^ 
4]),x]
 

Output:

$Aborted
 

Defintions of rubi rules used

rule 2250

Int[(Px_)*((f_.)*(x_))^(m_.)*((d_) + (e_.)*(x_)^2)^(q_.)*((a_) + (b_.)*(x_) 
^2 + (c_.)*(x_)^4)^(p_.), x_Symbol] :> Unintegrable[Px*(f*x)^m*(d + e*x^2)^ 
q*(a + b*x^2 + c*x^4)^p, x] /; FreeQ[{a, b, c, d, e, f, m, p, q}, x] && Pol 
yQ[Px, x]
 

Maple [F]

\[\int \frac {D x^{6}+C \,x^{4}+B \,x^{2}+A}{x^{4} \sqrt {e \,x^{2}+d}\, \sqrt {c \,x^{4}+b \,x^{2}+a}}d x\]

Input:

int((D*x^6+C*x^4+B*x^2+A)/x^4/(e*x^2+d)^(1/2)/(c*x^4+b*x^2+a)^(1/2),x)
 

Output:

int((D*x^6+C*x^4+B*x^2+A)/x^4/(e*x^2+d)^(1/2)/(c*x^4+b*x^2+a)^(1/2),x)
 

Fricas [F]

\[ \int \frac {A+B x^2+C x^4+D x^6}{x^4 \sqrt {d+e x^2} \sqrt {a+b x^2+c x^4}} \, dx=\int { \frac {D x^{6} + C x^{4} + B x^{2} + A}{\sqrt {c x^{4} + b x^{2} + a} \sqrt {e x^{2} + d} x^{4}} \,d x } \] Input:

integrate((D*x^6+C*x^4+B*x^2+A)/x^4/(e*x^2+d)^(1/2)/(c*x^4+b*x^2+a)^(1/2), 
x, algorithm="fricas")
 

Output:

integral((D*x^6 + C*x^4 + B*x^2 + A)*sqrt(c*x^4 + b*x^2 + a)*sqrt(e*x^2 + 
d)/(c*e*x^10 + (c*d + b*e)*x^8 + (b*d + a*e)*x^6 + a*d*x^4), x)
 

Sympy [F]

\[ \int \frac {A+B x^2+C x^4+D x^6}{x^4 \sqrt {d+e x^2} \sqrt {a+b x^2+c x^4}} \, dx=\int \frac {A + B x^{2} + C x^{4} + D x^{6}}{x^{4} \sqrt {d + e x^{2}} \sqrt {a + b x^{2} + c x^{4}}}\, dx \] Input:

integrate((D*x**6+C*x**4+B*x**2+A)/x**4/(e*x**2+d)**(1/2)/(c*x**4+b*x**2+a 
)**(1/2),x)
 

Output:

Integral((A + B*x**2 + C*x**4 + D*x**6)/(x**4*sqrt(d + e*x**2)*sqrt(a + b* 
x**2 + c*x**4)), x)
 

Maxima [F]

\[ \int \frac {A+B x^2+C x^4+D x^6}{x^4 \sqrt {d+e x^2} \sqrt {a+b x^2+c x^4}} \, dx=\int { \frac {D x^{6} + C x^{4} + B x^{2} + A}{\sqrt {c x^{4} + b x^{2} + a} \sqrt {e x^{2} + d} x^{4}} \,d x } \] Input:

integrate((D*x^6+C*x^4+B*x^2+A)/x^4/(e*x^2+d)^(1/2)/(c*x^4+b*x^2+a)^(1/2), 
x, algorithm="maxima")
 

Output:

integrate((D*x^6 + C*x^4 + B*x^2 + A)/(sqrt(c*x^4 + b*x^2 + a)*sqrt(e*x^2 
+ d)*x^4), x)
 

Giac [F]

\[ \int \frac {A+B x^2+C x^4+D x^6}{x^4 \sqrt {d+e x^2} \sqrt {a+b x^2+c x^4}} \, dx=\int { \frac {D x^{6} + C x^{4} + B x^{2} + A}{\sqrt {c x^{4} + b x^{2} + a} \sqrt {e x^{2} + d} x^{4}} \,d x } \] Input:

integrate((D*x^6+C*x^4+B*x^2+A)/x^4/(e*x^2+d)^(1/2)/(c*x^4+b*x^2+a)^(1/2), 
x, algorithm="giac")
 

Output:

integrate((D*x^6 + C*x^4 + B*x^2 + A)/(sqrt(c*x^4 + b*x^2 + a)*sqrt(e*x^2 
+ d)*x^4), x)
 

Mupad [F(-1)]

Timed out. \[ \int \frac {A+B x^2+C x^4+D x^6}{x^4 \sqrt {d+e x^2} \sqrt {a+b x^2+c x^4}} \, dx=\int \frac {A+B\,x^2+C\,x^4+x^6\,D}{x^4\,\sqrt {e\,x^2+d}\,\sqrt {c\,x^4+b\,x^2+a}} \,d x \] Input:

int((A + B*x^2 + C*x^4 + x^6*D)/(x^4*(d + e*x^2)^(1/2)*(a + b*x^2 + c*x^4) 
^(1/2)),x)
 

Output:

int((A + B*x^2 + C*x^4 + x^6*D)/(x^4*(d + e*x^2)^(1/2)*(a + b*x^2 + c*x^4) 
^(1/2)), x)
 

Reduce [F]

\[ \int \frac {A+B x^2+C x^4+D x^6}{x^4 \sqrt {d+e x^2} \sqrt {a+b x^2+c x^4}} \, dx =\text {Too large to display} \] Input:

int((D*x^6+C*x^4+B*x^2+A)/x^4/(e*x^2+d)^(1/2)/(c*x^4+b*x^2+a)^(1/2),x)
 

Output:

( - sqrt(d + e*x**2)*sqrt(a + b*x**2 + c*x**4)*b + 2*int((sqrt(d + e*x**2) 
*sqrt(a + b*x**2 + c*x**4)*x**2)/(a**2*d*e + a**2*e**2*x**2 + a*b*d**2 + 2 
*a*b*d*e*x**2 + a*b*e**2*x**4 + a*c*d*e*x**4 + a*c*e**2*x**6 + b**2*d**2*x 
**2 + b**2*d*e*x**4 + b*c*d**2*x**4 + b*c*d*e*x**6),x)*a**2*d*e**2*x**3 + 
4*int((sqrt(d + e*x**2)*sqrt(a + b*x**2 + c*x**4)*x**2)/(a**2*d*e + a**2*e 
**2*x**2 + a*b*d**2 + 2*a*b*d*e*x**2 + a*b*e**2*x**4 + a*c*d*e*x**4 + a*c* 
e**2*x**6 + b**2*d**2*x**2 + b**2*d*e*x**4 + b*c*d**2*x**4 + b*c*d*e*x**6) 
,x)*a*b*d**2*e*x**3 + 2*int((sqrt(d + e*x**2)*sqrt(a + b*x**2 + c*x**4)*x* 
*2)/(a**2*d*e + a**2*e**2*x**2 + a*b*d**2 + 2*a*b*d*e*x**2 + a*b*e**2*x**4 
 + a*c*d*e*x**4 + a*c*e**2*x**6 + b**2*d**2*x**2 + b**2*d*e*x**4 + b*c*d** 
2*x**4 + b*c*d*e*x**6),x)*b**2*d**3*x**3 + 2*int((sqrt(d + e*x**2)*sqrt(a 
+ b*x**2 + c*x**4))/(a**2*d*e*x**4 + a**2*e**2*x**6 + a*b*d**2*x**4 + 2*a* 
b*d*e*x**6 + a*b*e**2*x**8 + a*c*d*e*x**8 + a*c*e**2*x**10 + b**2*d**2*x** 
6 + b**2*d*e*x**8 + b*c*d**2*x**8 + b*c*d*e*x**10),x)*a**3*e**2*x**3 + int 
((sqrt(d + e*x**2)*sqrt(a + b*x**2 + c*x**4))/(a**2*d*e*x**4 + a**2*e**2*x 
**6 + a*b*d**2*x**4 + 2*a*b*d*e*x**6 + a*b*e**2*x**8 + a*c*d*e*x**8 + a*c* 
e**2*x**10 + b**2*d**2*x**6 + b**2*d*e*x**8 + b*c*d**2*x**8 + b*c*d*e*x**1 
0),x)*a**2*b*d*e*x**3 - int((sqrt(d + e*x**2)*sqrt(a + b*x**2 + c*x**4))/( 
a**2*d*e*x**4 + a**2*e**2*x**6 + a*b*d**2*x**4 + 2*a*b*d*e*x**6 + a*b*e**2 
*x**8 + a*c*d*e*x**8 + a*c*e**2*x**10 + b**2*d**2*x**6 + b**2*d*e*x**8 ...