\(\int \frac {\text {arctanh}(a+b x)^2}{x^3} \, dx\) [7]

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

Optimal result

Integrand size = 12, antiderivative size = 370 \[ \int \frac {\text {arctanh}(a+b x)^2}{x^3} \, dx=-\frac {b \text {arctanh}(a+b x)}{\left (1-a^2\right ) x}-\frac {\text {arctanh}(a+b x)^2}{2 x^2}+\frac {b^2 \log (x)}{\left (1-a^2\right )^2}+\frac {b^2 \text {arctanh}(a+b x) \log \left (\frac {2}{1-a-b x}\right )}{2 (1-a)^2}-\frac {b^2 \log (1-a-b x)}{2 (1-a)^2 (1+a)}-\frac {b^2 \text {arctanh}(a+b x) \log \left (\frac {2}{1+a+b x}\right )}{2 (1+a)^2}-\frac {2 a b^2 \text {arctanh}(a+b x) \log \left (\frac {2}{1+a+b x}\right )}{\left (1-a^2\right )^2}+\frac {2 a b^2 \text {arctanh}(a+b x) \log \left (\frac {2 b x}{(1-a) (1+a+b x)}\right )}{\left (1-a^2\right )^2}-\frac {b^2 \log (1+a+b x)}{2 (1-a) (1+a)^2}+\frac {b^2 \operatorname {PolyLog}\left (2,-\frac {1+a+b x}{1-a-b x}\right )}{4 (1-a)^2}+\frac {b^2 \operatorname {PolyLog}\left (2,1-\frac {2}{1+a+b x}\right )}{4 (1+a)^2}+\frac {a b^2 \operatorname {PolyLog}\left (2,1-\frac {2}{1+a+b x}\right )}{\left (1-a^2\right )^2}-\frac {a b^2 \operatorname {PolyLog}\left (2,1-\frac {2 b x}{(1-a) (1+a+b x)}\right )}{\left (1-a^2\right )^2} \] Output:

-b*arctanh(b*x+a)/(-a^2+1)/x-1/2*arctanh(b*x+a)^2/x^2+b^2*ln(x)/(-a^2+1)^2 
+1/2*b^2*arctanh(b*x+a)*ln(2/(-b*x-a+1))/(1-a)^2-1/2*b^2*ln(-b*x-a+1)/(1-a 
)^2/(1+a)-1/2*b^2*arctanh(b*x+a)*ln(2/(b*x+a+1))/(1+a)^2-2*a*b^2*arctanh(b 
*x+a)*ln(2/(b*x+a+1))/(-a^2+1)^2+2*a*b^2*arctanh(b*x+a)*ln(2*b*x/(1-a)/(b* 
x+a+1))/(-a^2+1)^2-1/2*b^2*ln(b*x+a+1)/(1-a)/(1+a)^2+1/4*b^2*polylog(2,-(b 
*x+a+1)/(-b*x-a+1))/(1-a)^2+1/4*b^2*polylog(2,1-2/(b*x+a+1))/(1+a)^2+a*b^2 
*polylog(2,1-2/(b*x+a+1))/(-a^2+1)^2-a*b^2*polylog(2,1-2*b*x/(1-a)/(b*x+a+ 
1))/(-a^2+1)^2
 

Mathematica [C] (warning: unable to verify)

Result contains complex when optimal does not.

Time = 1.59 (sec) , antiderivative size = 271, normalized size of antiderivative = 0.73 \[ \int \frac {\text {arctanh}(a+b x)^2}{x^3} \, dx=\frac {-\left (\left (1+a^4-b^2 \left (-1+2 \sqrt {1-a^2} e^{\text {arctanh}(a)}\right ) x^2-a^2 \left (2+b^2 x^2\right )\right ) \text {arctanh}(a+b x)^2\right )+2 b x \text {arctanh}(a+b x) \left (-1+a^2+a b x+i a b \pi x-2 a b x \text {arctanh}(a)+2 a b x \log \left (1-e^{2 \text {arctanh}(a)-2 \text {arctanh}(a+b x)}\right )\right )+2 b^2 x^2 \left (-i a \pi \log \left (1+e^{2 \text {arctanh}(a+b x)}\right )+i a \pi \log \left (\frac {1}{\sqrt {1-(a+b x)^2}}\right )+\log \left (-\frac {b x}{\sqrt {1-(a+b x)^2}}\right )-2 a \text {arctanh}(a) \left (\log \left (1-e^{2 \text {arctanh}(a)-2 \text {arctanh}(a+b x)}\right )-\log (-i \sinh (\text {arctanh}(a)-\text {arctanh}(a+b x)))\right )\right )-2 a b^2 x^2 \operatorname {PolyLog}\left (2,e^{2 \text {arctanh}(a)-2 \text {arctanh}(a+b x)}\right )}{2 \left (-1+a^2\right )^2 x^2} \] Input:

Integrate[ArcTanh[a + b*x]^2/x^3,x]
 

Output:

(-((1 + a^4 - b^2*(-1 + 2*Sqrt[1 - a^2]*E^ArcTanh[a])*x^2 - a^2*(2 + b^2*x 
^2))*ArcTanh[a + b*x]^2) + 2*b*x*ArcTanh[a + b*x]*(-1 + a^2 + a*b*x + I*a* 
b*Pi*x - 2*a*b*x*ArcTanh[a] + 2*a*b*x*Log[1 - E^(2*ArcTanh[a] - 2*ArcTanh[ 
a + b*x])]) + 2*b^2*x^2*((-I)*a*Pi*Log[1 + E^(2*ArcTanh[a + b*x])] + I*a*P 
i*Log[1/Sqrt[1 - (a + b*x)^2]] + Log[-((b*x)/Sqrt[1 - (a + b*x)^2])] - 2*a 
*ArcTanh[a]*(Log[1 - E^(2*ArcTanh[a] - 2*ArcTanh[a + b*x])] - Log[(-I)*Sin 
h[ArcTanh[a] - ArcTanh[a + b*x]]])) - 2*a*b^2*x^2*PolyLog[2, E^(2*ArcTanh[ 
a] - 2*ArcTanh[a + b*x])])/(2*(-1 + a^2)^2*x^2)
 

Rubi [A] (verified)

Time = 1.11 (sec) , antiderivative size = 347, normalized size of antiderivative = 0.94, number of steps used = 7, number of rules used = 6, \(\frac {\text {number of rules}}{\text {integrand size}}\) = 0.500, Rules used = {6659, 7292, 6671, 27, 7276, 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 {\text {arctanh}(a+b x)^2}{x^3} \, dx\)

\(\Big \downarrow \) 6659

\(\displaystyle b \int \frac {\text {arctanh}(a+b x)}{x^2 \left (1-(a+b x)^2\right )}dx-\frac {\text {arctanh}(a+b x)^2}{2 x^2}\)

\(\Big \downarrow \) 7292

\(\displaystyle b \int \frac {\text {arctanh}(a+b x)}{x^2 \left (-a^2-2 b x a-b^2 x^2+1\right )}dx-\frac {\text {arctanh}(a+b x)^2}{2 x^2}\)

\(\Big \downarrow \) 6671

\(\displaystyle \int \frac {\text {arctanh}(a+b x)}{x^2 \left (1-(a+b x)^2\right )}d(a+b x)-\frac {\text {arctanh}(a+b x)^2}{2 x^2}\)

\(\Big \downarrow \) 27

\(\displaystyle b^2 \int \frac {\text {arctanh}(a+b x)}{b^2 x^2 \left (1-(a+b x)^2\right )}d(a+b x)-\frac {\text {arctanh}(a+b x)^2}{2 x^2}\)

\(\Big \downarrow \) 7276

\(\displaystyle b^2 \int \left (\frac {2 a \text {arctanh}(a+b x)}{\left (a^2-1\right )^2 b x}-\frac {\text {arctanh}(a+b x)}{2 (a-1)^2 (a+b x-1)}+\frac {\text {arctanh}(a+b x)}{2 (a+1)^2 (a+b x+1)}-\frac {\text {arctanh}(a+b x)}{\left (a^2-1\right ) b^2 x^2}\right )d(a+b x)-\frac {\text {arctanh}(a+b x)^2}{2 x^2}\)

\(\Big \downarrow \) 2009

\(\displaystyle b^2 \left (-\frac {\text {arctanh}(a+b x)}{\left (1-a^2\right ) b x}-\frac {2 a \text {arctanh}(a+b x) \log \left (\frac {2}{a+b x+1}\right )}{\left (1-a^2\right )^2}+\frac {2 a \text {arctanh}(a+b x) \log \left (\frac {2 b x}{(1-a) (a+b x+1)}\right )}{\left (1-a^2\right )^2}+\frac {a \operatorname {PolyLog}\left (2,1-\frac {2}{a+b x+1}\right )}{\left (1-a^2\right )^2}-\frac {a \operatorname {PolyLog}\left (2,1-\frac {2 b x}{(1-a) (a+b x+1)}\right )}{\left (1-a^2\right )^2}+\frac {\log (-b x)}{\left (1-a^2\right )^2}+\frac {\text {arctanh}(a+b x) \log \left (\frac {2}{-a-b x+1}\right )}{2 (1-a)^2}-\frac {\text {arctanh}(a+b x) \log \left (\frac {2}{a+b x+1}\right )}{2 (a+1)^2}+\frac {\operatorname {PolyLog}\left (2,-\frac {a+b x+1}{-a-b x+1}\right )}{4 (1-a)^2}+\frac {\operatorname {PolyLog}\left (2,1-\frac {2}{a+b x+1}\right )}{4 (a+1)^2}-\frac {\log (-a-b x+1)}{2 (1-a)^2 (a+1)}-\frac {\log (a+b x+1)}{2 (1-a) (a+1)^2}\right )-\frac {\text {arctanh}(a+b x)^2}{2 x^2}\)

Input:

Int[ArcTanh[a + b*x]^2/x^3,x]
 

Output:

-1/2*ArcTanh[a + b*x]^2/x^2 + b^2*(-(ArcTanh[a + b*x]/((1 - a^2)*b*x)) + L 
og[-(b*x)]/(1 - a^2)^2 + (ArcTanh[a + b*x]*Log[2/(1 - a - b*x)])/(2*(1 - a 
)^2) - Log[1 - a - b*x]/(2*(1 - a)^2*(1 + a)) - (ArcTanh[a + b*x]*Log[2/(1 
 + a + b*x)])/(2*(1 + a)^2) - (2*a*ArcTanh[a + b*x]*Log[2/(1 + a + b*x)])/ 
(1 - a^2)^2 + (2*a*ArcTanh[a + b*x]*Log[(2*b*x)/((1 - a)*(1 + a + b*x))])/ 
(1 - a^2)^2 - Log[1 + a + b*x]/(2*(1 - a)*(1 + a)^2) + PolyLog[2, -((1 + a 
 + b*x)/(1 - a - b*x))]/(4*(1 - a)^2) + PolyLog[2, 1 - 2/(1 + a + b*x)]/(4 
*(1 + a)^2) + (a*PolyLog[2, 1 - 2/(1 + a + b*x)])/(1 - a^2)^2 - (a*PolyLog 
[2, 1 - (2*b*x)/((1 - a)*(1 + a + b*x))])/(1 - a^2)^2)
 

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 2009
Int[u_, x_Symbol] :> Simp[IntSum[u, x], x] /; SumQ[u]
 

rule 6659
Int[((a_.) + ArcTanh[(c_) + (d_.)*(x_)]*(b_.))^(p_.)*((e_.) + (f_.)*(x_))^( 
m_), x_Symbol] :> Simp[(e + f*x)^(m + 1)*((a + b*ArcTanh[c + d*x])^p/(f*(m 
+ 1))), x] - Simp[b*d*(p/(f*(m + 1)))   Int[(e + f*x)^(m + 1)*((a + b*ArcTa 
nh[c + d*x])^(p - 1)/(1 - (c + d*x)^2)), x], x] /; FreeQ[{a, b, c, d, e, f} 
, x] && IGtQ[p, 0] && ILtQ[m, -1]
 

rule 6671
Int[((a_.) + ArcTanh[(c_) + (d_.)*(x_)]*(b_.))^(p_.)*((e_.) + (f_.)*(x_))^( 
m_.)*((A_.) + (B_.)*(x_) + (C_.)*(x_)^2)^(q_.), x_Symbol] :> Simp[1/d   Sub 
st[Int[((d*e - c*f)/d + f*(x/d))^m*(-C/d^2 + (C/d^2)*x^2)^q*(a + b*ArcTanh[ 
x])^p, x], x, c + d*x], x] /; FreeQ[{a, b, c, d, e, f, A, B, C, m, p, q}, x 
] && EqQ[B*(1 - c^2) + 2*A*c*d, 0] && EqQ[2*c*C - B*d, 0]
 

rule 7276
Int[(u_)/((a_) + (b_.)*(x_)^(n_)), x_Symbol] :> With[{v = RationalFunctionE 
xpand[u/(a + b*x^n), x]}, Int[v, x] /; SumQ[v]] /; FreeQ[{a, b}, x] && IGtQ 
[n, 0]
 

rule 7292
Int[u_, x_Symbol] :> With[{v = NormalizeIntegrand[u, x]}, Int[v, x] /; v =! 
= u]
 
Maple [A] (verified)

Time = 0.27 (sec) , antiderivative size = 385, normalized size of antiderivative = 1.04

method result size
parts \(-\frac {\operatorname {arctanh}\left (b x +a \right )^{2}}{2 x^{2}}+b^{2} \left (\frac {\operatorname {arctanh}\left (b x +a \right )}{\left (a -1\right ) \left (a +1\right ) b x}+\frac {2 \,\operatorname {arctanh}\left (b x +a \right ) a \ln \left (-b x \right )}{\left (a -1\right )^{2} \left (a +1\right )^{2}}+\frac {\operatorname {arctanh}\left (b x +a \right ) \ln \left (b x +a +1\right )}{2 \left (a +1\right )^{2}}-\frac {\operatorname {arctanh}\left (b x +a \right ) \ln \left (b x +a -1\right )}{2 \left (a -1\right )^{2}}-\frac {\frac {\ln \left (b x +a -1\right )^{2}}{4}-\frac {\operatorname {dilog}\left (\frac {b x}{2}+\frac {a}{2}+\frac {1}{2}\right )}{2}-\frac {\ln \left (b x +a -1\right ) \ln \left (\frac {b x}{2}+\frac {a}{2}+\frac {1}{2}\right )}{2}}{2 \left (a -1\right )^{2}}+\frac {-\frac {\ln \left (b x +a +1\right )^{2}}{4}+\frac {\left (\ln \left (b x +a +1\right )-\ln \left (\frac {b x}{2}+\frac {a}{2}+\frac {1}{2}\right )\right ) \ln \left (-\frac {b x}{2}-\frac {a}{2}+\frac {1}{2}\right )}{2}-\frac {\operatorname {dilog}\left (\frac {b x}{2}+\frac {a}{2}+\frac {1}{2}\right )}{2}}{2 \left (a +1\right )^{2}}-\frac {-\frac {\ln \left (-b x \right )}{\left (a -1\right ) \left (a +1\right )}-\frac {\ln \left (b x +a +1\right )}{2 a +2}+\frac {\ln \left (b x +a -1\right )}{2 a -2}}{\left (a -1\right ) \left (a +1\right )}+\frac {2 a \left (-\frac {\operatorname {dilog}\left (\frac {-b x -a -1}{-1-a}\right )}{2}-\frac {\ln \left (-b x \right ) \ln \left (\frac {-b x -a -1}{-1-a}\right )}{2}+\frac {\operatorname {dilog}\left (\frac {-b x -a +1}{1-a}\right )}{2}+\frac {\ln \left (-b x \right ) \ln \left (\frac {-b x -a +1}{1-a}\right )}{2}\right )}{\left (a -1\right )^{2} \left (a +1\right )^{2}}\right )\) \(385\)
derivativedivides \(b^{2} \left (-\frac {\operatorname {arctanh}\left (b x +a \right )^{2}}{2 b^{2} x^{2}}+\frac {\operatorname {arctanh}\left (b x +a \right )}{\left (a -1\right ) \left (a +1\right ) b x}+\frac {2 \,\operatorname {arctanh}\left (b x +a \right ) a \ln \left (-b x \right )}{\left (a -1\right )^{2} \left (a +1\right )^{2}}+\frac {\operatorname {arctanh}\left (b x +a \right ) \ln \left (b x +a +1\right )}{2 \left (a +1\right )^{2}}-\frac {\operatorname {arctanh}\left (b x +a \right ) \ln \left (b x +a -1\right )}{2 \left (a -1\right )^{2}}-\frac {\frac {\ln \left (b x +a -1\right )^{2}}{4}-\frac {\operatorname {dilog}\left (\frac {b x}{2}+\frac {a}{2}+\frac {1}{2}\right )}{2}-\frac {\ln \left (b x +a -1\right ) \ln \left (\frac {b x}{2}+\frac {a}{2}+\frac {1}{2}\right )}{2}}{2 \left (a -1\right )^{2}}+\frac {-\frac {\ln \left (b x +a +1\right )^{2}}{4}+\frac {\left (\ln \left (b x +a +1\right )-\ln \left (\frac {b x}{2}+\frac {a}{2}+\frac {1}{2}\right )\right ) \ln \left (-\frac {b x}{2}-\frac {a}{2}+\frac {1}{2}\right )}{2}-\frac {\operatorname {dilog}\left (\frac {b x}{2}+\frac {a}{2}+\frac {1}{2}\right )}{2}}{2 \left (a +1\right )^{2}}-\frac {-\frac {\ln \left (-b x \right )}{\left (a -1\right ) \left (a +1\right )}-\frac {\ln \left (b x +a +1\right )}{2 a +2}+\frac {\ln \left (b x +a -1\right )}{2 a -2}}{\left (a -1\right ) \left (a +1\right )}+\frac {2 a \left (-\frac {\operatorname {dilog}\left (\frac {-b x -a -1}{-1-a}\right )}{2}-\frac {\ln \left (-b x \right ) \ln \left (\frac {-b x -a -1}{-1-a}\right )}{2}+\frac {\operatorname {dilog}\left (\frac {-b x -a +1}{1-a}\right )}{2}+\frac {\ln \left (-b x \right ) \ln \left (\frac {-b x -a +1}{1-a}\right )}{2}\right )}{\left (a -1\right )^{2} \left (a +1\right )^{2}}\right )\) \(387\)
default \(b^{2} \left (-\frac {\operatorname {arctanh}\left (b x +a \right )^{2}}{2 b^{2} x^{2}}+\frac {\operatorname {arctanh}\left (b x +a \right )}{\left (a -1\right ) \left (a +1\right ) b x}+\frac {2 \,\operatorname {arctanh}\left (b x +a \right ) a \ln \left (-b x \right )}{\left (a -1\right )^{2} \left (a +1\right )^{2}}+\frac {\operatorname {arctanh}\left (b x +a \right ) \ln \left (b x +a +1\right )}{2 \left (a +1\right )^{2}}-\frac {\operatorname {arctanh}\left (b x +a \right ) \ln \left (b x +a -1\right )}{2 \left (a -1\right )^{2}}-\frac {\frac {\ln \left (b x +a -1\right )^{2}}{4}-\frac {\operatorname {dilog}\left (\frac {b x}{2}+\frac {a}{2}+\frac {1}{2}\right )}{2}-\frac {\ln \left (b x +a -1\right ) \ln \left (\frac {b x}{2}+\frac {a}{2}+\frac {1}{2}\right )}{2}}{2 \left (a -1\right )^{2}}+\frac {-\frac {\ln \left (b x +a +1\right )^{2}}{4}+\frac {\left (\ln \left (b x +a +1\right )-\ln \left (\frac {b x}{2}+\frac {a}{2}+\frac {1}{2}\right )\right ) \ln \left (-\frac {b x}{2}-\frac {a}{2}+\frac {1}{2}\right )}{2}-\frac {\operatorname {dilog}\left (\frac {b x}{2}+\frac {a}{2}+\frac {1}{2}\right )}{2}}{2 \left (a +1\right )^{2}}-\frac {-\frac {\ln \left (-b x \right )}{\left (a -1\right ) \left (a +1\right )}-\frac {\ln \left (b x +a +1\right )}{2 a +2}+\frac {\ln \left (b x +a -1\right )}{2 a -2}}{\left (a -1\right ) \left (a +1\right )}+\frac {2 a \left (-\frac {\operatorname {dilog}\left (\frac {-b x -a -1}{-1-a}\right )}{2}-\frac {\ln \left (-b x \right ) \ln \left (\frac {-b x -a -1}{-1-a}\right )}{2}+\frac {\operatorname {dilog}\left (\frac {-b x -a +1}{1-a}\right )}{2}+\frac {\ln \left (-b x \right ) \ln \left (\frac {-b x -a +1}{1-a}\right )}{2}\right )}{\left (a -1\right )^{2} \left (a +1\right )^{2}}\right )\) \(387\)

Input:

int(arctanh(b*x+a)^2/x^3,x,method=_RETURNVERBOSE)
 

Output:

-1/2*arctanh(b*x+a)^2/x^2+b^2*(arctanh(b*x+a)/(a-1)/(a+1)/b/x+2*arctanh(b* 
x+a)*a/(a-1)^2/(a+1)^2*ln(-b*x)+1/2*arctanh(b*x+a)/(a+1)^2*ln(b*x+a+1)-1/2 
*arctanh(b*x+a)/(a-1)^2*ln(b*x+a-1)-1/2/(a-1)^2*(1/4*ln(b*x+a-1)^2-1/2*dil 
og(1/2*b*x+1/2*a+1/2)-1/2*ln(b*x+a-1)*ln(1/2*b*x+1/2*a+1/2))+1/2/(a+1)^2*( 
-1/4*ln(b*x+a+1)^2+1/2*(ln(b*x+a+1)-ln(1/2*b*x+1/2*a+1/2))*ln(-1/2*b*x-1/2 
*a+1/2)-1/2*dilog(1/2*b*x+1/2*a+1/2))-1/(a-1)/(a+1)*(-1/(a-1)/(a+1)*ln(-b* 
x)-1/(2*a+2)*ln(b*x+a+1)+1/(2*a-2)*ln(b*x+a-1))+2*a/(a-1)^2/(a+1)^2*(-1/2* 
dilog((-b*x-a-1)/(-1-a))-1/2*ln(-b*x)*ln((-b*x-a-1)/(-1-a))+1/2*dilog(1/(1 
-a)*(-b*x-a+1))+1/2*ln(-b*x)*ln(1/(1-a)*(-b*x-a+1))))
 

Fricas [F]

\[ \int \frac {\text {arctanh}(a+b x)^2}{x^3} \, dx=\int { \frac {\operatorname {artanh}\left (b x + a\right )^{2}}{x^{3}} \,d x } \] Input:

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

Output:

integral(arctanh(b*x + a)^2/x^3, x)
 

Sympy [F]

\[ \int \frac {\text {arctanh}(a+b x)^2}{x^3} \, dx=\int \frac {\operatorname {atanh}^{2}{\left (a + b x \right )}}{x^{3}}\, dx \] Input:

integrate(atanh(b*x+a)**2/x**3,x)
 

Output:

Integral(atanh(a + b*x)**2/x**3, x)
 

Maxima [A] (verification not implemented)

Time = 0.05 (sec) , antiderivative size = 360, normalized size of antiderivative = 0.97 \[ \int \frac {\text {arctanh}(a+b x)^2}{x^3} \, dx=\frac {1}{8} \, {\left (\frac {8 \, {\left (\log \left (b x + a - 1\right ) \log \left (\frac {1}{2} \, b x + \frac {1}{2} \, a + \frac {1}{2}\right ) + {\rm Li}_2\left (-\frac {1}{2} \, b x - \frac {1}{2} \, a + \frac {1}{2}\right )\right )} a}{a^{4} - 2 \, a^{2} + 1} - \frac {8 \, {\left (\log \left (\frac {b x}{a + 1} + 1\right ) \log \left (x\right ) + {\rm Li}_2\left (-\frac {b x}{a + 1}\right )\right )} a}{a^{4} - 2 \, a^{2} + 1} + \frac {8 \, {\left (\log \left (\frac {b x}{a - 1} + 1\right ) \log \left (x\right ) + {\rm Li}_2\left (-\frac {b x}{a - 1}\right )\right )} a}{a^{4} - 2 \, a^{2} + 1} - \frac {{\left (a^{2} - 2 \, a + 1\right )} \log \left (b x + a + 1\right )^{2} - 2 \, {\left (a^{2} - 2 \, a + 1\right )} \log \left (b x + a + 1\right ) \log \left (b x + a - 1\right ) + {\left (a^{2} + 2 \, a + 1\right )} \log \left (b x + a - 1\right )^{2}}{a^{4} - 2 \, a^{2} + 1} + \frac {4 \, \log \left (b x + a + 1\right )}{a^{3} + a^{2} - a - 1} - \frac {4 \, \log \left (b x + a - 1\right )}{a^{3} - a^{2} - a + 1} + \frac {8 \, \log \left (x\right )}{a^{4} - 2 \, a^{2} + 1}\right )} b^{2} + \frac {1}{2} \, {\left (\frac {4 \, a b \log \left (x\right )}{a^{4} - 2 \, a^{2} + 1} + \frac {b \log \left (b x + a + 1\right )}{a^{2} + 2 \, a + 1} - \frac {b \log \left (b x + a - 1\right )}{a^{2} - 2 \, a + 1} + \frac {2}{{\left (a^{2} - 1\right )} x}\right )} b \operatorname {artanh}\left (b x + a\right ) - \frac {\operatorname {artanh}\left (b x + a\right )^{2}}{2 \, x^{2}} \] Input:

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

Output:

1/8*(8*(log(b*x + a - 1)*log(1/2*b*x + 1/2*a + 1/2) + dilog(-1/2*b*x - 1/2 
*a + 1/2))*a/(a^4 - 2*a^2 + 1) - 8*(log(b*x/(a + 1) + 1)*log(x) + dilog(-b 
*x/(a + 1)))*a/(a^4 - 2*a^2 + 1) + 8*(log(b*x/(a - 1) + 1)*log(x) + dilog( 
-b*x/(a - 1)))*a/(a^4 - 2*a^2 + 1) - ((a^2 - 2*a + 1)*log(b*x + a + 1)^2 - 
 2*(a^2 - 2*a + 1)*log(b*x + a + 1)*log(b*x + a - 1) + (a^2 + 2*a + 1)*log 
(b*x + a - 1)^2)/(a^4 - 2*a^2 + 1) + 4*log(b*x + a + 1)/(a^3 + a^2 - a - 1 
) - 4*log(b*x + a - 1)/(a^3 - a^2 - a + 1) + 8*log(x)/(a^4 - 2*a^2 + 1))*b 
^2 + 1/2*(4*a*b*log(x)/(a^4 - 2*a^2 + 1) + b*log(b*x + a + 1)/(a^2 + 2*a + 
 1) - b*log(b*x + a - 1)/(a^2 - 2*a + 1) + 2/((a^2 - 1)*x))*b*arctanh(b*x 
+ a) - 1/2*arctanh(b*x + a)^2/x^2
 

Giac [F]

\[ \int \frac {\text {arctanh}(a+b x)^2}{x^3} \, dx=\int { \frac {\operatorname {artanh}\left (b x + a\right )^{2}}{x^{3}} \,d x } \] Input:

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

Output:

integrate(arctanh(b*x + a)^2/x^3, x)
 

Mupad [F(-1)]

Timed out. \[ \int \frac {\text {arctanh}(a+b x)^2}{x^3} \, dx=\int \frac {{\mathrm {atanh}\left (a+b\,x\right )}^2}{x^3} \,d x \] Input:

int(atanh(a + b*x)^2/x^3,x)
 

Output:

int(atanh(a + b*x)^2/x^3, x)
 

Reduce [F]

\[ \int \frac {\text {arctanh}(a+b x)^2}{x^3} \, dx=\frac {5 \mathit {atanh} \left (b x +a \right )^{2} a^{4}-\mathit {atanh} \left (b x +a \right )^{2} a^{2}-4 \mathit {atanh} \left (b x +a \right ) a^{3}+2 \mathit {atanh} \left (b x +a \right ) a -3 \mathit {atanh} \left (b x +a \right )^{2} a^{6}+2 \mathit {atanh} \left (b x +a \right ) a^{5}-\mathit {atanh} \left (b x +a \right )^{2}-2 \mathit {atanh} \left (b x +a \right ) b x +2 a b x +\mathit {atanh} \left (b x +a \right )^{2} a^{4} b^{2} x^{2}+\mathit {atanh} \left (b x +a \right )^{2} b^{2} x^{2}+4 \mathit {atanh} \left (b x +a \right ) a^{2} b x -2 \mathit {atanh} \left (b x +a \right )^{2} a^{2} b^{2} x^{2}-2 \mathit {atanh} \left (b x +a \right ) a^{4} b x -4 \mathit {atanh} \left (b x +a \right ) a^{3} b^{2} x^{2}+6 \mathit {atanh} \left (b x +a \right ) a^{2} b^{2} x^{2}+6 \,\mathrm {log}\left (b x +a -1\right ) a^{2} b^{2} x^{2}-6 \,\mathrm {log}\left (x \right ) a^{2} b^{2} x^{2}+12 \left (\int \frac {\mathit {atanh} \left (b x +a \right )}{3 a^{2} b^{2} x^{5}+6 a^{3} b \,x^{4}+3 a^{4} x^{3}+b^{2} x^{5}+2 a b \,x^{4}-2 a^{2} x^{3}-x^{3}}d x \right ) a^{9} x^{2}-32 \left (\int \frac {\mathit {atanh} \left (b x +a \right )}{3 a^{2} b^{2} x^{5}+6 a^{3} b \,x^{4}+3 a^{4} x^{3}+b^{2} x^{5}+2 a b \,x^{4}-2 a^{2} x^{3}-x^{3}}d x \right ) a^{7} x^{2}+24 \left (\int \frac {\mathit {atanh} \left (b x +a \right )}{3 a^{2} b^{2} x^{5}+6 a^{3} b \,x^{4}+3 a^{4} x^{3}+b^{2} x^{5}+2 a b \,x^{4}-2 a^{2} x^{3}-x^{3}}d x \right ) a^{5} x^{2}-4 \left (\int \frac {\mathit {atanh} \left (b x +a \right )}{3 a^{2} b^{2} x^{5}+6 a^{3} b \,x^{4}+3 a^{4} x^{3}+b^{2} x^{5}+2 a b \,x^{4}-2 a^{2} x^{3}-x^{3}}d x \right ) a \,x^{2}-2 \mathit {atanh} \left (b x +a \right ) b^{2} x^{2}-2 \,\mathrm {log}\left (b x +a -1\right ) b^{2} x^{2}+2 \,\mathrm {log}\left (x \right ) b^{2} x^{2}-2 a^{3} b x}{2 x^{2} \left (3 a^{6}-5 a^{4}+a^{2}+1\right )} \] Input:

int(atanh(b*x+a)^2/x^3,x)
 

Output:

( - 3*atanh(a + b*x)**2*a**6 + atanh(a + b*x)**2*a**4*b**2*x**2 + 5*atanh( 
a + b*x)**2*a**4 - 2*atanh(a + b*x)**2*a**2*b**2*x**2 - atanh(a + b*x)**2* 
a**2 + atanh(a + b*x)**2*b**2*x**2 - atanh(a + b*x)**2 + 2*atanh(a + b*x)* 
a**5 - 2*atanh(a + b*x)*a**4*b*x - 4*atanh(a + b*x)*a**3*b**2*x**2 - 4*ata 
nh(a + b*x)*a**3 + 6*atanh(a + b*x)*a**2*b**2*x**2 + 4*atanh(a + b*x)*a**2 
*b*x + 2*atanh(a + b*x)*a - 2*atanh(a + b*x)*b**2*x**2 - 2*atanh(a + b*x)* 
b*x + 12*int(atanh(a + b*x)/(3*a**4*x**3 + 6*a**3*b*x**4 + 3*a**2*b**2*x** 
5 - 2*a**2*x**3 + 2*a*b*x**4 + b**2*x**5 - x**3),x)*a**9*x**2 - 32*int(ata 
nh(a + b*x)/(3*a**4*x**3 + 6*a**3*b*x**4 + 3*a**2*b**2*x**5 - 2*a**2*x**3 
+ 2*a*b*x**4 + b**2*x**5 - x**3),x)*a**7*x**2 + 24*int(atanh(a + b*x)/(3*a 
**4*x**3 + 6*a**3*b*x**4 + 3*a**2*b**2*x**5 - 2*a**2*x**3 + 2*a*b*x**4 + b 
**2*x**5 - x**3),x)*a**5*x**2 - 4*int(atanh(a + b*x)/(3*a**4*x**3 + 6*a**3 
*b*x**4 + 3*a**2*b**2*x**5 - 2*a**2*x**3 + 2*a*b*x**4 + b**2*x**5 - x**3), 
x)*a*x**2 + 6*log(a + b*x - 1)*a**2*b**2*x**2 - 2*log(a + b*x - 1)*b**2*x* 
*2 - 6*log(x)*a**2*b**2*x**2 + 2*log(x)*b**2*x**2 - 2*a**3*b*x + 2*a*b*x)/ 
(2*x**2*(3*a**6 - 5*a**4 + a**2 + 1))