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

Optimal result

Integrand size = 29, antiderivative size = 221 \[ \int \frac {A+B x^2+C x^4+D x^6}{\sqrt [4]{a+b x^2}} \, dx=\frac {2 \left (195 A b^3-78 a b^2 B+52 a^2 b C-40 a^3 D\right ) x}{195 b^3 \sqrt [4]{a+b x^2}}+\frac {2 \left (39 b^2 B-26 a b C+20 a^2 D\right ) x \left (a+b x^2\right )^{3/4}}{195 b^3}+\frac {2 (13 b C-10 a D) x^3 \left (a+b x^2\right )^{3/4}}{117 b^2}+\frac {2 D x^5 \left (a+b x^2\right )^{3/4}}{13 b}-\frac {2 \sqrt {a} \left (195 A b^3-78 a b^2 B+52 a^2 b C-40 a^3 D\right ) \sqrt [4]{1+\frac {b x^2}{a}} E\left (\left .\frac {1}{2} \arctan \left (\frac {\sqrt {b} x}{\sqrt {a}}\right )\right |2\right )}{195 b^{7/2} \sqrt [4]{a+b x^2}} \] Output:

2/195*(195*A*b^3-78*B*a*b^2+52*C*a^2*b-40*D*a^3)*x/b^3/(b*x^2+a)^(1/4)+2/1 
95*(39*B*b^2-26*C*a*b+20*D*a^2)*x*(b*x^2+a)^(3/4)/b^3+2/117*(13*C*b-10*D*a 
)*x^3*(b*x^2+a)^(3/4)/b^2+2/13*D*x^5*(b*x^2+a)^(3/4)/b-2/195*a^(1/2)*(195* 
A*b^3-78*B*a*b^2+52*C*a^2*b-40*D*a^3)*(1+b*x^2/a)^(1/4)*EllipticE(sin(1/2* 
arctan(b^(1/2)*x/a^(1/2))),2^(1/2))/b^(7/2)/(b*x^2+a)^(1/4)
                                                                                    
                                                                                    
 

Mathematica [C] (verified)

Result contains higher order function than in optimal. Order 5 vs. order 4 in optimal.

Time = 10.11 (sec) , antiderivative size = 132, normalized size of antiderivative = 0.60 \[ \int \frac {A+B x^2+C x^4+D x^6}{\sqrt [4]{a+b x^2}} \, dx=\frac {x \left (2 \left (a+b x^2\right ) \left (60 a^2 D-2 a b \left (39 C+25 D x^2\right )+b^2 \left (117 B+65 C x^2+45 D x^4\right )\right )+3 \left (195 A b^3-78 a b^2 B+52 a^2 b C-40 a^3 D\right ) \sqrt [4]{1+\frac {b x^2}{a}} \operatorname {Hypergeometric2F1}\left (\frac {1}{4},\frac {1}{2},\frac {3}{2},-\frac {b x^2}{a}\right )\right )}{585 b^3 \sqrt [4]{a+b x^2}} \] Input:

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

Output:

(x*(2*(a + b*x^2)*(60*a^2*D - 2*a*b*(39*C + 25*D*x^2) + b^2*(117*B + 65*C* 
x^2 + 45*D*x^4)) + 3*(195*A*b^3 - 78*a*b^2*B + 52*a^2*b*C - 40*a^3*D)*(1 + 
 (b*x^2)/a)^(1/4)*Hypergeometric2F1[1/4, 1/2, 3/2, -((b*x^2)/a)]))/(585*b^ 
3*(a + b*x^2)^(1/4))
 

Rubi [A] (verified)

Time = 0.36 (sec) , antiderivative size = 214, normalized size of antiderivative = 0.97, number of steps used = 8, number of rules used = 8, \(\frac {\text {number of rules}}{\text {integrand size}}\) = 0.276, Rules used = {2346, 27, 1473, 27, 299, 227, 225, 212}

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}{\sqrt [4]{a+b x^2}} \, dx\)

\(\Big \downarrow \) 2346

\(\displaystyle \frac {2 \int \frac {(13 b C-10 a D) x^4+13 b B x^2+13 A b}{2 \sqrt [4]{b x^2+a}}dx}{13 b}+\frac {2 D x^5 \left (a+b x^2\right )^{3/4}}{13 b}\)

\(\Big \downarrow \) 27

\(\displaystyle \frac {\int \frac {(13 b C-10 a D) x^4+13 b B x^2+13 A b}{\sqrt [4]{b x^2+a}}dx}{13 b}+\frac {2 D x^5 \left (a+b x^2\right )^{3/4}}{13 b}\)

\(\Big \downarrow \) 1473

\(\displaystyle \frac {\frac {2 \int \frac {3 \left (39 A b^2+\left (20 D a^2-26 b C a+39 b^2 B\right ) x^2\right )}{2 \sqrt [4]{b x^2+a}}dx}{9 b}+\frac {2 x^3 \left (a+b x^2\right )^{3/4} (13 b C-10 a D)}{9 b}}{13 b}+\frac {2 D x^5 \left (a+b x^2\right )^{3/4}}{13 b}\)

\(\Big \downarrow \) 27

\(\displaystyle \frac {\frac {\int \frac {39 A b^2+\left (20 D a^2-26 b C a+39 b^2 B\right ) x^2}{\sqrt [4]{b x^2+a}}dx}{3 b}+\frac {2 x^3 \left (a+b x^2\right )^{3/4} (13 b C-10 a D)}{9 b}}{13 b}+\frac {2 D x^5 \left (a+b x^2\right )^{3/4}}{13 b}\)

\(\Big \downarrow \) 299

\(\displaystyle \frac {\frac {\frac {\left (-40 a^3 D+52 a^2 b C-78 a b^2 B+195 A b^3\right ) \int \frac {1}{\sqrt [4]{b x^2+a}}dx}{5 b}+\frac {2 x \left (a+b x^2\right )^{3/4} \left (20 a^2 D-26 a b C+39 b^2 B\right )}{5 b}}{3 b}+\frac {2 x^3 \left (a+b x^2\right )^{3/4} (13 b C-10 a D)}{9 b}}{13 b}+\frac {2 D x^5 \left (a+b x^2\right )^{3/4}}{13 b}\)

\(\Big \downarrow \) 227

\(\displaystyle \frac {\frac {\frac {\sqrt [4]{\frac {b x^2}{a}+1} \left (-40 a^3 D+52 a^2 b C-78 a b^2 B+195 A b^3\right ) \int \frac {1}{\sqrt [4]{\frac {b x^2}{a}+1}}dx}{5 b \sqrt [4]{a+b x^2}}+\frac {2 x \left (a+b x^2\right )^{3/4} \left (20 a^2 D-26 a b C+39 b^2 B\right )}{5 b}}{3 b}+\frac {2 x^3 \left (a+b x^2\right )^{3/4} (13 b C-10 a D)}{9 b}}{13 b}+\frac {2 D x^5 \left (a+b x^2\right )^{3/4}}{13 b}\)

\(\Big \downarrow \) 225

\(\displaystyle \frac {\frac {\frac {\sqrt [4]{\frac {b x^2}{a}+1} \left (-40 a^3 D+52 a^2 b C-78 a b^2 B+195 A b^3\right ) \left (\frac {2 x}{\sqrt [4]{\frac {b x^2}{a}+1}}-\int \frac {1}{\left (\frac {b x^2}{a}+1\right )^{5/4}}dx\right )}{5 b \sqrt [4]{a+b x^2}}+\frac {2 x \left (a+b x^2\right )^{3/4} \left (20 a^2 D-26 a b C+39 b^2 B\right )}{5 b}}{3 b}+\frac {2 x^3 \left (a+b x^2\right )^{3/4} (13 b C-10 a D)}{9 b}}{13 b}+\frac {2 D x^5 \left (a+b x^2\right )^{3/4}}{13 b}\)

\(\Big \downarrow \) 212

\(\displaystyle \frac {\frac {\frac {2 x \left (a+b x^2\right )^{3/4} \left (20 a^2 D-26 a b C+39 b^2 B\right )}{5 b}+\frac {\sqrt [4]{\frac {b x^2}{a}+1} \left (\frac {2 x}{\sqrt [4]{\frac {b x^2}{a}+1}}-\frac {2 \sqrt {a} E\left (\left .\frac {1}{2} \arctan \left (\frac {\sqrt {b} x}{\sqrt {a}}\right )\right |2\right )}{\sqrt {b}}\right ) \left (-40 a^3 D+52 a^2 b C-78 a b^2 B+195 A b^3\right )}{5 b \sqrt [4]{a+b x^2}}}{3 b}+\frac {2 x^3 \left (a+b x^2\right )^{3/4} (13 b C-10 a D)}{9 b}}{13 b}+\frac {2 D x^5 \left (a+b x^2\right )^{3/4}}{13 b}\)

Input:

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

Output:

(2*D*x^5*(a + b*x^2)^(3/4))/(13*b) + ((2*(13*b*C - 10*a*D)*x^3*(a + b*x^2) 
^(3/4))/(9*b) + ((2*(39*b^2*B - 26*a*b*C + 20*a^2*D)*x*(a + b*x^2)^(3/4))/ 
(5*b) + ((195*A*b^3 - 78*a*b^2*B + 52*a^2*b*C - 40*a^3*D)*(1 + (b*x^2)/a)^ 
(1/4)*((2*x)/(1 + (b*x^2)/a)^(1/4) - (2*Sqrt[a]*EllipticE[ArcTan[(Sqrt[b]* 
x)/Sqrt[a]]/2, 2])/Sqrt[b]))/(5*b*(a + b*x^2)^(1/4)))/(3*b))/(13*b)
 

Defintions of rubi rules used

rule 27

Int[(a_)*(Fx_), x_Symbol] :> Simp[a   Int[Fx, x], x] /; FreeQ[a, x] &&  !Ma 
tchQ[Fx, (b_)*(Gx_) /; FreeQ[b, x]]
 

rule 212

Int[((a_) + (b_.)*(x_)^2)^(-5/4), x_Symbol] :> Simp[(2/(a^(5/4)*Rt[b/a, 2]) 
)*EllipticE[(1/2)*ArcTan[Rt[b/a, 2]*x], 2], x] /; FreeQ[{a, b}, x] && GtQ[a 
, 0] && PosQ[b/a]
 

rule 225

Int[((a_) + (b_.)*(x_)^2)^(-1/4), x_Symbol] :> Simp[2*(x/(a + b*x^2)^(1/4)) 
, x] - Simp[a   Int[1/(a + b*x^2)^(5/4), x], x] /; FreeQ[{a, b}, x] && GtQ[ 
a, 0] && PosQ[b/a]
 

rule 227

Int[((a_) + (b_.)*(x_)^2)^(-1/4), x_Symbol] :> Simp[(1 + b*(x^2/a))^(1/4)/( 
a + b*x^2)^(1/4)   Int[1/(1 + b*(x^2/a))^(1/4), x], x] /; FreeQ[{a, b}, x] 
&& PosQ[a]
 

rule 299

Int[((a_) + (b_.)*(x_)^2)^(p_)*((c_) + (d_.)*(x_)^2), x_Symbol] :> Simp[d*x 
*((a + b*x^2)^(p + 1)/(b*(2*p + 3))), x] - Simp[(a*d - b*c*(2*p + 3))/(b*(2 
*p + 3))   Int[(a + b*x^2)^p, x], x] /; FreeQ[{a, b, c, d}, x] && NeQ[b*c - 
 a*d, 0] && NeQ[2*p + 3, 0]
 

rule 1473

Int[((d_) + (e_.)*(x_)^2)^(q_)*((a_) + (b_.)*(x_)^2 + (c_.)*(x_)^4)^(p_.), 
x_Symbol] :> Simp[c^p*x^(4*p - 1)*((d + e*x^2)^(q + 1)/(e*(4*p + 2*q + 1))) 
, x] + Simp[1/(e*(4*p + 2*q + 1))   Int[(d + e*x^2)^q*ExpandToSum[e*(4*p + 
2*q + 1)*(a + b*x^2 + c*x^4)^p - d*c^p*(4*p - 1)*x^(4*p - 2) - e*c^p*(4*p + 
 2*q + 1)*x^(4*p), x], x], x] /; FreeQ[{a, b, c, d, e, q}, x] && NeQ[b^2 - 
4*a*c, 0] && NeQ[c*d^2 - b*d*e + a*e^2, 0] && IGtQ[p, 0] &&  !LtQ[q, -1]
 

rule 2346

Int[(Pq_)*((a_) + (b_.)*(x_)^2)^(p_), x_Symbol] :> With[{q = Expon[Pq, x], 
e = Coeff[Pq, x, Expon[Pq, x]]}, Simp[e*x^(q - 1)*((a + b*x^2)^(p + 1)/(b*( 
q + 2*p + 1))), x] + Simp[1/(b*(q + 2*p + 1))   Int[(a + b*x^2)^p*ExpandToS 
um[b*(q + 2*p + 1)*Pq - a*e*(q - 1)*x^(q - 2) - b*e*(q + 2*p + 1)*x^q, x], 
x], x]] /; FreeQ[{a, b, p}, x] && PolyQ[Pq, x] &&  !LeQ[p, -1]
 

Maple [F]

\[\int \frac {D x^{6}+C \,x^{4}+x^{2} B +A}{\left (b \,x^{2}+a \right )^{\frac {1}{4}}}d x\]

Input:

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

Output:

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

Fricas [F]

\[ \int \frac {A+B x^2+C x^4+D x^6}{\sqrt [4]{a+b x^2}} \, dx=\int { \frac {D x^{6} + C x^{4} + B x^{2} + A}{{\left (b x^{2} + a\right )}^{\frac {1}{4}}} \,d x } \] Input:

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

Output:

integral((D*x^6 + C*x^4 + B*x^2 + A)/(b*x^2 + a)^(1/4), x)
 

Sympy [C] (verification not implemented)

Result contains complex when optimal does not.

Time = 1.55 (sec) , antiderivative size = 117, normalized size of antiderivative = 0.53 \[ \int \frac {A+B x^2+C x^4+D x^6}{\sqrt [4]{a+b x^2}} \, dx=\frac {A x {{}_{2}F_{1}\left (\begin {matrix} \frac {1}{4}, \frac {1}{2} \\ \frac {3}{2} \end {matrix}\middle | {\frac {b x^{2} e^{i \pi }}{a}} \right )}}{\sqrt [4]{a}} + \frac {B x^{3} {{}_{2}F_{1}\left (\begin {matrix} \frac {1}{4}, \frac {3}{2} \\ \frac {5}{2} \end {matrix}\middle | {\frac {b x^{2} e^{i \pi }}{a}} \right )}}{3 \sqrt [4]{a}} + \frac {C x^{5} {{}_{2}F_{1}\left (\begin {matrix} \frac {1}{4}, \frac {5}{2} \\ \frac {7}{2} \end {matrix}\middle | {\frac {b x^{2} e^{i \pi }}{a}} \right )}}{5 \sqrt [4]{a}} + \frac {D x^{7} {{}_{2}F_{1}\left (\begin {matrix} \frac {1}{4}, \frac {7}{2} \\ \frac {9}{2} \end {matrix}\middle | {\frac {b x^{2} e^{i \pi }}{a}} \right )}}{7 \sqrt [4]{a}} \] Input:

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

Output:

A*x*hyper((1/4, 1/2), (3/2,), b*x**2*exp_polar(I*pi)/a)/a**(1/4) + B*x**3* 
hyper((1/4, 3/2), (5/2,), b*x**2*exp_polar(I*pi)/a)/(3*a**(1/4)) + C*x**5* 
hyper((1/4, 5/2), (7/2,), b*x**2*exp_polar(I*pi)/a)/(5*a**(1/4)) + D*x**7* 
hyper((1/4, 7/2), (9/2,), b*x**2*exp_polar(I*pi)/a)/(7*a**(1/4))
 

Maxima [F]

\[ \int \frac {A+B x^2+C x^4+D x^6}{\sqrt [4]{a+b x^2}} \, dx=\int { \frac {D x^{6} + C x^{4} + B x^{2} + A}{{\left (b x^{2} + a\right )}^{\frac {1}{4}}} \,d x } \] Input:

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

Output:

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

Giac [F]

\[ \int \frac {A+B x^2+C x^4+D x^6}{\sqrt [4]{a+b x^2}} \, dx=\int { \frac {D x^{6} + C x^{4} + B x^{2} + A}{{\left (b x^{2} + a\right )}^{\frac {1}{4}}} \,d x } \] Input:

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

Output:

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

Mupad [F(-1)]

Timed out. \[ \int \frac {A+B x^2+C x^4+D x^6}{\sqrt [4]{a+b x^2}} \, dx=\int \frac {A+B\,x^2+C\,x^4+x^6\,D}{{\left (b\,x^2+a\right )}^{1/4}} \,d x \] Input:

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

Output:

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

Reduce [F]

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

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

Output:

int(x**6/(a + b*x**2)**(1/4),x)*d + int(x**4/(a + b*x**2)**(1/4),x)*c + in 
t(x**2/(a + b*x**2)**(1/4),x)*b + int(1/(a + b*x**2)**(1/4),x)*a