\(\int \frac {x^2}{\sqrt [4]{a x^2+b x^3}} \, dx\) [387]

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

Optimal result

Integrand size = 19, antiderivative size = 145 \[ \int \frac {x^2}{\sqrt [4]{a x^2+b x^3}} \, dx=\frac {8 a^2 x}{15 b^2 \sqrt [4]{a x^2+b x^3}}-\frac {4 a x^2}{45 b \sqrt [4]{a x^2+b x^3}}+\frac {4 x^3}{9 \sqrt [4]{a x^2+b x^3}}-\frac {16 a^{5/2} \sqrt {x} \sqrt [4]{\frac {a+b x}{a}} E\left (\left .\frac {1}{2} \arctan \left (\frac {\sqrt {b} \sqrt {x}}{\sqrt {a}}\right )\right |2\right )}{15 b^{5/2} \sqrt [4]{a x^2+b x^3}} \] Output:

8/15*a^2*x/b^2/(b*x^3+a*x^2)^(1/4)-4/45*a*x^2/b/(b*x^3+a*x^2)^(1/4)+4/9*x^ 
3/(b*x^3+a*x^2)^(1/4)-16/15*a^(5/2)*x^(1/2)*((b*x+a)/a)^(1/4)*EllipticE(si 
n(1/2*arctan(b^(1/2)*x^(1/2)/a^(1/2))),2^(1/2))/b^(5/2)/(b*x^3+a*x^2)^(1/4 
)
                                                                                    
                                                                                    
 

Mathematica [C] (verified)

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

Time = 10.02 (sec) , antiderivative size = 49, normalized size of antiderivative = 0.34 \[ \int \frac {x^2}{\sqrt [4]{a x^2+b x^3}} \, dx=\frac {2 x^3 \sqrt [4]{1+\frac {b x}{a}} \operatorname {Hypergeometric2F1}\left (\frac {1}{4},\frac {5}{2},\frac {7}{2},-\frac {b x}{a}\right )}{5 \sqrt [4]{x^2 (a+b x)}} \] Input:

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

Output:

(2*x^3*(1 + (b*x)/a)^(1/4)*Hypergeometric2F1[1/4, 5/2, 7/2, -((b*x)/a)])/( 
5*(x^2*(a + b*x))^(1/4))
 

Rubi [A] (verified)

Time = 0.74 (sec) , antiderivative size = 216, normalized size of antiderivative = 1.49, number of steps used = 12, number of rules used = 11, \(\frac {\text {number of rules}}{\text {integrand size}}\) = 0.579, Rules used = {1930, 1930, 1917, 73, 836, 27, 765, 762, 1390, 1389, 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 {x^2}{\sqrt [4]{a x^2+b x^3}} \, dx\)

\(\Big \downarrow \) 1930

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

\(\Big \downarrow \) 1930

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

\(\Big \downarrow \) 1917

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

\(\Big \downarrow \) 73

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

\(\Big \downarrow \) 836

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

\(\Big \downarrow \) 27

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

\(\Big \downarrow \) 765

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

\(\Big \downarrow \) 762

\(\displaystyle \frac {4 \left (a x^2+b x^3\right )^{3/4}}{9 b}-\frac {2 a \left (\frac {4 \left (a x^2+b x^3\right )^{3/4}}{5 b x}-\frac {8 a \sqrt {x} \sqrt [4]{a+b x} \left (\int \frac {\sqrt {a}+\sqrt {a+b x}}{\sqrt {\frac {a+b x}{b}-\frac {a}{b}}}d\sqrt [4]{a+b x}-\frac {a^{3/4} \sqrt {1-\frac {a+b x}{a}} \operatorname {EllipticF}\left (\arcsin \left (\frac {\sqrt [4]{a+b x}}{\sqrt [4]{a}}\right ),-1\right )}{\sqrt {\frac {a+b x}{b}-\frac {a}{b}}}\right )}{5 b^2 \sqrt [4]{a x^2+b x^3}}\right )}{3 b}\)

\(\Big \downarrow \) 1390

\(\displaystyle \frac {4 \left (a x^2+b x^3\right )^{3/4}}{9 b}-\frac {2 a \left (\frac {4 \left (a x^2+b x^3\right )^{3/4}}{5 b x}-\frac {8 a \sqrt {x} \sqrt [4]{a+b x} \left (\frac {\sqrt {1-\frac {a+b x}{a}} \int \frac {\sqrt {a}+\sqrt {a+b x}}{\sqrt {1-\frac {a+b x}{a}}}d\sqrt [4]{a+b x}}{\sqrt {\frac {a+b x}{b}-\frac {a}{b}}}-\frac {a^{3/4} \sqrt {1-\frac {a+b x}{a}} \operatorname {EllipticF}\left (\arcsin \left (\frac {\sqrt [4]{a+b x}}{\sqrt [4]{a}}\right ),-1\right )}{\sqrt {\frac {a+b x}{b}-\frac {a}{b}}}\right )}{5 b^2 \sqrt [4]{a x^2+b x^3}}\right )}{3 b}\)

\(\Big \downarrow \) 1389

\(\displaystyle \frac {4 \left (a x^2+b x^3\right )^{3/4}}{9 b}-\frac {2 a \left (\frac {4 \left (a x^2+b x^3\right )^{3/4}}{5 b x}-\frac {8 a \sqrt {x} \sqrt [4]{a+b x} \left (\frac {\sqrt {a} \sqrt {1-\frac {a+b x}{a}} \int \frac {\sqrt {\frac {\sqrt {a+b x}}{\sqrt {a}}+1}}{\sqrt {1-\frac {\sqrt {a+b x}}{\sqrt {a}}}}d\sqrt [4]{a+b x}}{\sqrt {\frac {a+b x}{b}-\frac {a}{b}}}-\frac {a^{3/4} \sqrt {1-\frac {a+b x}{a}} \operatorname {EllipticF}\left (\arcsin \left (\frac {\sqrt [4]{a+b x}}{\sqrt [4]{a}}\right ),-1\right )}{\sqrt {\frac {a+b x}{b}-\frac {a}{b}}}\right )}{5 b^2 \sqrt [4]{a x^2+b x^3}}\right )}{3 b}\)

\(\Big \downarrow \) 327

\(\displaystyle \frac {4 \left (a x^2+b x^3\right )^{3/4}}{9 b}-\frac {2 a \left (\frac {4 \left (a x^2+b x^3\right )^{3/4}}{5 b x}-\frac {8 a \sqrt {x} \sqrt [4]{a+b x} \left (\frac {a^{3/4} \sqrt {1-\frac {a+b x}{a}} E\left (\left .\arcsin \left (\frac {\sqrt [4]{a+b x}}{\sqrt [4]{a}}\right )\right |-1\right )}{\sqrt {\frac {a+b x}{b}-\frac {a}{b}}}-\frac {a^{3/4} \sqrt {1-\frac {a+b x}{a}} \operatorname {EllipticF}\left (\arcsin \left (\frac {\sqrt [4]{a+b x}}{\sqrt [4]{a}}\right ),-1\right )}{\sqrt {\frac {a+b x}{b}-\frac {a}{b}}}\right )}{5 b^2 \sqrt [4]{a x^2+b x^3}}\right )}{3 b}\)

Input:

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

Output:

(4*(a*x^2 + b*x^3)^(3/4))/(9*b) - (2*a*((4*(a*x^2 + b*x^3)^(3/4))/(5*b*x) 
- (8*a*Sqrt[x]*(a + b*x)^(1/4)*((a^(3/4)*Sqrt[1 - (a + b*x)/a]*EllipticE[A 
rcSin[(a + b*x)^(1/4)/a^(1/4)], -1])/Sqrt[-(a/b) + (a + b*x)/b] - (a^(3/4) 
*Sqrt[1 - (a + b*x)/a]*EllipticF[ArcSin[(a + b*x)^(1/4)/a^(1/4)], -1])/Sqr 
t[-(a/b) + (a + b*x)/b]))/(5*b^2*(a*x^2 + b*x^3)^(1/4))))/(3*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 73
Int[((a_.) + (b_.)*(x_))^(m_)*((c_.) + (d_.)*(x_))^(n_), x_Symbol] :> With[ 
{p = Denominator[m]}, Simp[p/b   Subst[Int[x^(p*(m + 1) - 1)*(c - a*(d/b) + 
 d*(x^p/b))^n, x], x, (a + b*x)^(1/p)], x]] /; FreeQ[{a, b, c, d}, x] && Lt 
Q[-1, m, 0] && LeQ[-1, n, 0] && LeQ[Denominator[n], Denominator[m]] && IntL 
inearQ[a, b, c, d, m, n, x]
 

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 762
Int[1/Sqrt[(a_) + (b_.)*(x_)^4], x_Symbol] :> Simp[(1/(Sqrt[a]*Rt[-b/a, 4]) 
)*EllipticF[ArcSin[Rt[-b/a, 4]*x], -1], x] /; FreeQ[{a, b}, x] && NegQ[b/a] 
 && GtQ[a, 0]
 

rule 765
Int[1/Sqrt[(a_) + (b_.)*(x_)^4], x_Symbol] :> Simp[Sqrt[1 + b*(x^4/a)]/Sqrt 
[a + b*x^4]   Int[1/Sqrt[1 + b*(x^4/a)], x], x] /; FreeQ[{a, b}, x] && NegQ 
[b/a] &&  !GtQ[a, 0]
 

rule 836
Int[(x_)^2/Sqrt[(a_) + (b_.)*(x_)^4], x_Symbol] :> With[{q = Rt[-b/a, 2]}, 
Simp[-q^(-1)   Int[1/Sqrt[a + b*x^4], x], x] + Simp[1/q   Int[(1 + q*x^2)/S 
qrt[a + b*x^4], x], x]] /; FreeQ[{a, b}, x] && NegQ[b/a]
 

rule 1389
Int[((d_) + (e_.)*(x_)^2)/Sqrt[(a_) + (c_.)*(x_)^4], x_Symbol] :> Simp[d/Sq 
rt[a]   Int[Sqrt[1 + e*(x^2/d)]/Sqrt[1 - e*(x^2/d)], x], x] /; FreeQ[{a, c, 
 d, e}, x] && EqQ[c*d^2 + a*e^2, 0] && NegQ[c/a] && GtQ[a, 0]
 

rule 1390
Int[((d_) + (e_.)*(x_)^2)/Sqrt[(a_) + (c_.)*(x_)^4], x_Symbol] :> Simp[Sqrt 
[1 + c*(x^4/a)]/Sqrt[a + c*x^4]   Int[(d + e*x^2)/Sqrt[1 + c*(x^4/a)], x], 
x] /; FreeQ[{a, c, d, e}, x] && EqQ[c*d^2 + a*e^2, 0] && NegQ[c/a] &&  !GtQ 
[a, 0] &&  !(LtQ[a, 0] && GtQ[c, 0])
 

rule 1917
Int[((a_.)*(x_)^(j_.) + (b_.)*(x_)^(n_.))^(p_), x_Symbol] :> Simp[(a*x^j + 
b*x^n)^FracPart[p]/(x^(j*FracPart[p])*(a + b*x^(n - j))^FracPart[p])   Int[ 
x^(j*p)*(a + b*x^(n - j))^p, x], x] /; FreeQ[{a, b, j, n, p}, x] &&  !Integ 
erQ[p] && NeQ[n, j] && PosQ[n - j]
 

rule 1930
Int[((c_.)*(x_))^(m_.)*((a_.)*(x_)^(j_.) + (b_.)*(x_)^(n_.))^(p_), x_Symbol 
] :> Simp[c^(n - 1)*(c*x)^(m - n + 1)*((a*x^j + b*x^n)^(p + 1)/(b*(m + n*p 
+ 1))), x] - Simp[a*c^(n - j)*((m + j*p - n + j + 1)/(b*(m + n*p + 1)))   I 
nt[(c*x)^(m - (n - j))*(a*x^j + b*x^n)^p, x], x] /; FreeQ[{a, b, c, m, p}, 
x] &&  !IntegerQ[p] && LtQ[0, j, n] && (IntegersQ[j, n] || GtQ[c, 0]) && Gt 
Q[m + j*p - n + j + 1, 0] && NeQ[m + n*p + 1, 0]
 
Maple [F]

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

Input:

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

Output:

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

Fricas [F]

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

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

Output:

integral((b*x^3 + a*x^2)^(3/4)/(b*x + a), x)
 

Sympy [F]

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

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

Output:

Integral(x**2/(x**2*(a + b*x))**(1/4), x)
 

Maxima [F]

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

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

Output:

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

Giac [F]

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

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

Output:

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

Mupad [F(-1)]

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

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

Output:

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

Reduce [F]

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

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

Output:

(4*( - 2*sqrt(x)*(a + b*x)**(1/4)*a**2 + sqrt(x)*(a + b*x)**(1/4)*a*b*x + 
7*sqrt(x)*(a + b*x)**(1/4)*b**2*x**2 + sqrt(a + b*x)*int((sqrt(x)*(a + b*x 
)**(3/4))/(a**2*x + 2*a*b*x**2 + b**2*x**3),x)*a**3))/(77*sqrt(a + b*x)*b* 
*2)