Integrand size = 13, antiderivative size = 61 \[ \int x^m \left (a+b x^3\right )^3 \, dx=\frac {a^3 x^{1+m}}{1+m}+\frac {3 a^2 b x^{4+m}}{4+m}+\frac {3 a b^2 x^{7+m}}{7+m}+\frac {b^3 x^{10+m}}{10+m} \] Output:
a^3*x^(1+m)/(1+m)+3*a^2*b*x^(4+m)/(4+m)+3*a*b^2*x^(7+m)/(7+m)+b^3*x^(10+m) /(10+m)
Time = 0.06 (sec) , antiderivative size = 56, normalized size of antiderivative = 0.92 \[ \int x^m \left (a+b x^3\right )^3 \, dx=x^{1+m} \left (\frac {a^3}{1+m}+\frac {3 a^2 b x^3}{4+m}+\frac {3 a b^2 x^6}{7+m}+\frac {b^3 x^9}{10+m}\right ) \] Input:
Integrate[x^m*(a + b*x^3)^3,x]
Output:
x^(1 + m)*(a^3/(1 + m) + (3*a^2*b*x^3)/(4 + m) + (3*a*b^2*x^6)/(7 + m) + ( b^3*x^9)/(10 + m))
Time = 0.32 (sec) , antiderivative size = 61, normalized size of antiderivative = 1.00, number of steps used = 2, number of rules used = 2, \(\frac {\text {number of rules}}{\text {integrand size}}\) = 0.154, Rules used = {802, 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 x^m \left (a+b x^3\right )^3 \, dx\) |
\(\Big \downarrow \) 802 |
\(\displaystyle \int \left (a^3 x^m+3 a^2 b x^{m+3}+3 a b^2 x^{m+6}+b^3 x^{m+9}\right )dx\) |
\(\Big \downarrow \) 2009 |
\(\displaystyle \frac {a^3 x^{m+1}}{m+1}+\frac {3 a^2 b x^{m+4}}{m+4}+\frac {3 a b^2 x^{m+7}}{m+7}+\frac {b^3 x^{m+10}}{m+10}\) |
Input:
Int[x^m*(a + b*x^3)^3,x]
Output:
(a^3*x^(1 + m))/(1 + m) + (3*a^2*b*x^(4 + m))/(4 + m) + (3*a*b^2*x^(7 + m) )/(7 + m) + (b^3*x^(10 + m))/(10 + m)
Int[((c_.)*(x_))^(m_.)*((a_) + (b_.)*(x_)^(n_))^(p_.), x_Symbol] :> Int[Exp andIntegrand[(c*x)^m*(a + b*x^n)^p, x], x] /; FreeQ[{a, b, c, m, n}, x] && IGtQ[p, 0]
Leaf count of result is larger than twice the leaf count of optimal. \(176\) vs. \(2(61)=122\).
Time = 0.56 (sec) , antiderivative size = 177, normalized size of antiderivative = 2.90
method | result | size |
risch | \(\frac {x \left (b^{3} m^{3} x^{9}+12 b^{3} m^{2} x^{9}+39 m \,x^{9} b^{3}+3 a \,b^{2} m^{3} x^{6}+28 b^{3} x^{9}+45 a \,b^{2} m^{2} x^{6}+162 m \,x^{6} a \,b^{2}+3 a^{2} b \,m^{3} x^{3}+120 a \,b^{2} x^{6}+54 a^{2} b \,m^{2} x^{3}+261 m \,x^{3} a^{2} b +a^{3} m^{3}+210 a^{2} b \,x^{3}+21 a^{3} m^{2}+138 m \,a^{3}+280 a^{3}\right ) x^{m}}{\left (10+m \right ) \left (7+m \right ) \left (4+m \right ) \left (1+m \right )}\) | \(177\) |
orering | \(\frac {x \left (b^{3} m^{3} x^{9}+12 b^{3} m^{2} x^{9}+39 m \,x^{9} b^{3}+3 a \,b^{2} m^{3} x^{6}+28 b^{3} x^{9}+45 a \,b^{2} m^{2} x^{6}+162 m \,x^{6} a \,b^{2}+3 a^{2} b \,m^{3} x^{3}+120 a \,b^{2} x^{6}+54 a^{2} b \,m^{2} x^{3}+261 m \,x^{3} a^{2} b +a^{3} m^{3}+210 a^{2} b \,x^{3}+21 a^{3} m^{2}+138 m \,a^{3}+280 a^{3}\right ) x^{m}}{\left (10+m \right ) \left (7+m \right ) \left (4+m \right ) \left (1+m \right )}\) | \(177\) |
gosper | \(\frac {x^{1+m} \left (b^{3} m^{3} x^{9}+12 b^{3} m^{2} x^{9}+39 m \,x^{9} b^{3}+3 a \,b^{2} m^{3} x^{6}+28 b^{3} x^{9}+45 a \,b^{2} m^{2} x^{6}+162 m \,x^{6} a \,b^{2}+3 a^{2} b \,m^{3} x^{3}+120 a \,b^{2} x^{6}+54 a^{2} b \,m^{2} x^{3}+261 m \,x^{3} a^{2} b +a^{3} m^{3}+210 a^{2} b \,x^{3}+21 a^{3} m^{2}+138 m \,a^{3}+280 a^{3}\right )}{\left (1+m \right ) \left (4+m \right ) \left (7+m \right ) \left (10+m \right )}\) | \(178\) |
parallelrisch | \(\frac {x^{10} x^{m} b^{3} m^{3}+12 x^{10} x^{m} b^{3} m^{2}+39 x^{10} x^{m} b^{3} m +28 x^{10} x^{m} b^{3}+3 x^{7} x^{m} a \,b^{2} m^{3}+45 x^{7} x^{m} a \,b^{2} m^{2}+162 x^{7} x^{m} a \,b^{2} m +120 x^{7} x^{m} a \,b^{2}+3 x^{4} x^{m} a^{2} b \,m^{3}+54 x^{4} x^{m} a^{2} b \,m^{2}+261 x^{4} x^{m} a^{2} b m +210 x^{4} x^{m} a^{2} b +x \,x^{m} a^{3} m^{3}+21 x \,x^{m} a^{3} m^{2}+138 x \,x^{m} a^{3} m +280 x \,x^{m} a^{3}}{\left (10+m \right ) \left (7+m \right ) \left (4+m \right ) \left (1+m \right )}\) | \(225\) |
Input:
int(x^m*(b*x^3+a)^3,x,method=_RETURNVERBOSE)
Output:
x*(b^3*m^3*x^9+12*b^3*m^2*x^9+39*b^3*m*x^9+3*a*b^2*m^3*x^6+28*b^3*x^9+45*a *b^2*m^2*x^6+162*a*b^2*m*x^6+3*a^2*b*m^3*x^3+120*a*b^2*x^6+54*a^2*b*m^2*x^ 3+261*a^2*b*m*x^3+a^3*m^3+210*a^2*b*x^3+21*a^3*m^2+138*a^3*m+280*a^3)*x^m/ (10+m)/(7+m)/(4+m)/(1+m)
Leaf count of result is larger than twice the leaf count of optimal. 157 vs. \(2 (61) = 122\).
Time = 0.08 (sec) , antiderivative size = 157, normalized size of antiderivative = 2.57 \[ \int x^m \left (a+b x^3\right )^3 \, dx=\frac {{\left ({\left (b^{3} m^{3} + 12 \, b^{3} m^{2} + 39 \, b^{3} m + 28 \, b^{3}\right )} x^{10} + 3 \, {\left (a b^{2} m^{3} + 15 \, a b^{2} m^{2} + 54 \, a b^{2} m + 40 \, a b^{2}\right )} x^{7} + 3 \, {\left (a^{2} b m^{3} + 18 \, a^{2} b m^{2} + 87 \, a^{2} b m + 70 \, a^{2} b\right )} x^{4} + {\left (a^{3} m^{3} + 21 \, a^{3} m^{2} + 138 \, a^{3} m + 280 \, a^{3}\right )} x\right )} x^{m}}{m^{4} + 22 \, m^{3} + 159 \, m^{2} + 418 \, m + 280} \] Input:
integrate(x^m*(b*x^3+a)^3,x, algorithm="fricas")
Output:
((b^3*m^3 + 12*b^3*m^2 + 39*b^3*m + 28*b^3)*x^10 + 3*(a*b^2*m^3 + 15*a*b^2 *m^2 + 54*a*b^2*m + 40*a*b^2)*x^7 + 3*(a^2*b*m^3 + 18*a^2*b*m^2 + 87*a^2*b *m + 70*a^2*b)*x^4 + (a^3*m^3 + 21*a^3*m^2 + 138*a^3*m + 280*a^3)*x)*x^m/( m^4 + 22*m^3 + 159*m^2 + 418*m + 280)
Leaf count of result is larger than twice the leaf count of optimal. 666 vs. \(2 (53) = 106\).
Time = 0.57 (sec) , antiderivative size = 666, normalized size of antiderivative = 10.92 \[ \int x^m \left (a+b x^3\right )^3 \, dx=\begin {cases} - \frac {a^{3}}{9 x^{9}} - \frac {a^{2} b}{2 x^{6}} - \frac {a b^{2}}{x^{3}} + b^{3} \log {\left (x \right )} & \text {for}\: m = -10 \\- \frac {a^{3}}{6 x^{6}} - \frac {a^{2} b}{x^{3}} + 3 a b^{2} \log {\left (x \right )} + \frac {b^{3} x^{3}}{3} & \text {for}\: m = -7 \\- \frac {a^{3}}{3 x^{3}} + 3 a^{2} b \log {\left (x \right )} + a b^{2} x^{3} + \frac {b^{3} x^{6}}{6} & \text {for}\: m = -4 \\a^{3} \log {\left (x \right )} + a^{2} b x^{3} + \frac {a b^{2} x^{6}}{2} + \frac {b^{3} x^{9}}{9} & \text {for}\: m = -1 \\\frac {a^{3} m^{3} x x^{m}}{m^{4} + 22 m^{3} + 159 m^{2} + 418 m + 280} + \frac {21 a^{3} m^{2} x x^{m}}{m^{4} + 22 m^{3} + 159 m^{2} + 418 m + 280} + \frac {138 a^{3} m x x^{m}}{m^{4} + 22 m^{3} + 159 m^{2} + 418 m + 280} + \frac {280 a^{3} x x^{m}}{m^{4} + 22 m^{3} + 159 m^{2} + 418 m + 280} + \frac {3 a^{2} b m^{3} x^{4} x^{m}}{m^{4} + 22 m^{3} + 159 m^{2} + 418 m + 280} + \frac {54 a^{2} b m^{2} x^{4} x^{m}}{m^{4} + 22 m^{3} + 159 m^{2} + 418 m + 280} + \frac {261 a^{2} b m x^{4} x^{m}}{m^{4} + 22 m^{3} + 159 m^{2} + 418 m + 280} + \frac {210 a^{2} b x^{4} x^{m}}{m^{4} + 22 m^{3} + 159 m^{2} + 418 m + 280} + \frac {3 a b^{2} m^{3} x^{7} x^{m}}{m^{4} + 22 m^{3} + 159 m^{2} + 418 m + 280} + \frac {45 a b^{2} m^{2} x^{7} x^{m}}{m^{4} + 22 m^{3} + 159 m^{2} + 418 m + 280} + \frac {162 a b^{2} m x^{7} x^{m}}{m^{4} + 22 m^{3} + 159 m^{2} + 418 m + 280} + \frac {120 a b^{2} x^{7} x^{m}}{m^{4} + 22 m^{3} + 159 m^{2} + 418 m + 280} + \frac {b^{3} m^{3} x^{10} x^{m}}{m^{4} + 22 m^{3} + 159 m^{2} + 418 m + 280} + \frac {12 b^{3} m^{2} x^{10} x^{m}}{m^{4} + 22 m^{3} + 159 m^{2} + 418 m + 280} + \frac {39 b^{3} m x^{10} x^{m}}{m^{4} + 22 m^{3} + 159 m^{2} + 418 m + 280} + \frac {28 b^{3} x^{10} x^{m}}{m^{4} + 22 m^{3} + 159 m^{2} + 418 m + 280} & \text {otherwise} \end {cases} \] Input:
integrate(x**m*(b*x**3+a)**3,x)
Output:
Piecewise((-a**3/(9*x**9) - a**2*b/(2*x**6) - a*b**2/x**3 + b**3*log(x), E q(m, -10)), (-a**3/(6*x**6) - a**2*b/x**3 + 3*a*b**2*log(x) + b**3*x**3/3, Eq(m, -7)), (-a**3/(3*x**3) + 3*a**2*b*log(x) + a*b**2*x**3 + b**3*x**6/6 , Eq(m, -4)), (a**3*log(x) + a**2*b*x**3 + a*b**2*x**6/2 + b**3*x**9/9, Eq (m, -1)), (a**3*m**3*x*x**m/(m**4 + 22*m**3 + 159*m**2 + 418*m + 280) + 21 *a**3*m**2*x*x**m/(m**4 + 22*m**3 + 159*m**2 + 418*m + 280) + 138*a**3*m*x *x**m/(m**4 + 22*m**3 + 159*m**2 + 418*m + 280) + 280*a**3*x*x**m/(m**4 + 22*m**3 + 159*m**2 + 418*m + 280) + 3*a**2*b*m**3*x**4*x**m/(m**4 + 22*m** 3 + 159*m**2 + 418*m + 280) + 54*a**2*b*m**2*x**4*x**m/(m**4 + 22*m**3 + 1 59*m**2 + 418*m + 280) + 261*a**2*b*m*x**4*x**m/(m**4 + 22*m**3 + 159*m**2 + 418*m + 280) + 210*a**2*b*x**4*x**m/(m**4 + 22*m**3 + 159*m**2 + 418*m + 280) + 3*a*b**2*m**3*x**7*x**m/(m**4 + 22*m**3 + 159*m**2 + 418*m + 280) + 45*a*b**2*m**2*x**7*x**m/(m**4 + 22*m**3 + 159*m**2 + 418*m + 280) + 16 2*a*b**2*m*x**7*x**m/(m**4 + 22*m**3 + 159*m**2 + 418*m + 280) + 120*a*b** 2*x**7*x**m/(m**4 + 22*m**3 + 159*m**2 + 418*m + 280) + b**3*m**3*x**10*x* *m/(m**4 + 22*m**3 + 159*m**2 + 418*m + 280) + 12*b**3*m**2*x**10*x**m/(m* *4 + 22*m**3 + 159*m**2 + 418*m + 280) + 39*b**3*m*x**10*x**m/(m**4 + 22*m **3 + 159*m**2 + 418*m + 280) + 28*b**3*x**10*x**m/(m**4 + 22*m**3 + 159*m **2 + 418*m + 280), True))
Time = 0.03 (sec) , antiderivative size = 61, normalized size of antiderivative = 1.00 \[ \int x^m \left (a+b x^3\right )^3 \, dx=\frac {b^{3} x^{m + 10}}{m + 10} + \frac {3 \, a b^{2} x^{m + 7}}{m + 7} + \frac {3 \, a^{2} b x^{m + 4}}{m + 4} + \frac {a^{3} x^{m + 1}}{m + 1} \] Input:
integrate(x^m*(b*x^3+a)^3,x, algorithm="maxima")
Output:
b^3*x^(m + 10)/(m + 10) + 3*a*b^2*x^(m + 7)/(m + 7) + 3*a^2*b*x^(m + 4)/(m + 4) + a^3*x^(m + 1)/(m + 1)
Leaf count of result is larger than twice the leaf count of optimal. 224 vs. \(2 (61) = 122\).
Time = 0.13 (sec) , antiderivative size = 224, normalized size of antiderivative = 3.67 \[ \int x^m \left (a+b x^3\right )^3 \, dx=\frac {b^{3} m^{3} x^{10} x^{m} + 12 \, b^{3} m^{2} x^{10} x^{m} + 39 \, b^{3} m x^{10} x^{m} + 3 \, a b^{2} m^{3} x^{7} x^{m} + 28 \, b^{3} x^{10} x^{m} + 45 \, a b^{2} m^{2} x^{7} x^{m} + 162 \, a b^{2} m x^{7} x^{m} + 3 \, a^{2} b m^{3} x^{4} x^{m} + 120 \, a b^{2} x^{7} x^{m} + 54 \, a^{2} b m^{2} x^{4} x^{m} + 261 \, a^{2} b m x^{4} x^{m} + a^{3} m^{3} x x^{m} + 210 \, a^{2} b x^{4} x^{m} + 21 \, a^{3} m^{2} x x^{m} + 138 \, a^{3} m x x^{m} + 280 \, a^{3} x x^{m}}{m^{4} + 22 \, m^{3} + 159 \, m^{2} + 418 \, m + 280} \] Input:
integrate(x^m*(b*x^3+a)^3,x, algorithm="giac")
Output:
(b^3*m^3*x^10*x^m + 12*b^3*m^2*x^10*x^m + 39*b^3*m*x^10*x^m + 3*a*b^2*m^3* x^7*x^m + 28*b^3*x^10*x^m + 45*a*b^2*m^2*x^7*x^m + 162*a*b^2*m*x^7*x^m + 3 *a^2*b*m^3*x^4*x^m + 120*a*b^2*x^7*x^m + 54*a^2*b*m^2*x^4*x^m + 261*a^2*b* m*x^4*x^m + a^3*m^3*x*x^m + 210*a^2*b*x^4*x^m + 21*a^3*m^2*x*x^m + 138*a^3 *m*x*x^m + 280*a^3*x*x^m)/(m^4 + 22*m^3 + 159*m^2 + 418*m + 280)
Time = 0.55 (sec) , antiderivative size = 167, normalized size of antiderivative = 2.74 \[ \int x^m \left (a+b x^3\right )^3 \, dx=x^m\,\left (\frac {a^3\,x\,\left (m^3+21\,m^2+138\,m+280\right )}{m^4+22\,m^3+159\,m^2+418\,m+280}+\frac {b^3\,x^{10}\,\left (m^3+12\,m^2+39\,m+28\right )}{m^4+22\,m^3+159\,m^2+418\,m+280}+\frac {3\,a\,b^2\,x^7\,\left (m^3+15\,m^2+54\,m+40\right )}{m^4+22\,m^3+159\,m^2+418\,m+280}+\frac {3\,a^2\,b\,x^4\,\left (m^3+18\,m^2+87\,m+70\right )}{m^4+22\,m^3+159\,m^2+418\,m+280}\right ) \] Input:
int(x^m*(a + b*x^3)^3,x)
Output:
x^m*((a^3*x*(138*m + 21*m^2 + m^3 + 280))/(418*m + 159*m^2 + 22*m^3 + m^4 + 280) + (b^3*x^10*(39*m + 12*m^2 + m^3 + 28))/(418*m + 159*m^2 + 22*m^3 + m^4 + 280) + (3*a*b^2*x^7*(54*m + 15*m^2 + m^3 + 40))/(418*m + 159*m^2 + 22*m^3 + m^4 + 280) + (3*a^2*b*x^4*(87*m + 18*m^2 + m^3 + 70))/(418*m + 15 9*m^2 + 22*m^3 + m^4 + 280))
Time = 0.21 (sec) , antiderivative size = 176, normalized size of antiderivative = 2.89 \[ \int x^m \left (a+b x^3\right )^3 \, dx=\frac {x^{m} x \left (b^{3} m^{3} x^{9}+12 b^{3} m^{2} x^{9}+39 b^{3} m \,x^{9}+3 a \,b^{2} m^{3} x^{6}+28 b^{3} x^{9}+45 a \,b^{2} m^{2} x^{6}+162 a \,b^{2} m \,x^{6}+3 a^{2} b \,m^{3} x^{3}+120 a \,b^{2} x^{6}+54 a^{2} b \,m^{2} x^{3}+261 a^{2} b m \,x^{3}+a^{3} m^{3}+210 a^{2} b \,x^{3}+21 a^{3} m^{2}+138 a^{3} m +280 a^{3}\right )}{m^{4}+22 m^{3}+159 m^{2}+418 m +280} \] Input:
int(x^m*(b*x^3+a)^3,x)
Output:
(x**m*x*(a**3*m**3 + 21*a**3*m**2 + 138*a**3*m + 280*a**3 + 3*a**2*b*m**3* x**3 + 54*a**2*b*m**2*x**3 + 261*a**2*b*m*x**3 + 210*a**2*b*x**3 + 3*a*b** 2*m**3*x**6 + 45*a*b**2*m**2*x**6 + 162*a*b**2*m*x**6 + 120*a*b**2*x**6 + b**3*m**3*x**9 + 12*b**3*m**2*x**9 + 39*b**3*m*x**9 + 28*b**3*x**9))/(m**4 + 22*m**3 + 159*m**2 + 418*m + 280)