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

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

Optimal result

Integrand size = 17, antiderivative size = 136 \[ \int \frac {x^3}{\sqrt [4]{a x+b x^2}} \, dx=\frac {11 a^2 \left (a x+b x^2\right )^{3/4}}{30 b^3}-\frac {11 a x \left (a x+b x^2\right )^{3/4}}{35 b^2}+\frac {2 x^2 \left (a x+b x^2\right )^{3/4}}{7 b}-\frac {11 a^4 \sqrt [4]{-\frac {b x}{a}-\frac {b^2 x^2}{a^2}} E\left (\left .\frac {1}{2} \arcsin \left (1+\frac {2 b x}{a}\right )\right |2\right )}{20 \sqrt {2} b^4 \sqrt [4]{a x+b x^2}} \] Output:

11/30*a^2*(b*x^2+a*x)^(3/4)/b^3-11/35*a*x*(b*x^2+a*x)^(3/4)/b^2+2/7*x^2*(b 
*x^2+a*x)^(3/4)/b-11/40*a^4*(-b*x/a-b^2*x^2/a^2)^(1/4)*EllipticE(sin(1/2*a 
rcsin(1+2*b*x/a)),2^(1/2))*2^(1/2)/b^4/(b*x^2+a*x)^(1/4)
                                                                                    
                                                                                    
 

Mathematica [C] (verified)

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

Time = 10.01 (sec) , antiderivative size = 47, normalized size of antiderivative = 0.35 \[ \int \frac {x^3}{\sqrt [4]{a x+b x^2}} \, dx=\frac {4 x^4 \sqrt [4]{1+\frac {b x}{a}} \operatorname {Hypergeometric2F1}\left (\frac {1}{4},\frac {15}{4},\frac {19}{4},-\frac {b x}{a}\right )}{15 \sqrt [4]{x (a+b x)}} \] Input:

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

Output:

(4*x^4*(1 + (b*x)/a)^(1/4)*Hypergeometric2F1[1/4, 15/4, 19/4, -((b*x)/a)]) 
/(15*(x*(a + b*x))^(1/4))
 

Rubi [A] (warning: unable to verify)

Time = 0.56 (sec) , antiderivative size = 200, normalized size of antiderivative = 1.47, number of steps used = 11, number of rules used = 10, \(\frac {\text {number of rules}}{\text {integrand size}}\) = 0.588, Rules used = {1137, 60, 60, 60, 73, 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^3}{\sqrt [4]{a x+b x^2}} \, dx\)

\(\Big \downarrow \) 1137

\(\displaystyle \frac {\sqrt [4]{x} \sqrt [4]{a+b x} \int \frac {x^{11/4}}{\sqrt [4]{a+b x}}dx}{\sqrt [4]{a x+b x^2}}\)

\(\Big \downarrow \) 60

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

\(\Big \downarrow \) 60

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

\(\Big \downarrow \) 60

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

\(\Big \downarrow \) 73

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

\(\Big \downarrow \) 839

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

\(\Big \downarrow \) 813

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

\(\Big \downarrow \) 858

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

\(\Big \downarrow \) 807

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

\(\Big \downarrow \) 212

\(\displaystyle \frac {\sqrt [4]{x} \sqrt [4]{a+b x} \left (\frac {2 x^{11/4} (a+b x)^{3/4}}{7 b}-\frac {11 a \left (\frac {2 x^{7/4} (a+b x)^{3/4}}{5 b}-\frac {7 a \left (\frac {2 x^{3/4} (a+b x)^{3/4}}{3 b}-\frac {2 a \left (\frac {\sqrt {a} \sqrt [4]{x} \sqrt [4]{\frac {a}{b x}+1} E\left (\left .\frac {1}{2} \arctan \left (\frac {\sqrt {a} \sqrt {x}}{\sqrt {b}}\right )\right |2\right )}{2 \sqrt {b} \sqrt [4]{a+b x}}+\frac {x^{3/4}}{2 \sqrt [4]{a+b x}}\right )}{b}\right )}{10 b}\right )}{14 b}\right )}{\sqrt [4]{a x+b x^2}}\)

Input:

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

Output:

(x^(1/4)*(a + b*x)^(1/4)*((2*x^(11/4)*(a + b*x)^(3/4))/(7*b) - (11*a*((2*x 
^(7/4)*(a + b*x)^(3/4))/(5*b) - (7*a*((2*x^(3/4)*(a + b*x)^(3/4))/(3*b) - 
(2*a*(x^(3/4)/(2*(a + b*x)^(1/4)) + (Sqrt[a]*(1 + a/(b*x))^(1/4)*x^(1/4)*E 
llipticE[ArcTan[(Sqrt[a]*Sqrt[x])/Sqrt[b]]/2, 2])/(2*Sqrt[b]*(a + b*x)^(1/ 
4))))/b))/(10*b)))/(14*b)))/(a*x + b*x^2)^(1/4)
 

Defintions of rubi rules used

rule 60
Int[((a_.) + (b_.)*(x_))^(m_)*((c_.) + (d_.)*(x_))^(n_), x_Symbol] :> Simp[ 
(a + b*x)^(m + 1)*((c + d*x)^n/(b*(m + n + 1))), x] + Simp[n*((b*c - a*d)/( 
b*(m + n + 1)))   Int[(a + b*x)^m*(c + d*x)^(n - 1), x], x] /; FreeQ[{a, b, 
 c, d}, x] && GtQ[n, 0] && NeQ[m + n + 1, 0] &&  !(IGtQ[m, 0] && ( !Integer 
Q[n] || (GtQ[m, 0] && LtQ[m - n, 0]))) &&  !ILtQ[m + n + 2, 0] && IntLinear 
Q[a, b, c, d, m, n, 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 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 1137
Int[((e_.)*(x_))^(m_)*((b_.)*(x_) + (c_.)*(x_)^2)^(p_), x_Symbol] :> Simp[( 
e*x)^m*((b*x + c*x^2)^p/(x^(m + p)*(b + c*x)^p))   Int[x^(m + p)*(b + c*x)^ 
p, x], x] /; FreeQ[{b, c, e, m}, x]
 
Maple [F]

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

Input:

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

Output:

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

Fricas [F]

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

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

Output:

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

Sympy [F]

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

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

Output:

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

Maxima [F]

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

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

Output:

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

Giac [F]

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

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

Output:

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

Mupad [F(-1)]

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

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

Output:

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

Reduce [F]

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

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

Output:

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