\(\int \frac {a+b \log (c x^n)}{x^2 (d+e \log (f x^m))} \, dx\) [180]

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

Optimal result

Integrand size = 26, antiderivative size = 133 \[ \int \frac {a+b \log \left (c x^n\right )}{x^2 \left (d+e \log \left (f x^m\right )\right )} \, dx=-\frac {b n}{e m x}-\frac {b e^{\frac {d}{e m}} n \left (f x^m\right )^{\frac {1}{m}} \operatorname {ExpIntegralEi}\left (-\frac {d+e \log \left (f x^m\right )}{e m}\right ) \left (d+e \log \left (f x^m\right )\right )}{e^2 m^2 x}+\frac {e^{\frac {d}{e m}} \left (f x^m\right )^{\frac {1}{m}} \operatorname {ExpIntegralEi}\left (-\frac {d+e \log \left (f x^m\right )}{e m}\right ) \left (a+b \log \left (c x^n\right )\right )}{e m x} \] Output:

-b*n/e/m/x-b*exp(d/e/m)*n*(f*x^m)^(1/m)*Ei(-(d+e*ln(f*x^m))/e/m)*(d+e*ln(f 
*x^m))/e^2/m^2/x+exp(d/e/m)*(f*x^m)^(1/m)*Ei(-(d+e*ln(f*x^m))/e/m)*(a+b*ln 
(c*x^n))/e/m/x
 

Mathematica [A] (verified)

Time = 0.15 (sec) , antiderivative size = 87, normalized size of antiderivative = 0.65 \[ \int \frac {a+b \log \left (c x^n\right )}{x^2 \left (d+e \log \left (f x^m\right )\right )} \, dx=\frac {-b e m n+e^{\frac {d}{e m}} \left (f x^m\right )^{\frac {1}{m}} \operatorname {ExpIntegralEi}\left (-\frac {d+e \log \left (f x^m\right )}{e m}\right ) \left (a e m-b d n-b e n \log \left (f x^m\right )+b e m \log \left (c x^n\right )\right )}{e^2 m^2 x} \] Input:

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

Output:

(-(b*e*m*n) + E^(d/(e*m))*(f*x^m)^m^(-1)*ExpIntegralEi[-((d + e*Log[f*x^m] 
)/(e*m))]*(a*e*m - b*d*n - b*e*n*Log[f*x^m] + b*e*m*Log[c*x^n]))/(e^2*m^2* 
x)
 

Rubi [A] (warning: unable to verify)

Time = 0.64 (sec) , antiderivative size = 141, normalized size of antiderivative = 1.06, number of steps used = 7, number of rules used = 6, \(\frac {\text {number of rules}}{\text {integrand size}}\) = 0.231, Rules used = {2813, 27, 31, 3039, 7281, 7036}

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^2 \left (d+e \log \left (f x^m\right )\right )} \, dx\)

\(\Big \downarrow \) 2813

\(\displaystyle \frac {e^{\frac {d}{e m}} \left (f x^m\right )^{\frac {1}{m}} \left (a+b \log \left (c x^n\right )\right ) \operatorname {ExpIntegralEi}\left (-\frac {d+e \log \left (f x^m\right )}{e m}\right )}{e m x}-b n \int \frac {e^{\frac {d}{e m}} \left (f x^m\right )^{\frac {1}{m}} \operatorname {ExpIntegralEi}\left (-\frac {d+e \log \left (f x^m\right )}{e m}\right )}{e m x^2}dx\)

\(\Big \downarrow \) 27

\(\displaystyle \frac {e^{\frac {d}{e m}} \left (f x^m\right )^{\frac {1}{m}} \left (a+b \log \left (c x^n\right )\right ) \operatorname {ExpIntegralEi}\left (-\frac {d+e \log \left (f x^m\right )}{e m}\right )}{e m x}-\frac {b n e^{\frac {d}{e m}} \int \frac {\left (f x^m\right )^{\frac {1}{m}} \operatorname {ExpIntegralEi}\left (-\frac {d+e \log \left (f x^m\right )}{e m}\right )}{x^2}dx}{e m}\)

\(\Big \downarrow \) 31

\(\displaystyle \frac {e^{\frac {d}{e m}} \left (f x^m\right )^{\frac {1}{m}} \left (a+b \log \left (c x^n\right )\right ) \operatorname {ExpIntegralEi}\left (-\frac {d+e \log \left (f x^m\right )}{e m}\right )}{e m x}-\frac {b n e^{\frac {d}{e m}} \left (f x^m\right )^{\frac {1}{m}} \int \frac {\operatorname {ExpIntegralEi}\left (-\frac {d+e \log \left (f x^m\right )}{e m}\right )}{x}dx}{e m x}\)

\(\Big \downarrow \) 3039

\(\displaystyle \frac {e^{\frac {d}{e m}} \left (f x^m\right )^{\frac {1}{m}} \left (a+b \log \left (c x^n\right )\right ) \operatorname {ExpIntegralEi}\left (-\frac {d+e \log \left (f x^m\right )}{e m}\right )}{e m x}-\frac {b n e^{\frac {d}{e m}} \left (f x^m\right )^{\frac {1}{m}} \int \operatorname {ExpIntegralEi}\left (-\frac {d+e \log \left (f x^m\right )}{e m}\right )d\log \left (f x^m\right )}{e m^2 x}\)

\(\Big \downarrow \) 7281

\(\displaystyle \frac {b n e^{\frac {d}{e m}} \left (f x^m\right )^{\frac {1}{m}} \int \operatorname {ExpIntegralEi}\left (-\frac {d}{e m}-\frac {\log \left (f x^m\right )}{m}\right )d\left (-\frac {d}{e m}-\frac {\log \left (f x^m\right )}{m}\right )}{e m x}+\frac {e^{\frac {d}{e m}} \left (f x^m\right )^{\frac {1}{m}} \left (a+b \log \left (c x^n\right )\right ) \operatorname {ExpIntegralEi}\left (-\frac {d+e \log \left (f x^m\right )}{e m}\right )}{e m x}\)

\(\Big \downarrow \) 7036

\(\displaystyle \frac {e^{\frac {d}{e m}} \left (f x^m\right )^{\frac {1}{m}} \left (a+b \log \left (c x^n\right )\right ) \operatorname {ExpIntegralEi}\left (-\frac {d+e \log \left (f x^m\right )}{e m}\right )}{e m x}+\frac {b n e^{\frac {d}{e m}} \left (f x^m\right )^{\frac {1}{m}} \left (\left (-\frac {d}{e m}-\frac {\log \left (f x^m\right )}{m}\right ) \operatorname {ExpIntegralEi}\left (-\frac {d}{e m}-\frac {\log \left (f x^m\right )}{m}\right )-f x^m\right )}{e m x}\)

Input:

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

Output:

(b*E^(d/(e*m))*n*(f*x^m)^m^(-1)*(-(f*x^m) + ExpIntegralEi[-(d/(e*m)) - Log 
[f*x^m]/m]*(-(d/(e*m)) - Log[f*x^m]/m)))/(e*m*x) + (E^(d/(e*m))*(f*x^m)^m^ 
(-1)*ExpIntegralEi[-((d + e*Log[f*x^m])/(e*m))]*(a + b*Log[c*x^n]))/(e*m*x 
)
 

Defintions of rubi rules used

rule 27
Int[(a_)*(Fx_), x_Symbol] :> Simp[a   Int[Fx, x], x] /; FreeQ[a, x] &&  !Ma 
tchQ[Fx, (b_)*(Gx_) /; FreeQ[b, x]]
 

rule 31
Int[(u_.)*((a_.)*(x_))^(m_.)*((b_.)*(x_)^(i_.))^(p_), x_Symbol] :> Simp[(b* 
x^i)^p/(a*x)^(i*p)   Int[u*(a*x)^(m + i*p), x], x] /; FreeQ[{a, b, i, m, p} 
, x] &&  !IntegerQ[p]
 

rule 2813
Int[((a_.) + Log[(c_.)*(x_)^(n_.)]*(b_.))^(p_.)*((d_.) + Log[(f_.)*(x_)^(r_ 
.)]*(e_.))*((g_.)*(x_))^(m_.), x_Symbol] :> With[{u = IntHide[(g*x)^m*(a + 
b*Log[c*x^n])^p, x]}, Simp[(d + e*Log[f*x^r])   u, x] - Simp[e*r   Int[Simp 
lifyIntegrand[u/x, x], x], x]] /; FreeQ[{a, b, c, d, e, f, g, m, n, p, r}, 
x] &&  !(EqQ[p, 1] && EqQ[a, 0] && NeQ[d, 0])
 

rule 3039
Int[u_, x_Symbol] :> With[{lst = FunctionOfLog[Cancel[x*u], x]}, Simp[1/lst 
[[3]]   Subst[Int[lst[[1]], x], x, Log[lst[[2]]]], x] /;  !FalseQ[lst]] /; 
NonsumQ[u]
 

rule 7036
Int[ExpIntegralEi[(a_.) + (b_.)*(x_)], x_Symbol] :> Simp[(a + b*x)*(ExpInte 
gralEi[a + b*x]/b), x] - Simp[E^(a + b*x)/b, x] /; FreeQ[{a, b}, x]
 

rule 7281
Int[u_, x_Symbol] :> With[{lst = FunctionOfLinear[u, x]}, Simp[1/lst[[3]] 
 Subst[Int[lst[[1]], x], x, lst[[2]] + lst[[3]]*x], x] /;  !FalseQ[lst]]
 
Maple [C] (warning: unable to verify)

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

Time = 2.64 (sec) , antiderivative size = 2296, normalized size of antiderivative = 17.26

method result size
risch \(\text {Expression too large to display}\) \(2296\)

Input:

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

Output:

-1/2*(I*Pi*b*csgn(I*x^n)*csgn(I*c*x^n)^2-I*Pi*b*csgn(I*x^n)*csgn(I*c*x^n)* 
csgn(I*c)-I*Pi*b*csgn(I*c*x^n)^3+I*Pi*b*csgn(I*c*x^n)^2*csgn(I*c)+2*b*ln(c 
)+2*a)/m/e/x*f^(1/m)*(x^m)^(1/m)*exp(1/2*(-I*Pi*csgn(I*f)*csgn(I*x^m)*csgn 
(I*f*x^m)*e+I*Pi*csgn(I*f)*csgn(I*f*x^m)^2*e+I*Pi*csgn(I*x^m)*csgn(I*f*x^m 
)^2*e-I*Pi*csgn(I*f*x^m)^3*e+2*d)/m/e)*Ei(1,ln(x)-1/2*I*(e*Pi*csgn(I*f)*cs 
gn(I*x^m)*csgn(I*f*x^m)-e*Pi*csgn(I*f)*csgn(I*f*x^m)^2-e*Pi*csgn(I*x^m)*cs 
gn(I*f*x^m)^2+e*Pi*csgn(I*f*x^m)^3+2*I*e*ln(f)+2*I*e*(ln(x^m)-m*ln(x))+2*I 
*d)/m/e)-b/m/e/x*f^(1/m)*(x^m)^(1/m)*exp(1/2*(-I*Pi*csgn(I*f)*csgn(I*x^m)* 
csgn(I*f*x^m)*e+I*Pi*csgn(I*f)*csgn(I*f*x^m)^2*e+I*Pi*csgn(I*x^m)*csgn(I*f 
*x^m)^2*e-I*Pi*csgn(I*f*x^m)^3*e+2*d)/m/e)*Ei(1,ln(x)-1/2*I*(e*Pi*csgn(I*f 
)*csgn(I*x^m)*csgn(I*f*x^m)-e*Pi*csgn(I*f)*csgn(I*f*x^m)^2-e*Pi*csgn(I*x^m 
)*csgn(I*f*x^m)^2+e*Pi*csgn(I*f*x^m)^3+2*I*e*ln(f)+2*I*e*(ln(x^m)-m*ln(x)) 
+2*I*d)/m/e)*ln(x^n)-b*n/e/m/x-1/2*I*b*n/m^2/e/x*f^(1/m)*(x^m)^(1/m)*exp(1 
/2*(-I*Pi*csgn(I*f)*csgn(I*x^m)*csgn(I*f*x^m)*e+I*Pi*csgn(I*f)*csgn(I*f*x^ 
m)^2*e+I*Pi*csgn(I*x^m)*csgn(I*f*x^m)^2*e-I*Pi*csgn(I*f*x^m)^3*e+2*d)/m/e) 
*Ei(1,ln(x)-1/2*I*(e*Pi*csgn(I*f)*csgn(I*x^m)*csgn(I*f*x^m)-e*Pi*csgn(I*f) 
*csgn(I*f*x^m)^2-e*Pi*csgn(I*x^m)*csgn(I*f*x^m)^2+e*Pi*csgn(I*f*x^m)^3+2*I 
*e*ln(f)+2*I*e*(ln(x^m)-m*ln(x))+2*I*d)/m/e)*Pi*csgn(I*f)*csgn(I*x^m)*csgn 
(I*f*x^m)+1/2*I*b*n/m^2/e/x*f^(1/m)*(x^m)^(1/m)*exp(1/2*(-I*Pi*csgn(I*f)*c 
sgn(I*x^m)*csgn(I*f*x^m)*e+I*Pi*csgn(I*f)*csgn(I*f*x^m)^2*e+I*Pi*csgn(I...
 

Fricas [A] (verification not implemented)

Time = 0.07 (sec) , antiderivative size = 81, normalized size of antiderivative = 0.61 \[ \int \frac {a+b \log \left (c x^n\right )}{x^2 \left (d+e \log \left (f x^m\right )\right )} \, dx=-\frac {b e m n - {\left (b e m x \log \left (c\right ) - b e n x \log \left (f\right ) + {\left (a e m - b d n\right )} x\right )} e^{\left (\frac {e \log \left (f\right ) + d}{e m}\right )} \operatorname {log\_integral}\left (\frac {e^{\left (-\frac {e \log \left (f\right ) + d}{e m}\right )}}{x}\right )}{e^{2} m^{2} x} \] Input:

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

Output:

-(b*e*m*n - (b*e*m*x*log(c) - b*e*n*x*log(f) + (a*e*m - b*d*n)*x)*e^((e*lo 
g(f) + d)/(e*m))*log_integral(e^(-(e*log(f) + d)/(e*m))/x))/(e^2*m^2*x)
 

Sympy [F]

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

integrate((a+b*ln(c*x**n))/x**2/(d+e*ln(f*x**m)),x)
 

Output:

Integral((a + b*log(c*x**n))/(x**2*(d + e*log(f*x**m))), x)
 

Maxima [F]

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

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

Output:

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

Giac [F]

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

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

Output:

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

Mupad [F(-1)]

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

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

Output:

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

Reduce [F]

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

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

Output:

int(log(x**n*c)/(log(x**m*f)*e*x**2 + d*x**2),x)*b + int(1/(log(x**m*f)*e* 
x**2 + d*x**2),x)*a