3.105 \(\int \frac{x}{a+b (c+d x)^3} \, dx\)

Optimal. Leaf size=180 \[ -\frac{\left (\sqrt [3]{a}+\sqrt [3]{b} c\right ) \log \left (\sqrt [3]{a}+\sqrt [3]{b} (c+d x)\right )}{3 a^{2/3} b^{2/3} d^2}+\frac{\left (\sqrt [3]{a}+\sqrt [3]{b} c\right ) \log \left (a^{2/3}-\sqrt [3]{a} \sqrt [3]{b} (c+d x)+b^{2/3} (c+d x)^2\right )}{6 a^{2/3} b^{2/3} d^2}-\frac{\left (\sqrt [3]{a}-\sqrt [3]{b} c\right ) \tan ^{-1}\left (\frac{\sqrt [3]{a}-2 \sqrt [3]{b} (c+d x)}{\sqrt{3} \sqrt [3]{a}}\right )}{\sqrt{3} a^{2/3} b^{2/3} d^2} \]

[Out]

-(((a^(1/3) - b^(1/3)*c)*ArcTan[(a^(1/3) - 2*b^(1/3)*(c + d*x))/(Sqrt[3]*a^(1/3))])/(Sqrt[3]*a^(2/3)*b^(2/3)*d
^2)) - ((a^(1/3) + b^(1/3)*c)*Log[a^(1/3) + b^(1/3)*(c + d*x)])/(3*a^(2/3)*b^(2/3)*d^2) + ((a^(1/3) + b^(1/3)*
c)*Log[a^(2/3) - a^(1/3)*b^(1/3)*(c + d*x) + b^(2/3)*(c + d*x)^2])/(6*a^(2/3)*b^(2/3)*d^2)

________________________________________________________________________________________

Rubi [A]  time = 0.156257, antiderivative size = 180, normalized size of antiderivative = 1., number of steps used = 7, number of rules used = 7, integrand size = 15, \(\frac{\text{number of rules}}{\text{integrand size}}\) = 0.467, Rules used = {371, 1860, 31, 634, 617, 204, 628} \[ -\frac{\left (\sqrt [3]{a}+\sqrt [3]{b} c\right ) \log \left (\sqrt [3]{a}+\sqrt [3]{b} (c+d x)\right )}{3 a^{2/3} b^{2/3} d^2}+\frac{\left (\sqrt [3]{a}+\sqrt [3]{b} c\right ) \log \left (a^{2/3}-\sqrt [3]{a} \sqrt [3]{b} (c+d x)+b^{2/3} (c+d x)^2\right )}{6 a^{2/3} b^{2/3} d^2}-\frac{\left (\sqrt [3]{a}-\sqrt [3]{b} c\right ) \tan ^{-1}\left (\frac{\sqrt [3]{a}-2 \sqrt [3]{b} (c+d x)}{\sqrt{3} \sqrt [3]{a}}\right )}{\sqrt{3} a^{2/3} b^{2/3} d^2} \]

Antiderivative was successfully verified.

[In]

Int[x/(a + b*(c + d*x)^3),x]

[Out]

-(((a^(1/3) - b^(1/3)*c)*ArcTan[(a^(1/3) - 2*b^(1/3)*(c + d*x))/(Sqrt[3]*a^(1/3))])/(Sqrt[3]*a^(2/3)*b^(2/3)*d
^2)) - ((a^(1/3) + b^(1/3)*c)*Log[a^(1/3) + b^(1/3)*(c + d*x)])/(3*a^(2/3)*b^(2/3)*d^2) + ((a^(1/3) + b^(1/3)*
c)*Log[a^(2/3) - a^(1/3)*b^(1/3)*(c + d*x) + b^(2/3)*(c + d*x)^2])/(6*a^(2/3)*b^(2/3)*d^2)

Rule 371

Int[((a_) + (b_.)*(v_)^(n_))^(p_.)*(x_)^(m_.), x_Symbol] :> With[{c = Coefficient[v, x, 0], d = Coefficient[v,
 x, 1]}, Dist[1/d^(m + 1), Subst[Int[SimplifyIntegrand[(x - c)^m*(a + b*x^n)^p, x], x], x, v], x] /; NeQ[c, 0]
] /; FreeQ[{a, b, n, p}, x] && LinearQ[v, x] && IntegerQ[m]

Rule 1860

Int[((A_) + (B_.)*(x_))/((a_) + (b_.)*(x_)^3), x_Symbol] :> With[{r = Numerator[Rt[a/b, 3]], s = Denominator[R
t[a/b, 3]]}, -Dist[(r*(B*r - A*s))/(3*a*s), Int[1/(r + s*x), x], x] + Dist[r/(3*a*s), Int[(r*(B*r + 2*A*s) + s
*(B*r - A*s)*x)/(r^2 - r*s*x + s^2*x^2), x], x]] /; FreeQ[{a, b, A, B}, x] && NeQ[a*B^3 - b*A^3, 0] && PosQ[a/
b]

Rule 31

Int[((a_) + (b_.)*(x_))^(-1), x_Symbol] :> Simp[Log[RemoveContent[a + b*x, x]]/b, x] /; FreeQ[{a, b}, x]

Rule 634

Int[((d_.) + (e_.)*(x_))/((a_) + (b_.)*(x_) + (c_.)*(x_)^2), x_Symbol] :> Dist[(2*c*d - b*e)/(2*c), Int[1/(a +
 b*x + c*x^2), x], x] + Dist[e/(2*c), Int[(b + 2*c*x)/(a + b*x + c*x^2), x], x] /; FreeQ[{a, b, c, d, e}, x] &
& NeQ[2*c*d - b*e, 0] && NeQ[b^2 - 4*a*c, 0] &&  !NiceSqrtQ[b^2 - 4*a*c]

Rule 617

Int[((a_) + (b_.)*(x_) + (c_.)*(x_)^2)^(-1), x_Symbol] :> With[{q = 1 - 4*Simplify[(a*c)/b^2]}, Dist[-2/b, Sub
st[Int[1/(q - x^2), x], x, 1 + (2*c*x)/b], x] /; RationalQ[q] && (EqQ[q^2, 1] ||  !RationalQ[b^2 - 4*a*c])] /;
 FreeQ[{a, b, c}, x] && NeQ[b^2 - 4*a*c, 0]

Rule 204

Int[((a_) + (b_.)*(x_)^2)^(-1), x_Symbol] :> -Simp[ArcTan[(Rt[-b, 2]*x)/Rt[-a, 2]]/(Rt[-a, 2]*Rt[-b, 2]), x] /
; FreeQ[{a, b}, x] && PosQ[a/b] && (LtQ[a, 0] || LtQ[b, 0])

Rule 628

Int[((d_) + (e_.)*(x_))/((a_.) + (b_.)*(x_) + (c_.)*(x_)^2), x_Symbol] :> Simp[(d*Log[RemoveContent[a + b*x +
c*x^2, x]])/b, x] /; FreeQ[{a, b, c, d, e}, x] && EqQ[2*c*d - b*e, 0]

Rubi steps

\begin{align*} \int \frac{x}{a+b (c+d x)^3} \, dx &=\frac{\operatorname{Subst}\left (\int \frac{-c+x}{a+b x^3} \, dx,x,c+d x\right )}{d^2}\\ &=\frac{\operatorname{Subst}\left (\int \frac{\sqrt [3]{a} \left (\sqrt [3]{a}-2 \sqrt [3]{b} c\right )+\sqrt [3]{b} \left (\sqrt [3]{a}+\sqrt [3]{b} c\right ) x}{a^{2/3}-\sqrt [3]{a} \sqrt [3]{b} x+b^{2/3} x^2} \, dx,x,c+d x\right )}{3 a^{2/3} \sqrt [3]{b} d^2}-\frac{\left (\frac{\sqrt [3]{a}}{\sqrt [3]{b}}+c\right ) \operatorname{Subst}\left (\int \frac{1}{\sqrt [3]{a}+\sqrt [3]{b} x} \, dx,x,c+d x\right )}{3 a^{2/3} d^2}\\ &=-\frac{\left (\sqrt [3]{a}+\sqrt [3]{b} c\right ) \log \left (\sqrt [3]{a}+\sqrt [3]{b} (c+d x)\right )}{3 a^{2/3} b^{2/3} d^2}+\frac{\left (\frac{1}{\sqrt [3]{b}}-\frac{c}{\sqrt [3]{a}}\right ) \operatorname{Subst}\left (\int \frac{1}{a^{2/3}-\sqrt [3]{a} \sqrt [3]{b} x+b^{2/3} x^2} \, dx,x,c+d x\right )}{2 d^2}+\frac{\left (\sqrt [3]{a}+\sqrt [3]{b} c\right ) \operatorname{Subst}\left (\int \frac{-\sqrt [3]{a} \sqrt [3]{b}+2 b^{2/3} x}{a^{2/3}-\sqrt [3]{a} \sqrt [3]{b} x+b^{2/3} x^2} \, dx,x,c+d x\right )}{6 a^{2/3} b^{2/3} d^2}\\ &=-\frac{\left (\sqrt [3]{a}+\sqrt [3]{b} c\right ) \log \left (\sqrt [3]{a}+\sqrt [3]{b} (c+d x)\right )}{3 a^{2/3} b^{2/3} d^2}+\frac{\left (\sqrt [3]{a}+\sqrt [3]{b} c\right ) \log \left (a^{2/3}-\sqrt [3]{a} \sqrt [3]{b} (c+d x)+b^{2/3} (c+d x)^2\right )}{6 a^{2/3} b^{2/3} d^2}+\frac{\left (\sqrt [3]{a}-\sqrt [3]{b} c\right ) \operatorname{Subst}\left (\int \frac{1}{-3-x^2} \, dx,x,1-\frac{2 \sqrt [3]{b} (c+d x)}{\sqrt [3]{a}}\right )}{a^{2/3} b^{2/3} d^2}\\ &=-\frac{\left (\sqrt [3]{a}-\sqrt [3]{b} c\right ) \tan ^{-1}\left (\frac{1-\frac{2 \sqrt [3]{b} (c+d x)}{\sqrt [3]{a}}}{\sqrt{3}}\right )}{\sqrt{3} a^{2/3} b^{2/3} d^2}-\frac{\left (\sqrt [3]{a}+\sqrt [3]{b} c\right ) \log \left (\sqrt [3]{a}+\sqrt [3]{b} (c+d x)\right )}{3 a^{2/3} b^{2/3} d^2}+\frac{\left (\sqrt [3]{a}+\sqrt [3]{b} c\right ) \log \left (a^{2/3}-\sqrt [3]{a} \sqrt [3]{b} (c+d x)+b^{2/3} (c+d x)^2\right )}{6 a^{2/3} b^{2/3} d^2}\\ \end{align*}

Mathematica [C]  time = 0.0215535, size = 79, normalized size = 0.44 \[ \frac{\text{RootSum}\left [3 \text{$\#$1}^2 b c d^2+\text{$\#$1}^3 b d^3+3 \text{$\#$1} b c^2 d+a+b c^3\& ,\frac{\text{$\#$1} \log (x-\text{$\#$1})}{\text{$\#$1}^2 d^2+2 \text{$\#$1} c d+c^2}\& \right ]}{3 b d} \]

Antiderivative was successfully verified.

[In]

Integrate[x/(a + b*(c + d*x)^3),x]

[Out]

RootSum[a + b*c^3 + 3*b*c^2*d*#1 + 3*b*c*d^2*#1^2 + b*d^3*#1^3 & , (Log[x - #1]*#1)/(c^2 + 2*c*d*#1 + d^2*#1^2
) & ]/(3*b*d)

________________________________________________________________________________________

Maple [C]  time = 0.001, size = 72, normalized size = 0.4 \begin{align*}{\frac{1}{3\,bd}\sum _{{\it \_R}={\it RootOf} \left ({{\it \_Z}}^{3}b{d}^{3}+3\,{{\it \_Z}}^{2}bc{d}^{2}+3\,{\it \_Z}\,b{c}^{2}d+b{c}^{3}+a \right ) }{\frac{{\it \_R}\,\ln \left ( x-{\it \_R} \right ) }{{d}^{2}{{\it \_R}}^{2}+2\,cd{\it \_R}+{c}^{2}}}} \end{align*}

Verification of antiderivative is not currently implemented for this CAS.

[In]

int(x/(a+b*(d*x+c)^3),x)

[Out]

1/3/b/d*sum(_R/(_R^2*d^2+2*_R*c*d+c^2)*ln(x-_R),_R=RootOf(_Z^3*b*d^3+3*_Z^2*b*c*d^2+3*_Z*b*c^2*d+b*c^3+a))

________________________________________________________________________________________

Maxima [F]  time = 0., size = 0, normalized size = 0. \begin{align*} \int \frac{x}{{\left (d x + c\right )}^{3} b + a}\,{d x} \end{align*}

Verification of antiderivative is not currently implemented for this CAS.

[In]

integrate(x/(a+b*(d*x+c)^3),x, algorithm="maxima")

[Out]

integrate(x/((d*x + c)^3*b + a), x)

________________________________________________________________________________________

Fricas [C]  time = 7.86129, size = 4685, normalized size = 26.03 \begin{align*} \text{result too large to display} \end{align*}

Verification of antiderivative is not currently implemented for this CAS.

[In]

integrate(x/(a+b*(d*x+c)^3),x, algorithm="fricas")

[Out]

1/36*(9*(I*sqrt(3) + 1)*(-1/54*(b*c^3 + a)/(a^2*b^2*d^6) + 1/54*(b*c^3 - a)/(a^2*b^2*d^6))^(1/3) - 3*sqrt(1/3)
*sqrt(-((9*(I*sqrt(3) + 1)*(-1/54*(b*c^3 + a)/(a^2*b^2*d^6) + 1/54*(b*c^3 - a)/(a^2*b^2*d^6))^(1/3) + c*(-I*sq
rt(3) + 1)/(a*b*d^4*(-1/54*(b*c^3 + a)/(a^2*b^2*d^6) + 1/54*(b*c^3 - a)/(a^2*b^2*d^6))^(1/3)))^2*a*b*d^4 - 144
*c)/(a*b*d^4)) + c*(-I*sqrt(3) + 1)/(a*b*d^4*(-1/54*(b*c^3 + a)/(a^2*b^2*d^6) + 1/54*(b*c^3 - a)/(a^2*b^2*d^6)
)^(1/3)))*log(1/36*(9*(I*sqrt(3) + 1)*(-1/54*(b*c^3 + a)/(a^2*b^2*d^6) + 1/54*(b*c^3 - a)/(a^2*b^2*d^6))^(1/3)
 + c*(-I*sqrt(3) + 1)/(a*b*d^4*(-1/54*(b*c^3 + a)/(a^2*b^2*d^6) + 1/54*(b*c^3 - a)/(a^2*b^2*d^6))^(1/3)))^2*a^
2*b*d^4 - 1/6*(9*(I*sqrt(3) + 1)*(-1/54*(b*c^3 + a)/(a^2*b^2*d^6) + 1/54*(b*c^3 - a)/(a^2*b^2*d^6))^(1/3) + c*
(-I*sqrt(3) + 1)/(a*b*d^4*(-1/54*(b*c^3 + a)/(a^2*b^2*d^6) + 1/54*(b*c^3 - a)/(a^2*b^2*d^6))^(1/3)))*a*b*c^2*d
^2 + 2*b*c^4 + 2*(b*c^3 - a)*d*x - 4*a*c + 1/12*sqrt(1/3)*((9*(I*sqrt(3) + 1)*(-1/54*(b*c^3 + a)/(a^2*b^2*d^6)
 + 1/54*(b*c^3 - a)/(a^2*b^2*d^6))^(1/3) + c*(-I*sqrt(3) + 1)/(a*b*d^4*(-1/54*(b*c^3 + a)/(a^2*b^2*d^6) + 1/54
*(b*c^3 - a)/(a^2*b^2*d^6))^(1/3)))*a^2*b*d^4 + 6*a*b*c^2*d^2)*sqrt(-((9*(I*sqrt(3) + 1)*(-1/54*(b*c^3 + a)/(a
^2*b^2*d^6) + 1/54*(b*c^3 - a)/(a^2*b^2*d^6))^(1/3) + c*(-I*sqrt(3) + 1)/(a*b*d^4*(-1/54*(b*c^3 + a)/(a^2*b^2*
d^6) + 1/54*(b*c^3 - a)/(a^2*b^2*d^6))^(1/3)))^2*a*b*d^4 - 144*c)/(a*b*d^4))) + 1/36*(9*(I*sqrt(3) + 1)*(-1/54
*(b*c^3 + a)/(a^2*b^2*d^6) + 1/54*(b*c^3 - a)/(a^2*b^2*d^6))^(1/3) + 3*sqrt(1/3)*sqrt(-((9*(I*sqrt(3) + 1)*(-1
/54*(b*c^3 + a)/(a^2*b^2*d^6) + 1/54*(b*c^3 - a)/(a^2*b^2*d^6))^(1/3) + c*(-I*sqrt(3) + 1)/(a*b*d^4*(-1/54*(b*
c^3 + a)/(a^2*b^2*d^6) + 1/54*(b*c^3 - a)/(a^2*b^2*d^6))^(1/3)))^2*a*b*d^4 - 144*c)/(a*b*d^4)) + c*(-I*sqrt(3)
 + 1)/(a*b*d^4*(-1/54*(b*c^3 + a)/(a^2*b^2*d^6) + 1/54*(b*c^3 - a)/(a^2*b^2*d^6))^(1/3)))*log(1/36*(9*(I*sqrt(
3) + 1)*(-1/54*(b*c^3 + a)/(a^2*b^2*d^6) + 1/54*(b*c^3 - a)/(a^2*b^2*d^6))^(1/3) + c*(-I*sqrt(3) + 1)/(a*b*d^4
*(-1/54*(b*c^3 + a)/(a^2*b^2*d^6) + 1/54*(b*c^3 - a)/(a^2*b^2*d^6))^(1/3)))^2*a^2*b*d^4 - 1/6*(9*(I*sqrt(3) +
1)*(-1/54*(b*c^3 + a)/(a^2*b^2*d^6) + 1/54*(b*c^3 - a)/(a^2*b^2*d^6))^(1/3) + c*(-I*sqrt(3) + 1)/(a*b*d^4*(-1/
54*(b*c^3 + a)/(a^2*b^2*d^6) + 1/54*(b*c^3 - a)/(a^2*b^2*d^6))^(1/3)))*a*b*c^2*d^2 + 2*b*c^4 + 2*(b*c^3 - a)*d
*x - 4*a*c - 1/12*sqrt(1/3)*((9*(I*sqrt(3) + 1)*(-1/54*(b*c^3 + a)/(a^2*b^2*d^6) + 1/54*(b*c^3 - a)/(a^2*b^2*d
^6))^(1/3) + c*(-I*sqrt(3) + 1)/(a*b*d^4*(-1/54*(b*c^3 + a)/(a^2*b^2*d^6) + 1/54*(b*c^3 - a)/(a^2*b^2*d^6))^(1
/3)))*a^2*b*d^4 + 6*a*b*c^2*d^2)*sqrt(-((9*(I*sqrt(3) + 1)*(-1/54*(b*c^3 + a)/(a^2*b^2*d^6) + 1/54*(b*c^3 - a)
/(a^2*b^2*d^6))^(1/3) + c*(-I*sqrt(3) + 1)/(a*b*d^4*(-1/54*(b*c^3 + a)/(a^2*b^2*d^6) + 1/54*(b*c^3 - a)/(a^2*b
^2*d^6))^(1/3)))^2*a*b*d^4 - 144*c)/(a*b*d^4))) - 1/18*(9*(I*sqrt(3) + 1)*(-1/54*(b*c^3 + a)/(a^2*b^2*d^6) + 1
/54*(b*c^3 - a)/(a^2*b^2*d^6))^(1/3) + c*(-I*sqrt(3) + 1)/(a*b*d^4*(-1/54*(b*c^3 + a)/(a^2*b^2*d^6) + 1/54*(b*
c^3 - a)/(a^2*b^2*d^6))^(1/3)))*log(-1/36*(9*(I*sqrt(3) + 1)*(-1/54*(b*c^3 + a)/(a^2*b^2*d^6) + 1/54*(b*c^3 -
a)/(a^2*b^2*d^6))^(1/3) + c*(-I*sqrt(3) + 1)/(a*b*d^4*(-1/54*(b*c^3 + a)/(a^2*b^2*d^6) + 1/54*(b*c^3 - a)/(a^2
*b^2*d^6))^(1/3)))^2*a^2*b*d^4 + 1/6*(9*(I*sqrt(3) + 1)*(-1/54*(b*c^3 + a)/(a^2*b^2*d^6) + 1/54*(b*c^3 - a)/(a
^2*b^2*d^6))^(1/3) + c*(-I*sqrt(3) + 1)/(a*b*d^4*(-1/54*(b*c^3 + a)/(a^2*b^2*d^6) + 1/54*(b*c^3 - a)/(a^2*b^2*
d^6))^(1/3)))*a*b*c^2*d^2 + b*c^4 + (b*c^3 - a)*d*x + a*c)

________________________________________________________________________________________

Sympy [A]  time = 0.659871, size = 83, normalized size = 0.46 \begin{align*} \operatorname{RootSum}{\left (27 t^{3} a^{2} b^{2} d^{6} - 9 t a b c d^{2} + a + b c^{3}, \left ( t \mapsto t \log{\left (x + \frac{9 t^{2} a^{2} b d^{4} + 3 t a b c^{2} d^{2} - a c - b c^{4}}{a d - b c^{3} d} \right )} \right )\right )} \end{align*}

Verification of antiderivative is not currently implemented for this CAS.

[In]

integrate(x/(a+b*(d*x+c)**3),x)

[Out]

RootSum(27*_t**3*a**2*b**2*d**6 - 9*_t*a*b*c*d**2 + a + b*c**3, Lambda(_t, _t*log(x + (9*_t**2*a**2*b*d**4 + 3
*_t*a*b*c**2*d**2 - a*c - b*c**4)/(a*d - b*c**3*d))))

________________________________________________________________________________________

Giac [F]  time = 0., size = 0, normalized size = 0. \begin{align*} \int \frac{x}{{\left (d x + c\right )}^{3} b + a}\,{d x} \end{align*}

Verification of antiderivative is not currently implemented for this CAS.

[In]

integrate(x/(a+b*(d*x+c)^3),x, algorithm="giac")

[Out]

integrate(x/((d*x + c)^3*b + a), x)