\(\int \frac {1}{(a x+b x^2)^{4/3}} \, dx\) [56]

Optimal result
Mathematica [A] (verified)
Rubi [B] (warning: unable to verify)
Maple [F]
Fricas [F]
Sympy [F]
Maxima [F]
Giac [F]
Mupad [B] (verification not implemented)
Reduce [F]

Optimal result

Integrand size = 13, antiderivative size = 33 \[ \int \frac {1}{\left (a x+b x^2\right )^{4/3}} \, dx=-\frac {3 \operatorname {Hypergeometric2F1}\left (-\frac {2}{3},1,\frac {2}{3},-\frac {b x}{a}\right )}{a \sqrt [3]{a x+b x^2}} \] Output:

-3*hypergeom([-2/3, 1],[2/3],-b*x/a)/a/(b*x^2+a*x)^(1/3)
 

Mathematica [A] (verified)

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

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

Output:

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

Rubi [B] (warning: unable to verify)

Leaf count is larger than twice the leaf count of optimal. \(684\) vs. \(2(33)=66\).

Time = 0.80 (sec) , antiderivative size = 684, normalized size of antiderivative = 20.73, number of steps used = 8, number of rules used = 7, \(\frac {\text {number of rules}}{\text {integrand size}}\) = 0.538, Rules used = {1089, 1093, 1090, 233, 833, 760, 2418}

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 {1}{\left (a x+b x^2\right )^{4/3}} \, dx\)

\(\Big \downarrow \) 1089

\(\displaystyle \frac {2 b \int \frac {1}{\sqrt [3]{b x^2+a x}}dx}{a^2}-\frac {3 (a+2 b x)}{a^2 \sqrt [3]{a x+b x^2}}\)

\(\Big \downarrow \) 1093

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

\(\Big \downarrow \) 1090

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

\(\Big \downarrow \) 233

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

\(\Big \downarrow \) 833

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

\(\Big \downarrow \) 760

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

\(\Big \downarrow \) 2418

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

Input:

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

Output:

(-3*(a + 2*b*x))/(a^2*(a*x + b*x^2)^(1/3)) + (3*b*Sqrt[-((a^2*(-(b/a) - (2 
*b^2*x)/a^2)^2)/b^2)]*(-((b*(a*x + b*x^2))/a^2))^(1/3)*((-2*Sqrt[-((a^2*(- 
(b/a) - (2*b^2*x)/a^2)^2)/b^2)])/(1 - Sqrt[3] + b/a + (2*b^2*x)/a^2) + (3^ 
(1/4)*Sqrt[2 + Sqrt[3]]*(1 + b/a + (2*b^2*x)/a^2)*Sqrt[(1 + (-(b/a) - (2*b 
^2*x)/a^2)^2 + (1 - (a^2*(-(b/a) - (2*b^2*x)/a^2)^2)/b^2)^(1/3))/(1 - Sqrt 
[3] + b/a + (2*b^2*x)/a^2)^2]*EllipticE[ArcSin[(1 + Sqrt[3] + b/a + (2*b^2 
*x)/a^2)/(1 - Sqrt[3] + b/a + (2*b^2*x)/a^2)], -7 + 4*Sqrt[3]])/(Sqrt[-((a 
^2*(-(b/a) - (2*b^2*x)/a^2)^2)/b^2)]*Sqrt[-((1 + b/a + (2*b^2*x)/a^2)/(1 - 
 Sqrt[3] + b/a + (2*b^2*x)/a^2)^2)]) - (2*Sqrt[2 - Sqrt[3]]*(1 + Sqrt[3])* 
(1 + b/a + (2*b^2*x)/a^2)*Sqrt[(1 + (-(b/a) - (2*b^2*x)/a^2)^2 + (1 - (a^2 
*(-(b/a) - (2*b^2*x)/a^2)^2)/b^2)^(1/3))/(1 - Sqrt[3] + b/a + (2*b^2*x)/a^ 
2)^2]*EllipticF[ArcSin[(1 + Sqrt[3] + b/a + (2*b^2*x)/a^2)/(1 - Sqrt[3] + 
b/a + (2*b^2*x)/a^2)], -7 + 4*Sqrt[3]])/(3^(1/4)*Sqrt[-((a^2*(-(b/a) - (2* 
b^2*x)/a^2)^2)/b^2)]*Sqrt[-((1 + b/a + (2*b^2*x)/a^2)/(1 - Sqrt[3] + b/a + 
 (2*b^2*x)/a^2)^2)])))/(2^(1/3)*a^2*(-(b/a) - (2*b^2*x)/a^2)*(a*x + b*x^2) 
^(1/3))
 

Defintions of rubi rules used

rule 233
Int[((a_) + (b_.)*(x_)^2)^(-1/3), x_Symbol] :> Simp[3*(Sqrt[b*x^2]/(2*b*x)) 
   Subst[Int[x/Sqrt[-a + x^3], x], x, (a + b*x^2)^(1/3)], x] /; FreeQ[{a, b 
}, x]
 

rule 760
Int[1/Sqrt[(a_) + (b_.)*(x_)^3], x_Symbol] :> With[{r = Numer[Rt[b/a, 3]], 
s = Denom[Rt[b/a, 3]]}, Simp[2*Sqrt[2 - Sqrt[3]]*(s + r*x)*(Sqrt[(s^2 - r*s 
*x + r^2*x^2)/((1 - Sqrt[3])*s + r*x)^2]/(3^(1/4)*r*Sqrt[a + b*x^3]*Sqrt[(- 
s)*((s + r*x)/((1 - Sqrt[3])*s + r*x)^2)]))*EllipticF[ArcSin[((1 + Sqrt[3]) 
*s + r*x)/((1 - Sqrt[3])*s + r*x)], -7 + 4*Sqrt[3]], x]] /; FreeQ[{a, b}, x 
] && NegQ[a]
 

rule 833
Int[(x_)/Sqrt[(a_) + (b_.)*(x_)^3], x_Symbol] :> With[{r = Numer[Rt[b/a, 3] 
], s = Denom[Rt[b/a, 3]]}, Simp[(-(1 + Sqrt[3]))*(s/r)   Int[1/Sqrt[a + b*x 
^3], x], x] + Simp[1/r   Int[((1 + Sqrt[3])*s + r*x)/Sqrt[a + b*x^3], x], x 
]] /; FreeQ[{a, b}, x] && NegQ[a]
 

rule 1089
Int[((a_.) + (b_.)*(x_) + (c_.)*(x_)^2)^(p_), x_Symbol] :> Simp[(b + 2*c*x) 
*((a + b*x + c*x^2)^(p + 1)/((p + 1)*(b^2 - 4*a*c))), x] - Simp[2*c*((2*p + 
 3)/((p + 1)*(b^2 - 4*a*c)))   Int[(a + b*x + c*x^2)^(p + 1), x], x] /; Fre 
eQ[{a, b, c}, x] && LtQ[p, -1] && (IntegerQ[4*p] || IntegerQ[3*p])
 

rule 1090
Int[((a_.) + (b_.)*(x_) + (c_.)*(x_)^2)^(p_), x_Symbol] :> Simp[1/(2*c*(-4* 
(c/(b^2 - 4*a*c)))^p)   Subst[Int[Simp[1 - x^2/(b^2 - 4*a*c), x]^p, x], x, 
b + 2*c*x], x] /; FreeQ[{a, b, c, p}, x] && GtQ[4*a - b^2/c, 0]
 

rule 1093
Int[((b_.)*(x_) + (c_.)*(x_)^2)^(p_), x_Symbol] :> Simp[(b*x + c*x^2)^p/((- 
c)*((b*x + c*x^2)/b^2))^p   Int[((-c)*(x/b) - c^2*(x^2/b^2))^p, x], x] /; F 
reeQ[{b, c}, x] && (IntegerQ[4*p] || IntegerQ[3*p])
 

rule 2418
Int[((c_) + (d_.)*(x_))/Sqrt[(a_) + (b_.)*(x_)^3], x_Symbol] :> With[{r = N 
umer[Simplify[(1 + Sqrt[3])*(d/c)]], s = Denom[Simplify[(1 + Sqrt[3])*(d/c) 
]]}, Simp[2*d*s^3*(Sqrt[a + b*x^3]/(a*r^2*((1 - Sqrt[3])*s + r*x))), x] + S 
imp[3^(1/4)*Sqrt[2 + Sqrt[3]]*d*s*(s + r*x)*(Sqrt[(s^2 - r*s*x + r^2*x^2)/( 
(1 - Sqrt[3])*s + r*x)^2]/(r^2*Sqrt[a + b*x^3]*Sqrt[(-s)*((s + r*x)/((1 - S 
qrt[3])*s + r*x)^2)]))*EllipticE[ArcSin[((1 + Sqrt[3])*s + r*x)/((1 - Sqrt[ 
3])*s + r*x)], -7 + 4*Sqrt[3]], x]] /; FreeQ[{a, b, c, d}, x] && NegQ[a] && 
 EqQ[b*c^3 - 2*(5 + 3*Sqrt[3])*a*d^3, 0]
 
Maple [F]

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

Input:

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

Output:

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

Fricas [F]

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

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

Output:

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

Sympy [F]

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

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

Output:

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

Maxima [F]

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

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

Output:

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

Giac [F]

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

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

Output:

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

Mupad [B] (verification not implemented)

Time = 9.44 (sec) , antiderivative size = 36, normalized size of antiderivative = 1.09 \[ \int \frac {1}{\left (a x+b x^2\right )^{4/3}} \, dx=-\frac {3\,x\,{\left (\frac {b\,x}{a}+1\right )}^{4/3}\,{{}}_2{\mathrm {F}}_1\left (-\frac {1}{3},\frac {4}{3};\ \frac {2}{3};\ -\frac {b\,x}{a}\right )}{{\left (b\,x^2+a\,x\right )}^{4/3}} \] Input:

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

Output:

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

Reduce [F]

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

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

Output:

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