\(\int \frac {a+b \log (c x^n)}{x^3 (d+e x)^3} \, dx\) [52]

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

Optimal result

Integrand size = 21, antiderivative size = 217 \[ \int \frac {a+b \log \left (c x^n\right )}{x^3 (d+e x)^3} \, dx=-\frac {b n}{4 d^3 x^2}+\frac {3 b e n}{d^4 x}-\frac {b e^2 n}{2 d^4 (d+e x)}-\frac {b e^2 n \log (x)}{2 d^5}-\frac {a+b \log \left (c x^n\right )}{2 d^3 x^2}+\frac {3 e \left (a+b \log \left (c x^n\right )\right )}{d^4 x}+\frac {e^2 \left (a+b \log \left (c x^n\right )\right )}{2 d^3 (d+e x)^2}-\frac {3 e^3 x \left (a+b \log \left (c x^n\right )\right )}{d^5 (d+e x)}-\frac {6 e^2 \log \left (1+\frac {d}{e x}\right ) \left (a+b \log \left (c x^n\right )\right )}{d^5}+\frac {7 b e^2 n \log (d+e x)}{2 d^5}+\frac {6 b e^2 n \operatorname {PolyLog}\left (2,-\frac {d}{e x}\right )}{d^5} \] Output:

-1/4*b*n/d^3/x^2+3*b*e*n/d^4/x-1/2*b*e^2*n/d^4/(e*x+d)-1/2*b*e^2*n*ln(x)/d 
^5-1/2*(a+b*ln(c*x^n))/d^3/x^2+3*e*(a+b*ln(c*x^n))/d^4/x+1/2*e^2*(a+b*ln(c 
*x^n))/d^3/(e*x+d)^2-3*e^3*x*(a+b*ln(c*x^n))/d^5/(e*x+d)-6*e^2*ln(1+d/e/x) 
*(a+b*ln(c*x^n))/d^5+7/2*b*e^2*n*ln(e*x+d)/d^5+6*b*e^2*n*polylog(2,-d/e/x) 
/d^5
 

Mathematica [A] (verified)

Time = 0.48 (sec) , antiderivative size = 227, normalized size of antiderivative = 1.05 \[ \int \frac {a+b \log \left (c x^n\right )}{x^3 (d+e x)^3} \, dx=-\frac {\frac {b d^2 n}{x^2}-\frac {12 b d e n}{x}+\frac {2 d^2 \left (a+b \log \left (c x^n\right )\right )}{x^2}-\frac {12 d e \left (a+b \log \left (c x^n\right )\right )}{x}-\frac {2 d^2 e^2 \left (a+b \log \left (c x^n\right )\right )}{(d+e x)^2}-\frac {12 d e^2 \left (a+b \log \left (c x^n\right )\right )}{d+e x}-\frac {12 e^2 \left (a+b \log \left (c x^n\right )\right )^2}{b n}+12 b e^2 n (\log (x)-\log (d+e x))+\frac {2 b e^2 n (d+(d+e x) \log (x)-(d+e x) \log (d+e x))}{d+e x}+24 e^2 \left (a+b \log \left (c x^n\right )\right ) \log \left (1+\frac {e x}{d}\right )+24 b e^2 n \operatorname {PolyLog}\left (2,-\frac {e x}{d}\right )}{4 d^5} \] Input:

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

Output:

-1/4*((b*d^2*n)/x^2 - (12*b*d*e*n)/x + (2*d^2*(a + b*Log[c*x^n]))/x^2 - (1 
2*d*e*(a + b*Log[c*x^n]))/x - (2*d^2*e^2*(a + b*Log[c*x^n]))/(d + e*x)^2 - 
 (12*d*e^2*(a + b*Log[c*x^n]))/(d + e*x) - (12*e^2*(a + b*Log[c*x^n])^2)/( 
b*n) + 12*b*e^2*n*(Log[x] - Log[d + e*x]) + (2*b*e^2*n*(d + (d + e*x)*Log[ 
x] - (d + e*x)*Log[d + e*x]))/(d + e*x) + 24*e^2*(a + b*Log[c*x^n])*Log[1 
+ (e*x)/d] + 24*b*e^2*n*PolyLog[2, -((e*x)/d)])/d^5
 

Rubi [A] (verified)

Time = 0.56 (sec) , antiderivative size = 217, 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.095, Rules used = {2793, 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 {a+b \log \left (c x^n\right )}{x^3 (d+e x)^3} \, dx\)

\(\Big \downarrow \) 2793

\(\displaystyle \int \left (-\frac {3 e^3 \left (a+b \log \left (c x^n\right )\right )}{d^4 (d+e x)^2}+\frac {6 e^2 \left (a+b \log \left (c x^n\right )\right )}{d^4 x (d+e x)}-\frac {3 e \left (a+b \log \left (c x^n\right )\right )}{d^4 x^2}-\frac {e^3 \left (a+b \log \left (c x^n\right )\right )}{d^3 (d+e x)^3}+\frac {a+b \log \left (c x^n\right )}{d^3 x^3}\right )dx\)

\(\Big \downarrow \) 2009

\(\displaystyle -\frac {3 e^3 x \left (a+b \log \left (c x^n\right )\right )}{d^5 (d+e x)}-\frac {6 e^2 \log \left (\frac {d}{e x}+1\right ) \left (a+b \log \left (c x^n\right )\right )}{d^5}+\frac {3 e \left (a+b \log \left (c x^n\right )\right )}{d^4 x}+\frac {e^2 \left (a+b \log \left (c x^n\right )\right )}{2 d^3 (d+e x)^2}-\frac {a+b \log \left (c x^n\right )}{2 d^3 x^2}+\frac {6 b e^2 n \operatorname {PolyLog}\left (2,-\frac {d}{e x}\right )}{d^5}-\frac {b e^2 n \log (x)}{2 d^5}+\frac {7 b e^2 n \log (d+e x)}{2 d^5}-\frac {b e^2 n}{2 d^4 (d+e x)}+\frac {3 b e n}{d^4 x}-\frac {b n}{4 d^3 x^2}\)

Input:

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

Output:

-1/4*(b*n)/(d^3*x^2) + (3*b*e*n)/(d^4*x) - (b*e^2*n)/(2*d^4*(d + e*x)) - ( 
b*e^2*n*Log[x])/(2*d^5) - (a + b*Log[c*x^n])/(2*d^3*x^2) + (3*e*(a + b*Log 
[c*x^n]))/(d^4*x) + (e^2*(a + b*Log[c*x^n]))/(2*d^3*(d + e*x)^2) - (3*e^3* 
x*(a + b*Log[c*x^n]))/(d^5*(d + e*x)) - (6*e^2*Log[1 + d/(e*x)]*(a + b*Log 
[c*x^n]))/d^5 + (7*b*e^2*n*Log[d + e*x])/(2*d^5) + (6*b*e^2*n*PolyLog[2, - 
(d/(e*x))])/d^5
 

Defintions of rubi rules used

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

rule 2793
Int[((a_.) + Log[(c_.)*(x_)^(n_.)]*(b_.))*((f_.)*(x_))^(m_.)*((d_) + (e_.)* 
(x_)^(r_.))^(q_.), x_Symbol] :> With[{u = ExpandIntegrand[a + b*Log[c*x^n], 
 (f*x)^m*(d + e*x^r)^q, x]}, Int[u, x] /; SumQ[u]] /; FreeQ[{a, b, c, d, e, 
 f, m, n, q, r}, x] && IntegerQ[q] && (GtQ[q, 0] || (IntegerQ[m] && Integer 
Q[r]))
 
Maple [C] (warning: unable to verify)

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

Time = 1.19 (sec) , antiderivative size = 386, normalized size of antiderivative = 1.78

method result size
risch \(-\frac {6 b \ln \left (x^{n}\right ) e^{2} \ln \left (e x +d \right )}{d^{5}}+\frac {3 b \ln \left (x^{n}\right ) e^{2}}{d^{4} \left (e x +d \right )}+\frac {b \ln \left (x^{n}\right ) e^{2}}{2 d^{3} \left (e x +d \right )^{2}}-\frac {b \ln \left (x^{n}\right )}{2 d^{3} x^{2}}+\frac {6 b \ln \left (x^{n}\right ) e^{2} \ln \left (x \right )}{d^{5}}+\frac {3 b \ln \left (x^{n}\right ) e}{d^{4} x}-\frac {b \,e^{2} n}{2 d^{4} \left (e x +d \right )}+\frac {7 b \,e^{2} n \ln \left (e x +d \right )}{2 d^{5}}-\frac {b n}{4 d^{3} x^{2}}+\frac {3 b e n}{d^{4} x}-\frac {7 b \,e^{2} n \ln \left (x \right )}{2 d^{5}}-\frac {3 b n \,e^{2} \ln \left (x \right )^{2}}{d^{5}}+\frac {6 b n \,e^{2} \ln \left (e x +d \right ) \ln \left (-\frac {e x}{d}\right )}{d^{5}}+\frac {6 b n \,e^{2} \operatorname {dilog}\left (-\frac {e x}{d}\right )}{d^{5}}+\left (\frac {i b \pi \,\operatorname {csgn}\left (i x^{n}\right ) \operatorname {csgn}\left (i c \,x^{n}\right )^{2}}{2}-\frac {i b \pi \,\operatorname {csgn}\left (i x^{n}\right ) \operatorname {csgn}\left (i c \,x^{n}\right ) \operatorname {csgn}\left (i c \right )}{2}-\frac {i b \pi \operatorname {csgn}\left (i c \,x^{n}\right )^{3}}{2}+\frac {i b \pi \operatorname {csgn}\left (i c \,x^{n}\right )^{2} \operatorname {csgn}\left (i c \right )}{2}+b \ln \left (c \right )+a \right ) \left (-\frac {6 e^{2} \ln \left (e x +d \right )}{d^{5}}+\frac {3 e^{2}}{d^{4} \left (e x +d \right )}+\frac {e^{2}}{2 d^{3} \left (e x +d \right )^{2}}-\frac {1}{2 d^{3} x^{2}}+\frac {6 e^{2} \ln \left (x \right )}{d^{5}}+\frac {3 e}{d^{4} x}\right )\) \(386\)

Input:

int((a+b*ln(c*x^n))/x^3/(e*x+d)^3,x,method=_RETURNVERBOSE)
 

Output:

-6*b*ln(x^n)/d^5*e^2*ln(e*x+d)+3*b*ln(x^n)/d^4*e^2/(e*x+d)+1/2*b*ln(x^n)/d 
^3*e^2/(e*x+d)^2-1/2*b*ln(x^n)/d^3/x^2+6*b*ln(x^n)/d^5*e^2*ln(x)+3*b*ln(x^ 
n)/d^4*e/x-1/2*b*e^2*n/d^4/(e*x+d)+7/2*b*e^2*n*ln(e*x+d)/d^5-1/4*b*n/d^3/x 
^2+3*b*e*n/d^4/x-7/2*b*e^2*n*ln(x)/d^5-3*b*n/d^5*e^2*ln(x)^2+6*b*n/d^5*e^2 
*ln(e*x+d)*ln(-e*x/d)+6*b*n/d^5*e^2*dilog(-e*x/d)+(1/2*I*Pi*b*csgn(I*x^n)* 
csgn(I*c*x^n)^2-1/2*I*Pi*b*csgn(I*x^n)*csgn(I*c*x^n)*csgn(I*c)-1/2*I*Pi*b* 
csgn(I*c*x^n)^3+1/2*I*Pi*b*csgn(I*c*x^n)^2*csgn(I*c)+b*ln(c)+a)*(-6/d^5*e^ 
2*ln(e*x+d)+3/d^4*e^2/(e*x+d)+1/2/d^3*e^2/(e*x+d)^2-1/2/d^3/x^2+6/d^5*e^2* 
ln(x)+3/d^4*e/x)
 

Fricas [F]

\[ \int \frac {a+b \log \left (c x^n\right )}{x^3 (d+e x)^3} \, dx=\int { \frac {b \log \left (c x^{n}\right ) + a}{{\left (e x + d\right )}^{3} x^{3}} \,d x } \] Input:

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

Output:

integral((b*log(c*x^n) + a)/(e^3*x^6 + 3*d*e^2*x^5 + 3*d^2*e*x^4 + d^3*x^3 
), x)
 

Sympy [A] (verification not implemented)

Time = 50.00 (sec) , antiderivative size = 496, normalized size of antiderivative = 2.29 \[ \int \frac {a+b \log \left (c x^n\right )}{x^3 (d+e x)^3} \, dx =\text {Too large to display} \] Input:

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

Output:

-a*e**3*Piecewise((x/d**3, Eq(e, 0)), (-1/(2*e*(d + e*x)**2), True))/d**3 
- a/(2*d**3*x**2) - 3*a*e**3*Piecewise((x/d**2, Eq(e, 0)), (-1/(d*e + e**2 
*x), True))/d**4 + 3*a*e/(d**4*x) - 6*a*e**3*Piecewise((x/d, Eq(e, 0)), (l 
og(d + e*x)/e, True))/d**5 + 6*a*e**2*log(x)/d**5 + b*e**3*n*Piecewise((x/ 
d**3, Eq(e, 0)), (-1/(2*d**2*e + 2*d*e**2*x) - log(x)/(2*d**2*e) + log(d/e 
 + x)/(2*d**2*e), True))/d**3 - b*e**3*Piecewise((x/d**3, Eq(e, 0)), (-1/( 
2*e*(d + e*x)**2), True))*log(c*x**n)/d**3 - b*n/(4*d**3*x**2) - b*log(c*x 
**n)/(2*d**3*x**2) + 3*b*e**3*n*Piecewise((x/d**2, Eq(e, 0)), (-log(x)/(d* 
e) + log(d/e + x)/(d*e), True))/d**4 - 3*b*e**3*Piecewise((x/d**2, Eq(e, 0 
)), (-1/(d*e + e**2*x), True))*log(c*x**n)/d**4 + 3*b*e*n/(d**4*x) + 3*b*e 
*log(c*x**n)/(d**4*x) + 6*b*e**3*n*Piecewise((x/d, Eq(e, 0)), (Piecewise(( 
-polylog(2, e*x*exp_polar(I*pi)/d), (Abs(x) < 1) & (1/Abs(x) < 1)), (log(d 
)*log(x) - polylog(2, e*x*exp_polar(I*pi)/d), Abs(x) < 1), (-log(d)*log(1/ 
x) - polylog(2, e*x*exp_polar(I*pi)/d), 1/Abs(x) < 1), (-meijerg(((), (1, 
1)), ((0, 0), ()), x)*log(d) + meijerg(((1, 1), ()), ((), (0, 0)), x)*log( 
d) - polylog(2, e*x*exp_polar(I*pi)/d), True))/e, True))/d**5 - 6*b*e**3*P 
iecewise((x/d, Eq(e, 0)), (log(d + e*x)/e, True))*log(c*x**n)/d**5 - 3*b*e 
**2*n*log(x)**2/d**5 + 6*b*e**2*log(x)*log(c*x**n)/d**5
 

Maxima [F]

\[ \int \frac {a+b \log \left (c x^n\right )}{x^3 (d+e x)^3} \, dx=\int { \frac {b \log \left (c x^{n}\right ) + a}{{\left (e x + d\right )}^{3} x^{3}} \,d x } \] Input:

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

Output:

1/2*a*((12*e^3*x^3 + 18*d*e^2*x^2 + 4*d^2*e*x - d^3)/(d^4*e^2*x^4 + 2*d^5* 
e*x^3 + d^6*x^2) - 12*e^2*log(e*x + d)/d^5 + 12*e^2*log(x)/d^5) + b*integr 
ate((log(c) + log(x^n))/(e^3*x^6 + 3*d*e^2*x^5 + 3*d^2*e*x^4 + d^3*x^3), x 
)
 

Giac [F]

\[ \int \frac {a+b \log \left (c x^n\right )}{x^3 (d+e x)^3} \, dx=\int { \frac {b \log \left (c x^{n}\right ) + a}{{\left (e x + d\right )}^{3} x^{3}} \,d x } \] Input:

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

Output:

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

Mupad [F(-1)]

Timed out. \[ \int \frac {a+b \log \left (c x^n\right )}{x^3 (d+e x)^3} \, dx=\int \frac {a+b\,\ln \left (c\,x^n\right )}{x^3\,{\left (d+e\,x\right )}^3} \,d x \] Input:

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

Output:

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

Reduce [F]

\[ \int \frac {a+b \log \left (c x^n\right )}{x^3 (d+e x)^3} \, dx=\frac {2 \left (\int \frac {\mathrm {log}\left (x^{n} c \right )}{e^{3} x^{6}+3 d \,e^{2} x^{5}+3 d^{2} e \,x^{4}+d^{3} x^{3}}d x \right ) b \,d^{7} x^{2}+4 \left (\int \frac {\mathrm {log}\left (x^{n} c \right )}{e^{3} x^{6}+3 d \,e^{2} x^{5}+3 d^{2} e \,x^{4}+d^{3} x^{3}}d x \right ) b \,d^{6} e \,x^{3}+2 \left (\int \frac {\mathrm {log}\left (x^{n} c \right )}{e^{3} x^{6}+3 d \,e^{2} x^{5}+3 d^{2} e \,x^{4}+d^{3} x^{3}}d x \right ) b \,d^{5} e^{2} x^{4}-12 \,\mathrm {log}\left (e x +d \right ) a \,d^{2} e^{2} x^{2}-24 \,\mathrm {log}\left (e x +d \right ) a d \,e^{3} x^{3}-12 \,\mathrm {log}\left (e x +d \right ) a \,e^{4} x^{4}+12 \,\mathrm {log}\left (x \right ) a \,d^{2} e^{2} x^{2}+24 \,\mathrm {log}\left (x \right ) a d \,e^{3} x^{3}+12 \,\mathrm {log}\left (x \right ) a \,e^{4} x^{4}-a \,d^{4}+4 a \,d^{3} e x +12 a \,d^{2} e^{2} x^{2}-6 a \,e^{4} x^{4}}{2 d^{5} x^{2} \left (e^{2} x^{2}+2 d e x +d^{2}\right )} \] Input:

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

Output:

(2*int(log(x**n*c)/(d**3*x**3 + 3*d**2*e*x**4 + 3*d*e**2*x**5 + e**3*x**6) 
,x)*b*d**7*x**2 + 4*int(log(x**n*c)/(d**3*x**3 + 3*d**2*e*x**4 + 3*d*e**2* 
x**5 + e**3*x**6),x)*b*d**6*e*x**3 + 2*int(log(x**n*c)/(d**3*x**3 + 3*d**2 
*e*x**4 + 3*d*e**2*x**5 + e**3*x**6),x)*b*d**5*e**2*x**4 - 12*log(d + e*x) 
*a*d**2*e**2*x**2 - 24*log(d + e*x)*a*d*e**3*x**3 - 12*log(d + e*x)*a*e**4 
*x**4 + 12*log(x)*a*d**2*e**2*x**2 + 24*log(x)*a*d*e**3*x**3 + 12*log(x)*a 
*e**4*x**4 - a*d**4 + 4*a*d**3*e*x + 12*a*d**2*e**2*x**2 - 6*a*e**4*x**4)/ 
(2*d**5*x**2*(d**2 + 2*d*e*x + e**2*x**2))