3.2.36 \(\int \frac {x^4}{27 a^3+27 a^2 b x^2+27 a^2 c x^3+9 a b^2 x^4+b^3 x^6} \, dx\) [136]

3.2.36.1 Optimal result
3.2.36.2 Mathematica [C] (verified)
3.2.36.3 Rubi [A] (verified)
3.2.36.4 Maple [C] (verified)
3.2.36.5 Fricas [F(-2)]
3.2.36.6 Sympy [F(-1)]
3.2.36.7 Maxima [F]
3.2.36.8 Giac [F]
3.2.36.9 Mupad [B] (verification not implemented)

3.2.36.1 Optimal result

Integrand size = 46, antiderivative size = 545 \[ \int \frac {x^4}{27 a^3+27 a^2 b x^2+27 a^2 c x^3+9 a b^2 x^4+b^3 x^6} \, dx=-\frac {\sqrt [3]{-1} \left (2 \sqrt [3]{-1} b+3 \sqrt [3]{a} c^{2/3}\right ) \arctan \left (\frac {3 \sqrt [3]{-1} a^{2/3} \sqrt [3]{c}-2 b x}{\sqrt {3} \sqrt {a} \sqrt {4 b-3 (-1)^{2/3} \sqrt [3]{a} c^{2/3}}}\right )}{3 \sqrt {3} \left (1+\sqrt [3]{-1}\right )^2 a^{5/6} b^2 \sqrt {4 b-3 (-1)^{2/3} \sqrt [3]{a} c^{2/3}} c^{2/3}}-\frac {\left (2 b-3 \sqrt [3]{a} c^{2/3}\right ) \arctan \left (\frac {3 a^{2/3} \sqrt [3]{c}+2 b x}{\sqrt {3} \sqrt {a} \sqrt {4 b-3 \sqrt [3]{a} c^{2/3}}}\right )}{9 \sqrt {3} a^{5/6} b^2 \sqrt {4 b-3 \sqrt [3]{a} c^{2/3}} c^{2/3}}-\frac {(-1)^{2/3} \left (2 b+3 \sqrt [3]{-1} \sqrt [3]{a} c^{2/3}\right ) \arctan \left (\frac {3 (-1)^{2/3} a^{2/3} \sqrt [3]{c}+2 b x}{\sqrt {3} \sqrt {a} \sqrt {4 b+3 \sqrt [3]{-1} \sqrt [3]{a} c^{2/3}}}\right )}{3 \sqrt {3} \left (1-\sqrt [3]{-1}\right ) \left (1+\sqrt [3]{-1}\right )^2 a^{5/6} b^2 \sqrt {4 b+3 \sqrt [3]{-1} \sqrt [3]{a} c^{2/3}} c^{2/3}}-\frac {\log \left (3 a+3 a^{2/3} \sqrt [3]{c} x+b x^2\right )}{18 a^{2/3} b^2 \sqrt [3]{c}}+\frac {\log \left (3 a-3 \sqrt [3]{-1} a^{2/3} \sqrt [3]{c} x+b x^2\right )}{6 \left (1+\sqrt [3]{-1}\right )^2 a^{2/3} b^2 \sqrt [3]{c}}+\frac {\sqrt [3]{-1} \log \left (3 a+3 (-1)^{2/3} a^{2/3} \sqrt [3]{c} x+b x^2\right )}{18 a^{2/3} b^2 \sqrt [3]{c}} \]

output
-1/18*ln(3*a+3*a^(2/3)*c^(1/3)*x+b*x^2)/a^(2/3)/b^2/c^(1/3)+1/6*ln(3*a-3*( 
-1)^(1/3)*a^(2/3)*c^(1/3)*x+b*x^2)/(1+(-1)^(1/3))^2/a^(2/3)/b^2/c^(1/3)+1/ 
18*(-1)^(1/3)*ln(3*a+3*(-1)^(2/3)*a^(2/3)*c^(1/3)*x+b*x^2)/a^(2/3)/b^2/c^( 
1/3)-1/27*(2*b-3*a^(1/3)*c^(2/3))*arctan(1/3*(3*a^(2/3)*c^(1/3)+2*b*x)*3^( 
1/2)/a^(1/2)/(4*b-3*a^(1/3)*c^(2/3))^(1/2))/a^(5/6)/b^2/c^(2/3)*3^(1/2)/(4 
*b-3*a^(1/3)*c^(2/3))^(1/2)-1/9*(-1)^(2/3)*(2*b+3*(-1)^(1/3)*a^(1/3)*c^(2/ 
3))*arctan(1/3*(3*(-1)^(2/3)*a^(2/3)*c^(1/3)+2*b*x)*3^(1/2)/a^(1/2)/(4*b+3 
*(-1)^(1/3)*a^(1/3)*c^(2/3))^(1/2))/(1-(-1)^(1/3))/(1+(-1)^(1/3))^2/a^(5/6 
)/b^2/c^(2/3)*3^(1/2)/(4*b+3*(-1)^(1/3)*a^(1/3)*c^(2/3))^(1/2)-1/9*(-1)^(1 
/3)*(2*(-1)^(1/3)*b+3*a^(1/3)*c^(2/3))*arctan(1/3*(3*(-1)^(1/3)*a^(2/3)*c^ 
(1/3)-2*b*x)*3^(1/2)/a^(1/2)/(4*b-3*(-1)^(2/3)*a^(1/3)*c^(2/3))^(1/2))/(1+ 
(-1)^(1/3))^2/a^(5/6)/b^2/c^(2/3)*3^(1/2)/(4*b-3*(-1)^(2/3)*a^(1/3)*c^(2/3 
))^(1/2)
 
3.2.36.2 Mathematica [C] (verified)

Result contains higher order function than in optimal. Order 9 vs. order 3 in optimal.

Time = 0.04 (sec) , antiderivative size = 99, normalized size of antiderivative = 0.18 \[ \int \frac {x^4}{27 a^3+27 a^2 b x^2+27 a^2 c x^3+9 a b^2 x^4+b^3 x^6} \, dx=\frac {1}{3} \text {RootSum}\left [27 a^3+27 a^2 b \text {$\#$1}^2+27 a^2 c \text {$\#$1}^3+9 a b^2 \text {$\#$1}^4+b^3 \text {$\#$1}^6\&,\frac {\log (x-\text {$\#$1}) \text {$\#$1}^3}{18 a^2 b+27 a^2 c \text {$\#$1}+12 a b^2 \text {$\#$1}^2+2 b^3 \text {$\#$1}^4}\&\right ] \]

input
Integrate[x^4/(27*a^3 + 27*a^2*b*x^2 + 27*a^2*c*x^3 + 9*a*b^2*x^4 + b^3*x^ 
6),x]
 
output
RootSum[27*a^3 + 27*a^2*b*#1^2 + 27*a^2*c*#1^3 + 9*a*b^2*#1^4 + b^3*#1^6 & 
 , (Log[x - #1]*#1^3)/(18*a^2*b + 27*a^2*c*#1 + 12*a*b^2*#1^2 + 2*b^3*#1^4 
) & ]/3
 
3.2.36.3 Rubi [A] (verified)

Time = 1.28 (sec) , antiderivative size = 530, normalized size of antiderivative = 0.97, number of steps used = 2, number of rules used = 2, \(\frac {\text {number of rules}}{\text {integrand size}}\) = 0.043, Rules used = {2466, 2009}

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^4}{27 a^3+27 a^2 b x^2+27 a^2 c x^3+9 a b^2 x^4+b^3 x^6} \, dx\)

\(\Big \downarrow \) 2466

\(\displaystyle 19683 a^6 \int \left (-\frac {\sqrt [3]{c} x+\sqrt [3]{a}}{177147 a^{20/3} b c^{2/3} \left (b x^2+3 a^{2/3} \sqrt [3]{c} x+3 a\right )}+\frac {\sqrt [3]{c} x+(-1)^{2/3} \sqrt [3]{a}}{59049 \left (1+\sqrt [3]{-1}\right )^2 a^{20/3} b c^{2/3} \left (b x^2-3 \sqrt [3]{-1} a^{2/3} \sqrt [3]{c} x+3 a\right )}-\frac {(-1)^{2/3} \left ((-1)^{2/3} \sqrt [3]{c} x+\sqrt [3]{a}\right )}{177147 a^{20/3} b c^{2/3} \left (b x^2+3 (-1)^{2/3} a^{2/3} \sqrt [3]{c} x+3 a\right )}\right )dx\)

\(\Big \downarrow \) 2009

\(\displaystyle 19683 a^6 \left (-\frac {\sqrt [3]{-1} \left (3 \sqrt [3]{a} c^{2/3}+2 \sqrt [3]{-1} b\right ) \arctan \left (\frac {3 \sqrt [3]{-1} a^{2/3} \sqrt [3]{c}-2 b x}{\sqrt {3} \sqrt {a} \sqrt {4 b-3 (-1)^{2/3} \sqrt [3]{a} c^{2/3}}}\right )}{59049 \sqrt {3} \left (1+\sqrt [3]{-1}\right )^2 a^{41/6} b^2 c^{2/3} \sqrt {4 b-3 (-1)^{2/3} \sqrt [3]{a} c^{2/3}}}-\frac {\left (2 b-3 \sqrt [3]{a} c^{2/3}\right ) \arctan \left (\frac {3 a^{2/3} \sqrt [3]{c}+2 b x}{\sqrt {3} \sqrt {a} \sqrt {4 b-3 \sqrt [3]{a} c^{2/3}}}\right )}{177147 \sqrt {3} a^{41/6} b^2 c^{2/3} \sqrt {4 b-3 \sqrt [3]{a} c^{2/3}}}-\frac {(-1)^{2/3} \left (3 \sqrt [3]{-1} \sqrt [3]{a} c^{2/3}+2 b\right ) \arctan \left (\frac {3 (-1)^{2/3} a^{2/3} \sqrt [3]{c}+2 b x}{\sqrt {3} \sqrt {a} \sqrt {3 \sqrt [3]{-1} \sqrt [3]{a} c^{2/3}+4 b}}\right )}{177147 \sqrt {3} a^{41/6} b^2 c^{2/3} \sqrt {3 \sqrt [3]{-1} \sqrt [3]{a} c^{2/3}+4 b}}-\frac {\log \left (3 a^{2/3} \sqrt [3]{c} x+3 a+b x^2\right )}{354294 a^{20/3} b^2 \sqrt [3]{c}}+\frac {\log \left (-3 \sqrt [3]{-1} a^{2/3} \sqrt [3]{c} x+3 a+b x^2\right )}{118098 \left (1+\sqrt [3]{-1}\right )^2 a^{20/3} b^2 \sqrt [3]{c}}+\frac {\sqrt [3]{-1} \log \left (3 (-1)^{2/3} a^{2/3} \sqrt [3]{c} x+3 a+b x^2\right )}{354294 a^{20/3} b^2 \sqrt [3]{c}}\right )\)

input
Int[x^4/(27*a^3 + 27*a^2*b*x^2 + 27*a^2*c*x^3 + 9*a*b^2*x^4 + b^3*x^6),x]
 
output
19683*a^6*(-1/59049*((-1)^(1/3)*(2*(-1)^(1/3)*b + 3*a^(1/3)*c^(2/3))*ArcTa 
n[(3*(-1)^(1/3)*a^(2/3)*c^(1/3) - 2*b*x)/(Sqrt[3]*Sqrt[a]*Sqrt[4*b - 3*(-1 
)^(2/3)*a^(1/3)*c^(2/3)])])/(Sqrt[3]*(1 + (-1)^(1/3))^2*a^(41/6)*b^2*Sqrt[ 
4*b - 3*(-1)^(2/3)*a^(1/3)*c^(2/3)]*c^(2/3)) - ((2*b - 3*a^(1/3)*c^(2/3))* 
ArcTan[(3*a^(2/3)*c^(1/3) + 2*b*x)/(Sqrt[3]*Sqrt[a]*Sqrt[4*b - 3*a^(1/3)*c 
^(2/3)])])/(177147*Sqrt[3]*a^(41/6)*b^2*Sqrt[4*b - 3*a^(1/3)*c^(2/3)]*c^(2 
/3)) - ((-1)^(2/3)*(2*b + 3*(-1)^(1/3)*a^(1/3)*c^(2/3))*ArcTan[(3*(-1)^(2/ 
3)*a^(2/3)*c^(1/3) + 2*b*x)/(Sqrt[3]*Sqrt[a]*Sqrt[4*b + 3*(-1)^(1/3)*a^(1/ 
3)*c^(2/3)])])/(177147*Sqrt[3]*a^(41/6)*b^2*Sqrt[4*b + 3*(-1)^(1/3)*a^(1/3 
)*c^(2/3)]*c^(2/3)) - Log[3*a + 3*a^(2/3)*c^(1/3)*x + b*x^2]/(354294*a^(20 
/3)*b^2*c^(1/3)) + Log[3*a - 3*(-1)^(1/3)*a^(2/3)*c^(1/3)*x + b*x^2]/(1180 
98*(1 + (-1)^(1/3))^2*a^(20/3)*b^2*c^(1/3)) + ((-1)^(1/3)*Log[3*a + 3*(-1) 
^(2/3)*a^(2/3)*c^(1/3)*x + b*x^2])/(354294*a^(20/3)*b^2*c^(1/3)))
 

3.2.36.3.1 Defintions of rubi rules used

rule 2009
Int[u_, x_Symbol] :> Simp[IntSum[u, x], x] /; SumQ[u]
 

rule 2466
Int[(u_.)*(Q6_)^(p_), x_Symbol] :> With[{a = Coeff[Q6, x, 0], b = Coeff[Q6, 
 x, 2], c = Coeff[Q6, x, 3], d = Coeff[Q6, x, 4], e = Coeff[Q6, x, 6]}, Sim 
p[1/(3^(3*p)*a^(2*p))   Int[ExpandIntegrand[u*(3*a + 3*Rt[a, 3]^2*Rt[c, 3]* 
x + b*x^2)^p*(3*a - 3*(-1)^(1/3)*Rt[a, 3]^2*Rt[c, 3]*x + b*x^2)^p*(3*a + 3* 
(-1)^(2/3)*Rt[a, 3]^2*Rt[c, 3]*x + b*x^2)^p, x], x], x] /; EqQ[b^2 - 3*a*d, 
 0] && EqQ[b^3 - 27*a^2*e, 0]] /; ILtQ[p, 0] && PolyQ[Q6, x, 6] && EqQ[Coef 
f[Q6, x, 1], 0] && EqQ[Coeff[Q6, x, 5], 0] && RationalFunctionQ[u, x]
 
3.2.36.4 Maple [C] (verified)

Result contains higher order function than in optimal. Order 9 vs. order 3.

Time = 0.08 (sec) , antiderivative size = 93, normalized size of antiderivative = 0.17

method result size
default \(\frac {\left (\munderset {\textit {\_R} =\operatorname {RootOf}\left (b^{3} \textit {\_Z}^{6}+9 b^{2} a \,\textit {\_Z}^{4}+27 c \,a^{2} \textit {\_Z}^{3}+27 a^{2} b \,\textit {\_Z}^{2}+27 a^{3}\right )}{\sum }\frac {\textit {\_R}^{4} \ln \left (x -\textit {\_R} \right )}{2 \textit {\_R}^{5} b^{3}+12 \textit {\_R}^{3} a \,b^{2}+27 \textit {\_R}^{2} a^{2} c +18 a^{2} b \textit {\_R}}\right )}{3}\) \(93\)
risch \(\frac {\left (\munderset {\textit {\_R} =\operatorname {RootOf}\left (b^{3} \textit {\_Z}^{6}+9 b^{2} a \,\textit {\_Z}^{4}+27 c \,a^{2} \textit {\_Z}^{3}+27 a^{2} b \,\textit {\_Z}^{2}+27 a^{3}\right )}{\sum }\frac {\textit {\_R}^{4} \ln \left (x -\textit {\_R} \right )}{2 \textit {\_R}^{5} b^{3}+12 \textit {\_R}^{3} a \,b^{2}+27 \textit {\_R}^{2} a^{2} c +18 a^{2} b \textit {\_R}}\right )}{3}\) \(93\)

input
int(x^4/(b^3*x^6+9*a*b^2*x^4+27*a^2*c*x^3+27*a^2*b*x^2+27*a^3),x,method=_R 
ETURNVERBOSE)
 
output
1/3*sum(_R^4/(2*_R^5*b^3+12*_R^3*a*b^2+27*_R^2*a^2*c+18*_R*a^2*b)*ln(x-_R) 
,_R=RootOf(_Z^6*b^3+9*_Z^4*a*b^2+27*_Z^3*a^2*c+27*_Z^2*a^2*b+27*a^3))
 
3.2.36.5 Fricas [F(-2)]

Exception generated. \[ \int \frac {x^4}{27 a^3+27 a^2 b x^2+27 a^2 c x^3+9 a b^2 x^4+b^3 x^6} \, dx=\text {Exception raised: RuntimeError} \]

input
integrate(x^4/(b^3*x^6+9*a*b^2*x^4+27*a^2*c*x^3+27*a^2*b*x^2+27*a^3),x, al 
gorithm="fricas")
 
output
Exception raised: RuntimeError >> no explicit roots found
 
3.2.36.6 Sympy [F(-1)]

Timed out. \[ \int \frac {x^4}{27 a^3+27 a^2 b x^2+27 a^2 c x^3+9 a b^2 x^4+b^3 x^6} \, dx=\text {Timed out} \]

input
integrate(x**4/(b**3*x**6+9*a*b**2*x**4+27*a**2*c*x**3+27*a**2*b*x**2+27*a 
**3),x)
 
output
Timed out
 
3.2.36.7 Maxima [F]

\[ \int \frac {x^4}{27 a^3+27 a^2 b x^2+27 a^2 c x^3+9 a b^2 x^4+b^3 x^6} \, dx=\int { \frac {x^{4}}{b^{3} x^{6} + 9 \, a b^{2} x^{4} + 27 \, a^{2} c x^{3} + 27 \, a^{2} b x^{2} + 27 \, a^{3}} \,d x } \]

input
integrate(x^4/(b^3*x^6+9*a*b^2*x^4+27*a^2*c*x^3+27*a^2*b*x^2+27*a^3),x, al 
gorithm="maxima")
 
output
integrate(x^4/(b^3*x^6 + 9*a*b^2*x^4 + 27*a^2*c*x^3 + 27*a^2*b*x^2 + 27*a^ 
3), x)
 
3.2.36.8 Giac [F]

\[ \int \frac {x^4}{27 a^3+27 a^2 b x^2+27 a^2 c x^3+9 a b^2 x^4+b^3 x^6} \, dx=\int { \frac {x^{4}}{b^{3} x^{6} + 9 \, a b^{2} x^{4} + 27 \, a^{2} c x^{3} + 27 \, a^{2} b x^{2} + 27 \, a^{3}} \,d x } \]

input
integrate(x^4/(b^3*x^6+9*a*b^2*x^4+27*a^2*c*x^3+27*a^2*b*x^2+27*a^3),x, al 
gorithm="giac")
 
output
integrate(x^4/(b^3*x^6 + 9*a*b^2*x^4 + 27*a^2*c*x^3 + 27*a^2*b*x^2 + 27*a^ 
3), x)
 
3.2.36.9 Mupad [B] (verification not implemented)

Time = 9.52 (sec) , antiderivative size = 1563, normalized size of antiderivative = 2.87 \[ \int \frac {x^4}{27 a^3+27 a^2 b x^2+27 a^2 c x^3+9 a b^2 x^4+b^3 x^6} \, dx=\text {Too large to display} \]

input
int(x^4/(27*a^3 + b^3*x^6 + 27*a^2*b*x^2 + 9*a*b^2*x^4 + 27*a^2*c*x^3),x)
 
output
symsum(log(-19683*a^8*b^3*(c*x - b + 6561*root(918330048*a^5*b^9*c^4*z^6 - 
 387420489*a^6*b^6*c^6*z^6 + 1594323*a^4*b^4*c^4*z^4 + 1023516*a^3*b^3*c^3 
*z^3 - 531441*a^4*c^5*z^3 + 32805*a^2*b^2*c^2*z^2 + 324*a*b*c*z + 1, z, k) 
^2*a^3*c^4 + 2*root(918330048*a^5*b^9*c^4*z^6 - 387420489*a^6*b^6*c^6*z^6 
+ 1594323*a^4*b^4*c^4*z^4 + 1023516*a^3*b^3*c^3*z^3 - 531441*a^4*c^5*z^3 + 
 32805*a^2*b^2*c^2*z^2 + 324*a*b*c*z + 1, z, k)*b^4*x - 198*root(918330048 
*a^5*b^9*c^4*z^6 - 387420489*a^6*b^6*c^6*z^6 + 1594323*a^4*b^4*c^4*z^4 + 1 
023516*a^3*b^3*c^3*z^3 - 531441*a^4*c^5*z^3 + 32805*a^2*b^2*c^2*z^2 + 324* 
a*b*c*z + 1, z, k)*a*b^2*c - 8991*root(918330048*a^5*b^9*c^4*z^6 - 3874204 
89*a^6*b^6*c^6*z^6 + 1594323*a^4*b^4*c^4*z^4 + 1023516*a^3*b^3*c^3*z^3 - 5 
31441*a^4*c^5*z^3 + 32805*a^2*b^2*c^2*z^2 + 324*a*b*c*z + 1, z, k)^2*a^2*b 
^3*c^2 - 19683*root(918330048*a^5*b^9*c^4*z^6 - 387420489*a^6*b^6*c^6*z^6 
+ 1594323*a^4*b^4*c^4*z^4 + 1023516*a^3*b^3*c^3*z^3 - 531441*a^4*c^5*z^3 + 
 32805*a^2*b^2*c^2*z^2 + 324*a*b*c*z + 1, z, k)^3*a^3*b^4*c^3 + 104976*roo 
t(918330048*a^5*b^9*c^4*z^6 - 387420489*a^6*b^6*c^6*z^6 + 1594323*a^4*b^4* 
c^4*z^4 + 1023516*a^3*b^3*c^3*z^3 - 531441*a^4*c^5*z^3 + 32805*a^2*b^2*c^2 
*z^2 + 324*a*b*c*z + 1, z, k)^4*a^3*b^8*c^2 - 8503056*root(918330048*a^5*b 
^9*c^4*z^6 - 387420489*a^6*b^6*c^6*z^6 + 1594323*a^4*b^4*c^4*z^4 + 1023516 
*a^3*b^3*c^3*z^3 - 531441*a^4*c^5*z^3 + 32805*a^2*b^2*c^2*z^2 + 324*a*b*c* 
z + 1, z, k)^5*a^4*b^9*c^3 + 4782969*root(918330048*a^5*b^9*c^4*z^6 - 3...