\(\int (1-a^2 x^2)^3 \text {arctanh}(a x) \, dx\) [224]

Optimal result
Mathematica [A] (verified)
Rubi [A] (verified)
Maple [A] (verified)
Fricas [A] (verification not implemented)
Sympy [A] (verification not implemented)
Maxima [A] (verification not implemented)
Giac [B] (verification not implemented)
Mupad [B] (verification not implemented)
Reduce [B] (verification not implemented)

Optimal result

Integrand size = 17, antiderivative size = 144 \[ \int \left (1-a^2 x^2\right )^3 \text {arctanh}(a x) \, dx=\frac {4 \left (1-a^2 x^2\right )}{35 a}+\frac {3 \left (1-a^2 x^2\right )^2}{70 a}+\frac {\left (1-a^2 x^2\right )^3}{42 a}+\frac {16}{35} x \text {arctanh}(a x)+\frac {8}{35} x \left (1-a^2 x^2\right ) \text {arctanh}(a x)+\frac {6}{35} x \left (1-a^2 x^2\right )^2 \text {arctanh}(a x)+\frac {1}{7} x \left (1-a^2 x^2\right )^3 \text {arctanh}(a x)+\frac {8 \log \left (1-a^2 x^2\right )}{35 a} \] Output:

4/35*(-a^2*x^2+1)/a+3/70*(-a^2*x^2+1)^2/a+1/42*(-a^2*x^2+1)^3/a+16/35*x*ar 
ctanh(a*x)+8/35*x*(-a^2*x^2+1)*arctanh(a*x)+6/35*x*(-a^2*x^2+1)^2*arctanh( 
a*x)+1/7*x*(-a^2*x^2+1)^3*arctanh(a*x)+8/35*ln(-a^2*x^2+1)/a
 

Mathematica [A] (verified)

Time = 0.03 (sec) , antiderivative size = 79, normalized size of antiderivative = 0.55 \[ \int \left (1-a^2 x^2\right )^3 \text {arctanh}(a x) \, dx=\frac {-57 a^2 x^2+24 a^4 x^4-5 a^6 x^6-6 a x \left (-35+35 a^2 x^2-21 a^4 x^4+5 a^6 x^6\right ) \text {arctanh}(a x)+48 \log \left (1-a^2 x^2\right )}{210 a} \] Input:

Integrate[(1 - a^2*x^2)^3*ArcTanh[a*x],x]
 

Output:

(-57*a^2*x^2 + 24*a^4*x^4 - 5*a^6*x^6 - 6*a*x*(-35 + 35*a^2*x^2 - 21*a^4*x 
^4 + 5*a^6*x^6)*ArcTanh[a*x] + 48*Log[1 - a^2*x^2])/(210*a)
 

Rubi [A] (verified)

Time = 0.52 (sec) , antiderivative size = 156, normalized size of antiderivative = 1.08, number of steps used = 5, number of rules used = 5, \(\frac {\text {number of rules}}{\text {integrand size}}\) = 0.294, Rules used = {6504, 6504, 6504, 6436, 240}

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 \left (1-a^2 x^2\right )^3 \text {arctanh}(a x) \, dx\)

\(\Big \downarrow \) 6504

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

\(\Big \downarrow \) 6504

\(\displaystyle \frac {6}{7} \left (\frac {4}{5} \int \left (1-a^2 x^2\right ) \text {arctanh}(a x)dx+\frac {1}{5} x \left (1-a^2 x^2\right )^2 \text {arctanh}(a x)+\frac {\left (1-a^2 x^2\right )^2}{20 a}\right )+\frac {1}{7} x \left (1-a^2 x^2\right )^3 \text {arctanh}(a x)+\frac {\left (1-a^2 x^2\right )^3}{42 a}\)

\(\Big \downarrow \) 6504

\(\displaystyle \frac {6}{7} \left (\frac {4}{5} \left (\frac {2}{3} \int \text {arctanh}(a x)dx+\frac {1}{3} x \left (1-a^2 x^2\right ) \text {arctanh}(a x)+\frac {1-a^2 x^2}{6 a}\right )+\frac {1}{5} x \left (1-a^2 x^2\right )^2 \text {arctanh}(a x)+\frac {\left (1-a^2 x^2\right )^2}{20 a}\right )+\frac {1}{7} x \left (1-a^2 x^2\right )^3 \text {arctanh}(a x)+\frac {\left (1-a^2 x^2\right )^3}{42 a}\)

\(\Big \downarrow \) 6436

\(\displaystyle \frac {6}{7} \left (\frac {4}{5} \left (\frac {2}{3} \left (x \text {arctanh}(a x)-a \int \frac {x}{1-a^2 x^2}dx\right )+\frac {1}{3} x \left (1-a^2 x^2\right ) \text {arctanh}(a x)+\frac {1-a^2 x^2}{6 a}\right )+\frac {1}{5} x \left (1-a^2 x^2\right )^2 \text {arctanh}(a x)+\frac {\left (1-a^2 x^2\right )^2}{20 a}\right )+\frac {1}{7} x \left (1-a^2 x^2\right )^3 \text {arctanh}(a x)+\frac {\left (1-a^2 x^2\right )^3}{42 a}\)

\(\Big \downarrow \) 240

\(\displaystyle \frac {1}{7} x \left (1-a^2 x^2\right )^3 \text {arctanh}(a x)+\frac {6}{7} \left (\frac {1}{5} x \left (1-a^2 x^2\right )^2 \text {arctanh}(a x)+\frac {4}{5} \left (\frac {1}{3} x \left (1-a^2 x^2\right ) \text {arctanh}(a x)+\frac {2}{3} \left (\frac {\log \left (1-a^2 x^2\right )}{2 a}+x \text {arctanh}(a x)\right )+\frac {1-a^2 x^2}{6 a}\right )+\frac {\left (1-a^2 x^2\right )^2}{20 a}\right )+\frac {\left (1-a^2 x^2\right )^3}{42 a}\)

Input:

Int[(1 - a^2*x^2)^3*ArcTanh[a*x],x]
 

Output:

(1 - a^2*x^2)^3/(42*a) + (x*(1 - a^2*x^2)^3*ArcTanh[a*x])/7 + (6*((1 - a^2 
*x^2)^2/(20*a) + (x*(1 - a^2*x^2)^2*ArcTanh[a*x])/5 + (4*((1 - a^2*x^2)/(6 
*a) + (x*(1 - a^2*x^2)*ArcTanh[a*x])/3 + (2*(x*ArcTanh[a*x] + Log[1 - a^2* 
x^2]/(2*a)))/3))/5))/7
 

Defintions of rubi rules used

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

rule 6436
Int[((a_.) + ArcTanh[(c_.)*(x_)^(n_.)]*(b_.))^(p_.), x_Symbol] :> Simp[x*(a 
 + b*ArcTanh[c*x^n])^p, x] - Simp[b*c*n*p   Int[x^n*((a + b*ArcTanh[c*x^n]) 
^(p - 1)/(1 - c^2*x^(2*n))), x], x] /; FreeQ[{a, b, c, n}, x] && IGtQ[p, 0] 
 && (EqQ[n, 1] || EqQ[p, 1])
 

rule 6504
Int[((a_.) + ArcTanh[(c_.)*(x_)]*(b_.))*((d_) + (e_.)*(x_)^2)^(q_.), x_Symb 
ol] :> Simp[b*((d + e*x^2)^q/(2*c*q*(2*q + 1))), x] + (Simp[x*(d + e*x^2)^q 
*((a + b*ArcTanh[c*x])/(2*q + 1)), x] + Simp[2*d*(q/(2*q + 1))   Int[(d + e 
*x^2)^(q - 1)*(a + b*ArcTanh[c*x]), x], x]) /; FreeQ[{a, b, c, d, e}, x] && 
 EqQ[c^2*d + e, 0] && GtQ[q, 0]
 
Maple [A] (verified)

Time = 0.35 (sec) , antiderivative size = 84, normalized size of antiderivative = 0.58

method result size
parts \(-\frac {\operatorname {arctanh}\left (a x \right ) a^{6} x^{7}}{7}+\frac {3 \,\operatorname {arctanh}\left (a x \right ) a^{4} x^{5}}{5}-x^{3} a^{2} \operatorname {arctanh}\left (a x \right )+x \,\operatorname {arctanh}\left (a x \right )-\frac {a \left (\frac {5 a^{4} x^{6}}{6}-4 a^{2} x^{4}+\frac {19 x^{2}}{2}-\frac {8 \ln \left (a^{2} x^{2}-1\right )}{a^{2}}\right )}{35}\) \(84\)
derivativedivides \(\frac {-\frac {\operatorname {arctanh}\left (a x \right ) a^{7} x^{7}}{7}+\frac {3 \,\operatorname {arctanh}\left (a x \right ) a^{5} x^{5}}{5}-a^{3} x^{3} \operatorname {arctanh}\left (a x \right )+a x \,\operatorname {arctanh}\left (a x \right )-\frac {a^{6} x^{6}}{42}+\frac {4 a^{4} x^{4}}{35}-\frac {19 a^{2} x^{2}}{70}+\frac {8 \ln \left (a x -1\right )}{35}+\frac {8 \ln \left (a x +1\right )}{35}}{a}\) \(89\)
default \(\frac {-\frac {\operatorname {arctanh}\left (a x \right ) a^{7} x^{7}}{7}+\frac {3 \,\operatorname {arctanh}\left (a x \right ) a^{5} x^{5}}{5}-a^{3} x^{3} \operatorname {arctanh}\left (a x \right )+a x \,\operatorname {arctanh}\left (a x \right )-\frac {a^{6} x^{6}}{42}+\frac {4 a^{4} x^{4}}{35}-\frac {19 a^{2} x^{2}}{70}+\frac {8 \ln \left (a x -1\right )}{35}+\frac {8 \ln \left (a x +1\right )}{35}}{a}\) \(89\)
parallelrisch \(-\frac {30 \,\operatorname {arctanh}\left (a x \right ) a^{7} x^{7}+5 a^{6} x^{6}-126 \,\operatorname {arctanh}\left (a x \right ) a^{5} x^{5}-24 a^{4} x^{4}+210 a^{3} x^{3} \operatorname {arctanh}\left (a x \right )+57 a^{2} x^{2}-210 a x \,\operatorname {arctanh}\left (a x \right )-96 \ln \left (a x -1\right )-96 \,\operatorname {arctanh}\left (a x \right )}{210 a}\) \(89\)
risch \(\left (-\frac {1}{14} a^{6} x^{7}+\frac {3}{10} a^{4} x^{5}-\frac {1}{2} a^{2} x^{3}+\frac {1}{2} x \right ) \ln \left (a x +1\right )+\frac {a^{6} x^{7} \ln \left (-a x +1\right )}{14}-\frac {a^{5} x^{6}}{42}-\frac {3 a^{4} x^{5} \ln \left (-a x +1\right )}{10}+\frac {4 a^{3} x^{4}}{35}+\frac {a^{2} x^{3} \ln \left (-a x +1\right )}{2}-\frac {19 a \,x^{2}}{70}-\frac {x \ln \left (-a x +1\right )}{2}+\frac {8 \ln \left (a^{2} x^{2}-1\right )}{35 a}\) \(129\)
meijerg \(-\frac {\frac {2 a^{2} x^{2} \left (\ln \left (1-\sqrt {a^{2} x^{2}}\right )-\ln \left (1+\sqrt {a^{2} x^{2}}\right )\right )}{\sqrt {a^{2} x^{2}}}-2 \ln \left (-a^{2} x^{2}+1\right )}{4 a}-\frac {\frac {x^{2} a^{2} \left (4 a^{4} x^{4}+6 a^{2} x^{2}+12\right )}{42}-\frac {2 x^{8} a^{8} \left (\ln \left (1-\sqrt {a^{2} x^{2}}\right )-\ln \left (1+\sqrt {a^{2} x^{2}}\right )\right )}{7 \sqrt {a^{2} x^{2}}}+\frac {2 \ln \left (-a^{2} x^{2}+1\right )}{7}}{4 a}-\frac {3 \left (-\frac {a^{2} x^{2} \left (3 a^{2} x^{2}+6\right )}{15}+\frac {2 a^{6} x^{6} \left (\ln \left (1-\sqrt {a^{2} x^{2}}\right )-\ln \left (1+\sqrt {a^{2} x^{2}}\right )\right )}{5 \sqrt {a^{2} x^{2}}}-\frac {2 \ln \left (-a^{2} x^{2}+1\right )}{5}\right )}{4 a}-\frac {3 \left (\frac {2 a^{2} x^{2}}{3}-\frac {2 a^{4} x^{4} \left (\ln \left (1-\sqrt {a^{2} x^{2}}\right )-\ln \left (1+\sqrt {a^{2} x^{2}}\right )\right )}{3 \sqrt {a^{2} x^{2}}}+\frac {2 \ln \left (-a^{2} x^{2}+1\right )}{3}\right )}{4 a}\) \(314\)

Input:

int((-a^2*x^2+1)^3*arctanh(a*x),x,method=_RETURNVERBOSE)
 

Output:

-1/7*arctanh(a*x)*a^6*x^7+3/5*arctanh(a*x)*a^4*x^5-x^3*a^2*arctanh(a*x)+x* 
arctanh(a*x)-1/35*a*(5/6*a^4*x^6-4*a^2*x^4+19/2*x^2-8/a^2*ln(a^2*x^2-1))
 

Fricas [A] (verification not implemented)

Time = 0.08 (sec) , antiderivative size = 88, normalized size of antiderivative = 0.61 \[ \int \left (1-a^2 x^2\right )^3 \text {arctanh}(a x) \, dx=-\frac {5 \, a^{6} x^{6} - 24 \, a^{4} x^{4} + 57 \, a^{2} x^{2} + 3 \, {\left (5 \, a^{7} x^{7} - 21 \, a^{5} x^{5} + 35 \, a^{3} x^{3} - 35 \, a x\right )} \log \left (-\frac {a x + 1}{a x - 1}\right ) - 48 \, \log \left (a^{2} x^{2} - 1\right )}{210 \, a} \] Input:

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

Output:

-1/210*(5*a^6*x^6 - 24*a^4*x^4 + 57*a^2*x^2 + 3*(5*a^7*x^7 - 21*a^5*x^5 + 
35*a^3*x^3 - 35*a*x)*log(-(a*x + 1)/(a*x - 1)) - 48*log(a^2*x^2 - 1))/a
 

Sympy [A] (verification not implemented)

Time = 0.46 (sec) , antiderivative size = 97, normalized size of antiderivative = 0.67 \[ \int \left (1-a^2 x^2\right )^3 \text {arctanh}(a x) \, dx=\begin {cases} - \frac {a^{6} x^{7} \operatorname {atanh}{\left (a x \right )}}{7} - \frac {a^{5} x^{6}}{42} + \frac {3 a^{4} x^{5} \operatorname {atanh}{\left (a x \right )}}{5} + \frac {4 a^{3} x^{4}}{35} - a^{2} x^{3} \operatorname {atanh}{\left (a x \right )} - \frac {19 a x^{2}}{70} + x \operatorname {atanh}{\left (a x \right )} + \frac {16 \log {\left (x - \frac {1}{a} \right )}}{35 a} + \frac {16 \operatorname {atanh}{\left (a x \right )}}{35 a} & \text {for}\: a \neq 0 \\0 & \text {otherwise} \end {cases} \] Input:

integrate((-a**2*x**2+1)**3*atanh(a*x),x)
 

Output:

Piecewise((-a**6*x**7*atanh(a*x)/7 - a**5*x**6/42 + 3*a**4*x**5*atanh(a*x) 
/5 + 4*a**3*x**4/35 - a**2*x**3*atanh(a*x) - 19*a*x**2/70 + x*atanh(a*x) + 
 16*log(x - 1/a)/(35*a) + 16*atanh(a*x)/(35*a), Ne(a, 0)), (0, True))
 

Maxima [A] (verification not implemented)

Time = 0.03 (sec) , antiderivative size = 82, normalized size of antiderivative = 0.57 \[ \int \left (1-a^2 x^2\right )^3 \text {arctanh}(a x) \, dx=-\frac {1}{210} \, {\left (5 \, a^{4} x^{6} - 24 \, a^{2} x^{4} + 57 \, x^{2} - \frac {48 \, \log \left (a x + 1\right )}{a^{2}} - \frac {48 \, \log \left (a x - 1\right )}{a^{2}}\right )} a - \frac {1}{35} \, {\left (5 \, a^{6} x^{7} - 21 \, a^{4} x^{5} + 35 \, a^{2} x^{3} - 35 \, x\right )} \operatorname {artanh}\left (a x\right ) \] Input:

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

Output:

-1/210*(5*a^4*x^6 - 24*a^2*x^4 + 57*x^2 - 48*log(a*x + 1)/a^2 - 48*log(a*x 
 - 1)/a^2)*a - 1/35*(5*a^6*x^7 - 21*a^4*x^5 + 35*a^2*x^3 - 35*x)*arctanh(a 
*x)
 

Giac [B] (verification not implemented)

Leaf count of result is larger than twice the leaf count of optimal. 303 vs. \(2 (122) = 244\).

Time = 0.14 (sec) , antiderivative size = 303, normalized size of antiderivative = 2.10 \[ \int \left (1-a^2 x^2\right )^3 \text {arctanh}(a x) \, dx=\frac {8}{105} \, a {\left (\frac {6 \, \log \left (\frac {{\left | -a x - 1 \right |}}{{\left | a x - 1 \right |}}\right )}{a^{2}} - \frac {6 \, \log \left ({\left | -\frac {a x + 1}{a x - 1} + 1 \right |}\right )}{a^{2}} - \frac {\frac {6 \, {\left (a x + 1\right )}^{5}}{{\left (a x - 1\right )}^{5}} - \frac {33 \, {\left (a x + 1\right )}^{4}}{{\left (a x - 1\right )}^{4}} + \frac {74 \, {\left (a x + 1\right )}^{3}}{{\left (a x - 1\right )}^{3}} - \frac {33 \, {\left (a x + 1\right )}^{2}}{{\left (a x - 1\right )}^{2}} + \frac {6 \, {\left (a x + 1\right )}}{a x - 1}}{a^{2} {\left (\frac {a x + 1}{a x - 1} - 1\right )}^{6}} - \frac {6 \, {\left (\frac {35 \, {\left (a x + 1\right )}^{3}}{{\left (a x - 1\right )}^{3}} - \frac {21 \, {\left (a x + 1\right )}^{2}}{{\left (a x - 1\right )}^{2}} + \frac {7 \, {\left (a x + 1\right )}}{a x - 1} - 1\right )} \log \left (-\frac {\frac {a {\left (\frac {a x + 1}{a x - 1} + 1\right )}}{\frac {{\left (a x + 1\right )} a}{a x - 1} - a} + 1}{\frac {a {\left (\frac {a x + 1}{a x - 1} + 1\right )}}{\frac {{\left (a x + 1\right )} a}{a x - 1} - a} - 1}\right )}{a^{2} {\left (\frac {a x + 1}{a x - 1} - 1\right )}^{7}}\right )} \] Input:

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

Output:

8/105*a*(6*log(abs(-a*x - 1)/abs(a*x - 1))/a^2 - 6*log(abs(-(a*x + 1)/(a*x 
 - 1) + 1))/a^2 - (6*(a*x + 1)^5/(a*x - 1)^5 - 33*(a*x + 1)^4/(a*x - 1)^4 
+ 74*(a*x + 1)^3/(a*x - 1)^3 - 33*(a*x + 1)^2/(a*x - 1)^2 + 6*(a*x + 1)/(a 
*x - 1))/(a^2*((a*x + 1)/(a*x - 1) - 1)^6) - 6*(35*(a*x + 1)^3/(a*x - 1)^3 
 - 21*(a*x + 1)^2/(a*x - 1)^2 + 7*(a*x + 1)/(a*x - 1) - 1)*log(-(a*((a*x + 
 1)/(a*x - 1) + 1)/((a*x + 1)*a/(a*x - 1) - a) + 1)/(a*((a*x + 1)/(a*x - 1 
) + 1)/((a*x + 1)*a/(a*x - 1) - a) - 1))/(a^2*((a*x + 1)/(a*x - 1) - 1)^7) 
)
 

Mupad [B] (verification not implemented)

Time = 3.58 (sec) , antiderivative size = 80, normalized size of antiderivative = 0.56 \[ \int \left (1-a^2 x^2\right )^3 \text {arctanh}(a x) \, dx=x\,\mathrm {atanh}\left (a\,x\right )-\frac {19\,a\,x^2}{70}+\frac {8\,\ln \left (a^2\,x^2-1\right )}{35\,a}+\frac {4\,a^3\,x^4}{35}-\frac {a^5\,x^6}{42}-a^2\,x^3\,\mathrm {atanh}\left (a\,x\right )+\frac {3\,a^4\,x^5\,\mathrm {atanh}\left (a\,x\right )}{5}-\frac {a^6\,x^7\,\mathrm {atanh}\left (a\,x\right )}{7} \] Input:

int(-atanh(a*x)*(a^2*x^2 - 1)^3,x)
 

Output:

x*atanh(a*x) - (19*a*x^2)/70 + (8*log(a^2*x^2 - 1))/(35*a) + (4*a^3*x^4)/3 
5 - (a^5*x^6)/42 - a^2*x^3*atanh(a*x) + (3*a^4*x^5*atanh(a*x))/5 - (a^6*x^ 
7*atanh(a*x))/7
 

Reduce [B] (verification not implemented)

Time = 0.17 (sec) , antiderivative size = 92, normalized size of antiderivative = 0.64 \[ \int \left (1-a^2 x^2\right )^3 \text {arctanh}(a x) \, dx=\frac {-30 \mathit {atanh} \left (a x \right ) a^{7} x^{7}+126 \mathit {atanh} \left (a x \right ) a^{5} x^{5}-210 \mathit {atanh} \left (a x \right ) a^{3} x^{3}+210 \mathit {atanh} \left (a x \right ) a x +96 \mathit {atanh} \left (a x \right )+96 \,\mathrm {log}\left (a^{2} x -a \right )-5 a^{6} x^{6}+24 a^{4} x^{4}-57 a^{2} x^{2}}{210 a} \] Input:

int((-a^2*x^2+1)^3*atanh(a*x),x)
 

Output:

( - 30*atanh(a*x)*a**7*x**7 + 126*atanh(a*x)*a**5*x**5 - 210*atanh(a*x)*a* 
*3*x**3 + 210*atanh(a*x)*a*x + 96*atanh(a*x) + 96*log(a**2*x - a) - 5*a**6 
*x**6 + 24*a**4*x**4 - 57*a**2*x**2)/(210*a)