\(\int \frac {x^2 (c+d x^4)}{\sqrt [4]{a+b x^4}} \, dx\) [82]

Optimal result
Mathematica [C] (verified)
Rubi [A] (verified)
Maple [F]
Fricas [F]
Sympy [C] (verification not implemented)
Maxima [F]
Giac [F]
Mupad [F(-1)]
Reduce [F]

Optimal result

Integrand size = 22, antiderivative size = 123 \[ \int \frac {x^2 \left (c+d x^4\right )}{\sqrt [4]{a+b x^4}} \, dx=\frac {(2 b c-a d) x^3}{4 b \sqrt [4]{a+b x^4}}+\frac {d x^3 \left (a+b x^4\right )^{3/4}}{6 b}+\frac {\sqrt {a} (2 b c-a d) \sqrt [4]{1+\frac {a}{b x^4}} x E\left (\left .\frac {1}{2} \cot ^{-1}\left (\frac {\sqrt {b} x^2}{\sqrt {a}}\right )\right |2\right )}{4 b^{3/2} \sqrt [4]{a+b x^4}} \] Output:

1/4*(-a*d+2*b*c)*x^3/b/(b*x^4+a)^(1/4)+1/6*d*x^3*(b*x^4+a)^(3/4)/b+1/4*a^( 
1/2)*(-a*d+2*b*c)*(1+a/b/x^4)^(1/4)*x*EllipticE(sin(1/2*arccot(b^(1/2)*x^2 
/a^(1/2))),2^(1/2))/b^(3/2)/(b*x^4+a)^(1/4)
                                                                                    
                                                                                    
 

Mathematica [C] (verified)

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

Time = 10.07 (sec) , antiderivative size = 74, normalized size of antiderivative = 0.60 \[ \int \frac {x^2 \left (c+d x^4\right )}{\sqrt [4]{a+b x^4}} \, dx=\frac {x^3 \left (d \left (a+b x^4\right )+(2 b c-a d) \sqrt [4]{1+\frac {b x^4}{a}} \operatorname {Hypergeometric2F1}\left (\frac {1}{4},\frac {3}{4},\frac {7}{4},-\frac {b x^4}{a}\right )\right )}{6 b \sqrt [4]{a+b x^4}} \] Input:

Integrate[(x^2*(c + d*x^4))/(a + b*x^4)^(1/4),x]
 

Output:

(x^3*(d*(a + b*x^4) + (2*b*c - a*d)*(1 + (b*x^4)/a)^(1/4)*Hypergeometric2F 
1[1/4, 3/4, 7/4, -((b*x^4)/a)]))/(6*b*(a + b*x^4)^(1/4))
 

Rubi [A] (verified)

Time = 0.45 (sec) , antiderivative size = 119, normalized size of antiderivative = 0.97, number of steps used = 7, number of rules used = 6, \(\frac {\text {number of rules}}{\text {integrand size}}\) = 0.273, Rules used = {959, 839, 813, 858, 807, 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 {x^2 \left (c+d x^4\right )}{\sqrt [4]{a+b x^4}} \, dx\)

\(\Big \downarrow \) 959

\(\displaystyle \frac {(2 b c-a d) \int \frac {x^2}{\sqrt [4]{b x^4+a}}dx}{2 b}+\frac {d x^3 \left (a+b x^4\right )^{3/4}}{6 b}\)

\(\Big \downarrow \) 839

\(\displaystyle \frac {(2 b c-a d) \left (\frac {x^3}{2 \sqrt [4]{a+b x^4}}-\frac {1}{2} a \int \frac {x^2}{\left (b x^4+a\right )^{5/4}}dx\right )}{2 b}+\frac {d x^3 \left (a+b x^4\right )^{3/4}}{6 b}\)

\(\Big \downarrow \) 813

\(\displaystyle \frac {(2 b c-a d) \left (\frac {x^3}{2 \sqrt [4]{a+b x^4}}-\frac {a x \sqrt [4]{\frac {a}{b x^4}+1} \int \frac {1}{\left (\frac {a}{b x^4}+1\right )^{5/4} x^3}dx}{2 b \sqrt [4]{a+b x^4}}\right )}{2 b}+\frac {d x^3 \left (a+b x^4\right )^{3/4}}{6 b}\)

\(\Big \downarrow \) 858

\(\displaystyle \frac {(2 b c-a d) \left (\frac {a x \sqrt [4]{\frac {a}{b x^4}+1} \int \frac {1}{\left (\frac {a}{b x^4}+1\right )^{5/4} x}d\frac {1}{x}}{2 b \sqrt [4]{a+b x^4}}+\frac {x^3}{2 \sqrt [4]{a+b x^4}}\right )}{2 b}+\frac {d x^3 \left (a+b x^4\right )^{3/4}}{6 b}\)

\(\Big \downarrow \) 807

\(\displaystyle \frac {(2 b c-a d) \left (\frac {a x \sqrt [4]{\frac {a}{b x^4}+1} \int \frac {1}{\left (\frac {a}{b x^2}+1\right )^{5/4}}d\frac {1}{x^2}}{4 b \sqrt [4]{a+b x^4}}+\frac {x^3}{2 \sqrt [4]{a+b x^4}}\right )}{2 b}+\frac {d x^3 \left (a+b x^4\right )^{3/4}}{6 b}\)

\(\Big \downarrow \) 212

\(\displaystyle \frac {(2 b c-a d) \left (\frac {\sqrt {a} x \sqrt [4]{\frac {a}{b x^4}+1} E\left (\left .\frac {1}{2} \arctan \left (\frac {\sqrt {a}}{\sqrt {b} x^2}\right )\right |2\right )}{2 \sqrt {b} \sqrt [4]{a+b x^4}}+\frac {x^3}{2 \sqrt [4]{a+b x^4}}\right )}{2 b}+\frac {d x^3 \left (a+b x^4\right )^{3/4}}{6 b}\)

Input:

Int[(x^2*(c + d*x^4))/(a + b*x^4)^(1/4),x]
 

Output:

(d*x^3*(a + b*x^4)^(3/4))/(6*b) + ((2*b*c - a*d)*(x^3/(2*(a + b*x^4)^(1/4) 
) + (Sqrt[a]*(1 + a/(b*x^4))^(1/4)*x*EllipticE[ArcTan[Sqrt[a]/(Sqrt[b]*x^2 
)]/2, 2])/(2*Sqrt[b]*(a + b*x^4)^(1/4))))/(2*b)
 

Defintions of rubi rules used

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 807
Int[(x_)^(m_.)*((a_) + (b_.)*(x_)^(n_))^(p_), x_Symbol] :> With[{k = GCD[m 
+ 1, n]}, Simp[1/k   Subst[Int[x^((m + 1)/k - 1)*(a + b*x^(n/k))^p, x], x, 
x^k], x] /; k != 1] /; FreeQ[{a, b, p}, x] && IGtQ[n, 0] && IntegerQ[m]
 

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

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

rule 858
Int[(x_)^(m_.)*((a_) + (b_.)*(x_)^(n_))^(p_), x_Symbol] :> -Subst[Int[(a + 
b/x^n)^p/x^(m + 2), x], x, 1/x] /; FreeQ[{a, b, p}, x] && ILtQ[n, 0] && Int 
egerQ[m]
 

rule 959
Int[((e_.)*(x_))^(m_.)*((a_) + (b_.)*(x_)^(n_))^(p_.)*((c_) + (d_.)*(x_)^(n 
_)), x_Symbol] :> Simp[d*(e*x)^(m + 1)*((a + b*x^n)^(p + 1)/(b*e*(m + n*(p 
+ 1) + 1))), x] - Simp[(a*d*(m + 1) - b*c*(m + n*(p + 1) + 1))/(b*(m + n*(p 
 + 1) + 1))   Int[(e*x)^m*(a + b*x^n)^p, x], x] /; FreeQ[{a, b, c, d, e, m, 
 n, p}, x] && NeQ[b*c - a*d, 0] && NeQ[m + n*(p + 1) + 1, 0]
 
Maple [F]

\[\int \frac {x^{2} \left (d \,x^{4}+c \right )}{\left (b \,x^{4}+a \right )^{\frac {1}{4}}}d x\]

Input:

int(x^2*(d*x^4+c)/(b*x^4+a)^(1/4),x)
 

Output:

int(x^2*(d*x^4+c)/(b*x^4+a)^(1/4),x)
 

Fricas [F]

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

integrate(x^2*(d*x^4+c)/(b*x^4+a)^(1/4),x, algorithm="fricas")
 

Output:

integral((d*x^6 + c*x^2)/(b*x^4 + a)^(1/4), x)
 

Sympy [C] (verification not implemented)

Result contains complex when optimal does not.

Time = 3.44 (sec) , antiderivative size = 80, normalized size of antiderivative = 0.65 \[ \int \frac {x^2 \left (c+d x^4\right )}{\sqrt [4]{a+b x^4}} \, dx=\frac {c x^{3} \Gamma \left (\frac {3}{4}\right ) {{}_{2}F_{1}\left (\begin {matrix} \frac {1}{4}, \frac {3}{4} \\ \frac {7}{4} \end {matrix}\middle | {\frac {b x^{4} e^{i \pi }}{a}} \right )}}{4 \sqrt [4]{a} \Gamma \left (\frac {7}{4}\right )} + \frac {d x^{7} \Gamma \left (\frac {7}{4}\right ) {{}_{2}F_{1}\left (\begin {matrix} \frac {1}{4}, \frac {7}{4} \\ \frac {11}{4} \end {matrix}\middle | {\frac {b x^{4} e^{i \pi }}{a}} \right )}}{4 \sqrt [4]{a} \Gamma \left (\frac {11}{4}\right )} \] Input:

integrate(x**2*(d*x**4+c)/(b*x**4+a)**(1/4),x)
 

Output:

c*x**3*gamma(3/4)*hyper((1/4, 3/4), (7/4,), b*x**4*exp_polar(I*pi)/a)/(4*a 
**(1/4)*gamma(7/4)) + d*x**7*gamma(7/4)*hyper((1/4, 7/4), (11/4,), b*x**4* 
exp_polar(I*pi)/a)/(4*a**(1/4)*gamma(11/4))
 

Maxima [F]

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

integrate(x^2*(d*x^4+c)/(b*x^4+a)^(1/4),x, algorithm="maxima")
 

Output:

integrate((d*x^4 + c)*x^2/(b*x^4 + a)^(1/4), x)
 

Giac [F]

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

integrate(x^2*(d*x^4+c)/(b*x^4+a)^(1/4),x, algorithm="giac")
 

Output:

integrate((d*x^4 + c)*x^2/(b*x^4 + a)^(1/4), x)
 

Mupad [F(-1)]

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

int((x^2*(c + d*x^4))/(a + b*x^4)^(1/4),x)
 

Output:

int((x^2*(c + d*x^4))/(a + b*x^4)^(1/4), x)
 

Reduce [F]

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

int(x^2*(d*x^4+c)/(b*x^4+a)^(1/4),x)
 

Output:

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