\(\int \frac {e^{-2 \text {arctanh}(a x)}}{(c-\frac {c}{a^2 x^2})^4} \, dx\) [718]

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 [A] (verification not implemented)
Mupad [B] (verification not implemented)
Reduce [B] (verification not implemented)

Optimal result

Integrand size = 22, antiderivative size = 144 \[ \int \frac {e^{-2 \text {arctanh}(a x)}}{\left (c-\frac {c}{a^2 x^2}\right )^4} \, dx=-\frac {x}{c^4}+\frac {1}{64 a c^4 (1-a x)^2}-\frac {11}{64 a c^4 (1-a x)}-\frac {1}{32 a c^4 (1+a x)^4}+\frac {13}{48 a c^4 (1+a x)^3}-\frac {35}{32 a c^4 (1+a x)^2}+\frac {99}{32 a c^4 (1+a x)}-\frac {47 \log (1-a x)}{128 a c^4}+\frac {303 \log (1+a x)}{128 a c^4} \] Output:

-x/c^4+1/64/a/c^4/(-a*x+1)^2-11/64/a/c^4/(-a*x+1)-1/32/a/c^4/(a*x+1)^4+13/ 
48/a/c^4/(a*x+1)^3-35/32/a/c^4/(a*x+1)^2+99/32/a/c^4/(a*x+1)-47/128*ln(-a* 
x+1)/a/c^4+303/128*ln(a*x+1)/a/c^4
                                                                                    
                                                                                    
 

Mathematica [A] (verified)

Time = 0.11 (sec) , antiderivative size = 121, normalized size of antiderivative = 0.84 \[ \int \frac {e^{-2 \text {arctanh}(a x)}}{\left (c-\frac {c}{a^2 x^2}\right )^4} \, dx=\frac {800+550 a x-2516 a^2 x^2-1732 a^3 x^3+2508 a^4 x^4+1638 a^5 x^5-768 a^6 x^6-384 a^7 x^7-141 (-1+a x)^2 (1+a x)^4 \log (1-a x)+909 (-1+a x)^2 (1+a x)^4 \log (1+a x)}{384 a (-1+a x)^2 (c+a c x)^4} \] Input:

Integrate[1/(E^(2*ArcTanh[a*x])*(c - c/(a^2*x^2))^4),x]
 

Output:

(800 + 550*a*x - 2516*a^2*x^2 - 1732*a^3*x^3 + 2508*a^4*x^4 + 1638*a^5*x^5 
 - 768*a^6*x^6 - 384*a^7*x^7 - 141*(-1 + a*x)^2*(1 + a*x)^4*Log[1 - a*x] + 
 909*(-1 + a*x)^2*(1 + a*x)^4*Log[1 + a*x])/(384*a*(-1 + a*x)^2*(c + a*c*x 
)^4)
 

Rubi [A] (verified)

Time = 0.77 (sec) , antiderivative size = 127, normalized size of antiderivative = 0.88, number of steps used = 4, number of rules used = 4, \(\frac {\text {number of rules}}{\text {integrand size}}\) = 0.182, Rules used = {6707, 6700, 99, 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 {e^{-2 \text {arctanh}(a x)}}{\left (c-\frac {c}{a^2 x^2}\right )^4} \, dx\)

\(\Big \downarrow \) 6707

\(\displaystyle \frac {a^8 \int \frac {e^{-2 \text {arctanh}(a x)} x^8}{\left (1-a^2 x^2\right )^4}dx}{c^4}\)

\(\Big \downarrow \) 6700

\(\displaystyle \frac {a^8 \int \frac {x^8}{(1-a x)^3 (a x+1)^5}dx}{c^4}\)

\(\Big \downarrow \) 99

\(\displaystyle \frac {a^8 \int \left (\frac {303}{128 a^8 (a x+1)}-\frac {99}{32 a^8 (a x+1)^2}+\frac {35}{16 a^8 (a x+1)^3}-\frac {13}{16 a^8 (a x+1)^4}+\frac {1}{8 a^8 (a x+1)^5}-\frac {1}{a^8}-\frac {47}{128 a^8 (a x-1)}-\frac {11}{64 a^8 (a x-1)^2}-\frac {1}{32 a^8 (a x-1)^3}\right )dx}{c^4}\)

\(\Big \downarrow \) 2009

\(\displaystyle \frac {a^8 \left (-\frac {11}{64 a^9 (1-a x)}+\frac {99}{32 a^9 (a x+1)}+\frac {1}{64 a^9 (1-a x)^2}-\frac {35}{32 a^9 (a x+1)^2}+\frac {13}{48 a^9 (a x+1)^3}-\frac {1}{32 a^9 (a x+1)^4}-\frac {47 \log (1-a x)}{128 a^9}+\frac {303 \log (a x+1)}{128 a^9}-\frac {x}{a^8}\right )}{c^4}\)

Input:

Int[1/(E^(2*ArcTanh[a*x])*(c - c/(a^2*x^2))^4),x]
 

Output:

(a^8*(-(x/a^8) + 1/(64*a^9*(1 - a*x)^2) - 11/(64*a^9*(1 - a*x)) - 1/(32*a^ 
9*(1 + a*x)^4) + 13/(48*a^9*(1 + a*x)^3) - 35/(32*a^9*(1 + a*x)^2) + 99/(3 
2*a^9*(1 + a*x)) - (47*Log[1 - a*x])/(128*a^9) + (303*Log[1 + a*x])/(128*a 
^9)))/c^4
 

Defintions of rubi rules used

rule 99
Int[((a_.) + (b_.)*(x_))^(m_)*((c_.) + (d_.)*(x_))^(n_)*((e_.) + (f_.)*(x_) 
)^(p_), x_] :> Int[ExpandIntegrand[(a + b*x)^m*(c + d*x)^n*(e + f*x)^p, x], 
 x] /; FreeQ[{a, b, c, d, e, f, p}, x] && IntegersQ[m, n] && (IntegerQ[p] | 
| (GtQ[m, 0] && GeQ[n, -1]))
 

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

rule 6700
Int[E^(ArcTanh[(a_.)*(x_)]*(n_.))*(x_)^(m_.)*((c_) + (d_.)*(x_)^2)^(p_.), x 
_Symbol] :> Simp[c^p   Int[x^m*(1 - a*x)^(p - n/2)*(1 + a*x)^(p + n/2), x], 
 x] /; FreeQ[{a, c, d, m, n, p}, x] && EqQ[a^2*c + d, 0] && (IntegerQ[p] || 
 GtQ[c, 0])
 

rule 6707
Int[E^(ArcTanh[(a_.)*(x_)]*(n_.))*(u_.)*((c_) + (d_.)/(x_)^2)^(p_.), x_Symb 
ol] :> Simp[d^p   Int[(u/x^(2*p))*(1 - a^2*x^2)^p*E^(n*ArcTanh[a*x]), x], x 
] /; FreeQ[{a, c, d, n}, x] && EqQ[c + a^2*d, 0] && IntegerQ[p]
 
Maple [A] (verified)

Time = 0.11 (sec) , antiderivative size = 109, normalized size of antiderivative = 0.76

method result size
default \(\frac {a^{8} \left (-\frac {x}{a^{8}}-\frac {1}{32 a^{9} \left (a x +1\right )^{4}}+\frac {13}{48 a^{9} \left (a x +1\right )^{3}}-\frac {35}{32 a^{9} \left (a x +1\right )^{2}}+\frac {99}{32 a^{9} \left (a x +1\right )}+\frac {303 \ln \left (a x +1\right )}{128 a^{9}}+\frac {1}{64 a^{9} \left (a x -1\right )^{2}}+\frac {11}{64 a^{9} \left (a x -1\right )}-\frac {47 \ln \left (a x -1\right )}{128 a^{9}}\right )}{c^{4}}\) \(109\)
risch \(-\frac {x}{c^{4}}+\frac {\frac {209 a^{4} c^{4} x^{5}}{64}+\frac {81 a^{3} c^{4} x^{4}}{32}-\frac {529 a^{2} c^{4} x^{3}}{96}-\frac {437 a \,c^{4} x^{2}}{96}+\frac {467 c^{4} x}{192}+\frac {25 c^{4}}{12 a}}{c^{8} \left (a x +1\right )^{2} \left (a^{2} x^{2}-1\right )^{2}}-\frac {47 \ln \left (a x -1\right )}{128 a \,c^{4}}+\frac {303 \ln \left (-a x -1\right )}{128 a \,c^{4}}\) \(116\)
norman \(\frac {-\frac {157 a^{3} x^{4}}{12 c}-\frac {a^{8} x^{9}}{c}+\frac {175 x}{64 c}+\frac {143 a \,x^{2}}{32 c}-\frac {1259 a^{2} x^{3}}{192 c}+\frac {715 a^{4} x^{5}}{192 c}+\frac {1219 a^{5} x^{6}}{96 c}+\frac {211 a^{6} x^{7}}{192 c}-\frac {49 a^{7} x^{8}}{12 c}}{c^{3} \left (a x +1\right )^{5} \left (a x -1\right )^{3}}-\frac {47 \ln \left (a x -1\right )}{128 a \,c^{4}}+\frac {303 \ln \left (a x +1\right )}{128 a \,c^{4}}\) \(141\)
parallelrisch \(\frac {1468 a^{3} x^{3}-1050 a x +3308 a^{4} x^{4}+38 a^{5} x^{5}-1568 x^{6} a^{6}+1818 \ln \left (a x +1\right ) x a -384 a^{7} x^{7}+909 \ln \left (a x +1\right )-1716 a^{2} x^{2}-141 \ln \left (a x -1\right )-909 \ln \left (a x +1\right ) x^{2} a^{2}+141 \ln \left (a x -1\right ) x^{4} a^{4}-909 \ln \left (a x +1\right ) x^{4} a^{4}-3636 \ln \left (a x +1\right ) x^{3} a^{3}-282 a \ln \left (a x -1\right ) x +564 a^{3} \ln \left (a x -1\right ) x^{3}-141 \ln \left (a x -1\right ) x^{6} a^{6}+909 \ln \left (a x +1\right ) x^{6} a^{6}-282 \ln \left (a x -1\right ) x^{5} a^{5}+1818 \ln \left (a x +1\right ) x^{5} a^{5}+141 a^{2} \ln \left (a x -1\right ) x^{2}}{384 c^{4} \left (a x +1\right )^{2} \left (a^{2} x^{2}-1\right )^{2} a}\) \(256\)

Input:

int(1/(a*x+1)^2*(-a^2*x^2+1)/(c-c/a^2/x^2)^4,x,method=_RETURNVERBOSE)
 

Output:

a^8/c^4*(-x/a^8-1/32/a^9/(a*x+1)^4+13/48/a^9/(a*x+1)^3-35/32/a^9/(a*x+1)^2 
+99/32/a^9/(a*x+1)+303/128/a^9*ln(a*x+1)+1/64/a^9/(a*x-1)^2+11/64/a^9/(a*x 
-1)-47/128/a^9*ln(a*x-1))
 

Fricas [A] (verification not implemented)

Time = 0.07 (sec) , antiderivative size = 233, normalized size of antiderivative = 1.62 \[ \int \frac {e^{-2 \text {arctanh}(a x)}}{\left (c-\frac {c}{a^2 x^2}\right )^4} \, dx=-\frac {384 \, a^{7} x^{7} + 768 \, a^{6} x^{6} - 1638 \, a^{5} x^{5} - 2508 \, a^{4} x^{4} + 1732 \, a^{3} x^{3} + 2516 \, a^{2} x^{2} - 550 \, a x - 909 \, {\left (a^{6} x^{6} + 2 \, a^{5} x^{5} - a^{4} x^{4} - 4 \, a^{3} x^{3} - a^{2} x^{2} + 2 \, a x + 1\right )} \log \left (a x + 1\right ) + 141 \, {\left (a^{6} x^{6} + 2 \, a^{5} x^{5} - a^{4} x^{4} - 4 \, a^{3} x^{3} - a^{2} x^{2} + 2 \, a x + 1\right )} \log \left (a x - 1\right ) - 800}{384 \, {\left (a^{7} c^{4} x^{6} + 2 \, a^{6} c^{4} x^{5} - a^{5} c^{4} x^{4} - 4 \, a^{4} c^{4} x^{3} - a^{3} c^{4} x^{2} + 2 \, a^{2} c^{4} x + a c^{4}\right )}} \] Input:

integrate(1/(a*x+1)^2*(-a^2*x^2+1)/(c-c/a^2/x^2)^4,x, algorithm="fricas")
 

Output:

-1/384*(384*a^7*x^7 + 768*a^6*x^6 - 1638*a^5*x^5 - 2508*a^4*x^4 + 1732*a^3 
*x^3 + 2516*a^2*x^2 - 550*a*x - 909*(a^6*x^6 + 2*a^5*x^5 - a^4*x^4 - 4*a^3 
*x^3 - a^2*x^2 + 2*a*x + 1)*log(a*x + 1) + 141*(a^6*x^6 + 2*a^5*x^5 - a^4* 
x^4 - 4*a^3*x^3 - a^2*x^2 + 2*a*x + 1)*log(a*x - 1) - 800)/(a^7*c^4*x^6 + 
2*a^6*c^4*x^5 - a^5*c^4*x^4 - 4*a^4*c^4*x^3 - a^3*c^4*x^2 + 2*a^2*c^4*x + 
a*c^4)
 

Sympy [A] (verification not implemented)

Time = 0.66 (sec) , antiderivative size = 158, normalized size of antiderivative = 1.10 \[ \int \frac {e^{-2 \text {arctanh}(a x)}}{\left (c-\frac {c}{a^2 x^2}\right )^4} \, dx=- a^{8} \left (\frac {- 627 a^{5} x^{5} - 486 a^{4} x^{4} + 1058 a^{3} x^{3} + 874 a^{2} x^{2} - 467 a x - 400}{192 a^{15} c^{4} x^{6} + 384 a^{14} c^{4} x^{5} - 192 a^{13} c^{4} x^{4} - 768 a^{12} c^{4} x^{3} - 192 a^{11} c^{4} x^{2} + 384 a^{10} c^{4} x + 192 a^{9} c^{4}} + \frac {x}{a^{8} c^{4}} + \frac {\frac {47 \log {\left (x - \frac {1}{a} \right )}}{128} - \frac {303 \log {\left (x + \frac {1}{a} \right )}}{128}}{a^{9} c^{4}}\right ) \] Input:

integrate(1/(a*x+1)**2*(-a**2*x**2+1)/(c-c/a**2/x**2)**4,x)
 

Output:

-a**8*((-627*a**5*x**5 - 486*a**4*x**4 + 1058*a**3*x**3 + 874*a**2*x**2 - 
467*a*x - 400)/(192*a**15*c**4*x**6 + 384*a**14*c**4*x**5 - 192*a**13*c**4 
*x**4 - 768*a**12*c**4*x**3 - 192*a**11*c**4*x**2 + 384*a**10*c**4*x + 192 
*a**9*c**4) + x/(a**8*c**4) + (47*log(x - 1/a)/128 - 303*log(x + 1/a)/128) 
/(a**9*c**4))
 

Maxima [A] (verification not implemented)

Time = 0.04 (sec) , antiderivative size = 146, normalized size of antiderivative = 1.01 \[ \int \frac {e^{-2 \text {arctanh}(a x)}}{\left (c-\frac {c}{a^2 x^2}\right )^4} \, dx=\frac {627 \, a^{5} x^{5} + 486 \, a^{4} x^{4} - 1058 \, a^{3} x^{3} - 874 \, a^{2} x^{2} + 467 \, a x + 400}{192 \, {\left (a^{7} c^{4} x^{6} + 2 \, a^{6} c^{4} x^{5} - a^{5} c^{4} x^{4} - 4 \, a^{4} c^{4} x^{3} - a^{3} c^{4} x^{2} + 2 \, a^{2} c^{4} x + a c^{4}\right )}} - \frac {x}{c^{4}} + \frac {303 \, \log \left (a x + 1\right )}{128 \, a c^{4}} - \frac {47 \, \log \left (a x - 1\right )}{128 \, a c^{4}} \] Input:

integrate(1/(a*x+1)^2*(-a^2*x^2+1)/(c-c/a^2/x^2)^4,x, algorithm="maxima")
 

Output:

1/192*(627*a^5*x^5 + 486*a^4*x^4 - 1058*a^3*x^3 - 874*a^2*x^2 + 467*a*x + 
400)/(a^7*c^4*x^6 + 2*a^6*c^4*x^5 - a^5*c^4*x^4 - 4*a^4*c^4*x^3 - a^3*c^4* 
x^2 + 2*a^2*c^4*x + a*c^4) - x/c^4 + 303/128*log(a*x + 1)/(a*c^4) - 47/128 
*log(a*x - 1)/(a*c^4)
 

Giac [A] (verification not implemented)

Time = 0.13 (sec) , antiderivative size = 164, normalized size of antiderivative = 1.14 \[ \int \frac {e^{-2 \text {arctanh}(a x)}}{\left (c-\frac {c}{a^2 x^2}\right )^4} \, dx=-\frac {2 \, \log \left (\frac {{\left | a x + 1 \right |}}{{\left (a x + 1\right )}^{2} {\left | a \right |}}\right )}{a c^{4}} - \frac {47 \, \log \left ({\left | -\frac {2}{a x + 1} + 1 \right |}\right )}{128 \, a c^{4}} + \frac {{\left (a x + 1\right )} {\left (\frac {1045}{a x + 1} - \frac {1064}{{\left (a x + 1\right )}^{2}} - 256\right )}}{256 \, a c^{4} {\left (\frac {2}{a x + 1} - 1\right )}^{2}} + \frac {\frac {297 \, a^{19} c^{12}}{a x + 1} - \frac {105 \, a^{19} c^{12}}{{\left (a x + 1\right )}^{2}} + \frac {26 \, a^{19} c^{12}}{{\left (a x + 1\right )}^{3}} - \frac {3 \, a^{19} c^{12}}{{\left (a x + 1\right )}^{4}}}{96 \, a^{20} c^{16}} \] Input:

integrate(1/(a*x+1)^2*(-a^2*x^2+1)/(c-c/a^2/x^2)^4,x, algorithm="giac")
 

Output:

-2*log(abs(a*x + 1)/((a*x + 1)^2*abs(a)))/(a*c^4) - 47/128*log(abs(-2/(a*x 
 + 1) + 1))/(a*c^4) + 1/256*(a*x + 1)*(1045/(a*x + 1) - 1064/(a*x + 1)^2 - 
 256)/(a*c^4*(2/(a*x + 1) - 1)^2) + 1/96*(297*a^19*c^12/(a*x + 1) - 105*a^ 
19*c^12/(a*x + 1)^2 + 26*a^19*c^12/(a*x + 1)^3 - 3*a^19*c^12/(a*x + 1)^4)/ 
(a^20*c^16)
 

Mupad [B] (verification not implemented)

Time = 0.14 (sec) , antiderivative size = 142, normalized size of antiderivative = 0.99 \[ \int \frac {e^{-2 \text {arctanh}(a x)}}{\left (c-\frac {c}{a^2 x^2}\right )^4} \, dx=\frac {\frac {467\,x}{192}-\frac {437\,a\,x^2}{96}+\frac {25}{12\,a}-\frac {529\,a^2\,x^3}{96}+\frac {81\,a^3\,x^4}{32}+\frac {209\,a^4\,x^5}{64}}{a^6\,c^4\,x^6+2\,a^5\,c^4\,x^5-a^4\,c^4\,x^4-4\,a^3\,c^4\,x^3-a^2\,c^4\,x^2+2\,a\,c^4\,x+c^4}-\frac {x}{c^4}-\frac {47\,\ln \left (a\,x-1\right )}{128\,a\,c^4}+\frac {303\,\ln \left (a\,x+1\right )}{128\,a\,c^4} \] Input:

int(-(a^2*x^2 - 1)/((c - c/(a^2*x^2))^4*(a*x + 1)^2),x)
 

Output:

((467*x)/192 - (437*a*x^2)/96 + 25/(12*a) - (529*a^2*x^3)/96 + (81*a^3*x^4 
)/32 + (209*a^4*x^5)/64)/(c^4 - a^2*c^4*x^2 - 4*a^3*c^4*x^3 - a^4*c^4*x^4 
+ 2*a^5*c^4*x^5 + a^6*c^4*x^6 + 2*a*c^4*x) - x/c^4 - (47*log(a*x - 1))/(12 
8*a*c^4) + (303*log(a*x + 1))/(128*a*c^4)
 

Reduce [B] (verification not implemented)

Time = 0.14 (sec) , antiderivative size = 277, normalized size of antiderivative = 1.92 \[ \int \frac {e^{-2 \text {arctanh}(a x)}}{\left (c-\frac {c}{a^2 x^2}\right )^4} \, dx=\frac {-141 \,\mathrm {log}\left (a x -1\right ) a^{6} x^{6}-282 \,\mathrm {log}\left (a x -1\right ) a^{5} x^{5}+141 \,\mathrm {log}\left (a x -1\right ) a^{4} x^{4}+564 \,\mathrm {log}\left (a x -1\right ) a^{3} x^{3}+141 \,\mathrm {log}\left (a x -1\right ) a^{2} x^{2}-282 \,\mathrm {log}\left (a x -1\right ) a x -141 \,\mathrm {log}\left (a x -1\right )+909 \,\mathrm {log}\left (a x +1\right ) a^{6} x^{6}+1818 \,\mathrm {log}\left (a x +1\right ) a^{5} x^{5}-909 \,\mathrm {log}\left (a x +1\right ) a^{4} x^{4}-3636 \,\mathrm {log}\left (a x +1\right ) a^{3} x^{3}-909 \,\mathrm {log}\left (a x +1\right ) a^{2} x^{2}+1818 \,\mathrm {log}\left (a x +1\right ) a x +909 \,\mathrm {log}\left (a x +1\right )-384 a^{7} x^{7}-1587 a^{6} x^{6}+3327 a^{4} x^{4}+1544 a^{3} x^{3}-1697 a^{2} x^{2}-1088 a x -19}{384 a \,c^{4} \left (a^{6} x^{6}+2 a^{5} x^{5}-a^{4} x^{4}-4 a^{3} x^{3}-a^{2} x^{2}+2 a x +1\right )} \] Input:

int(1/(a*x+1)^2*(-a^2*x^2+1)/(c-c/a^2/x^2)^4,x)
 

Output:

( - 141*log(a*x - 1)*a**6*x**6 - 282*log(a*x - 1)*a**5*x**5 + 141*log(a*x 
- 1)*a**4*x**4 + 564*log(a*x - 1)*a**3*x**3 + 141*log(a*x - 1)*a**2*x**2 - 
 282*log(a*x - 1)*a*x - 141*log(a*x - 1) + 909*log(a*x + 1)*a**6*x**6 + 18 
18*log(a*x + 1)*a**5*x**5 - 909*log(a*x + 1)*a**4*x**4 - 3636*log(a*x + 1) 
*a**3*x**3 - 909*log(a*x + 1)*a**2*x**2 + 1818*log(a*x + 1)*a*x + 909*log( 
a*x + 1) - 384*a**7*x**7 - 1587*a**6*x**6 + 3327*a**4*x**4 + 1544*a**3*x** 
3 - 1697*a**2*x**2 - 1088*a*x - 19)/(384*a*c**4*(a**6*x**6 + 2*a**5*x**5 - 
 a**4*x**4 - 4*a**3*x**3 - a**2*x**2 + 2*a*x + 1))