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

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 = 146 \[ \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}{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 {1}{64 a c^4 (1+a x)^2}+\frac {11}{64 a c^4 (1+a x)}-\frac {303 \log (1-a x)}{128 a c^4}+\frac {47 \log (1+a x)}{128 a c^4} \] Output:

-x/c^4+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)-1/64/a/c^4/(a*x+1)^2+11/64/a/c^4/(a*x+1)-303/128*ln( 
-a*x+1)/a/c^4+47/128*ln(a*x+1)/a/c^4
                                                                                    
                                                                                    
 

Mathematica [A] (verified)

Time = 0.12 (sec) , antiderivative size = 123, 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-909 (-1+a x)^4 (1+a x)^2 \log (1-a x)+141 (-1+a x)^4 (1+a x)^2 \log (1+a x)}{384 a c^4 (-1+a x)^4 (1+a x)^2} \] Input:

Integrate[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 - 909*(-1 + a*x)^4*(1 + a*x)^2*Log[1 - a*x] 
+ 141*(-1 + a*x)^4*(1 + a*x)^2*Log[1 + a*x])/(384*a*c^4*(-1 + a*x)^4*(1 + 
a*x)^2)
 

Rubi [A] (verified)

Time = 0.79 (sec) , antiderivative size = 129, 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)^5 (a x+1)^3}dx}{c^4}\)

\(\Big \downarrow \) 99

\(\displaystyle \frac {a^8 \int \left (\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}-\frac {1}{a^8}-\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}\right )dx}{c^4}\)

\(\Big \downarrow \) 2009

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

Input:

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

Output:

(a^8*(-(x/a^8) + 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/(32*a^9*(1 - a*x)) - 1/(64*a^9*(1 + a*x)^2) + 11/(6 
4*a^9*(1 + a*x)) - (303*Log[1 - a*x])/(128*a^9) + (47*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.75

method result size
default \(\frac {a^{8} \left (-\frac {x}{a^{8}}-\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}}+\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}}\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 {175 x}{64 c}-\frac {2555 a^{4} x^{5}}{192 c}+\frac {465 a^{6} x^{7}}{64 c}-\frac {a^{8} x^{9}}{c}+\frac {1925 a^{2} x^{3}}{192 c}-\frac {a \,x^{2}}{c}+\frac {7 a^{3} x^{4}}{2 c}-\frac {13 a^{5} x^{6}}{3 c}+\frac {25 a^{7} x^{8}}{12 c}}{\left (a^{2} x^{2}-1\right )^{4} c^{3}}-\frac {303 \ln \left (a x -1\right )}{128 a \,c^{4}}+\frac {47 \ln \left (a x +1\right )}{128 a \,c^{4}}\) \(138\)
parallelrisch \(\frac {1468 a^{3} x^{3}-1050 a x -3308 a^{4} x^{4}+38 a^{5} x^{5}-909 \ln \left (a x -1\right ) x^{6} a^{6}+141 \ln \left (a x +1\right ) x^{6} a^{6}+1818 \ln \left (a x -1\right ) x^{5} a^{5}-282 \ln \left (a x +1\right ) x^{5} a^{5}+1568 x^{6} a^{6}-282 \ln \left (a x +1\right ) x a -384 a^{7} x^{7}+141 \ln \left (a x +1\right )+1716 a^{2} x^{2}-909 \ln \left (a x -1\right )-141 \ln \left (a x +1\right ) x^{2} a^{2}+909 \ln \left (a x -1\right ) x^{4} a^{4}-141 \ln \left (a x +1\right ) x^{4} a^{4}+564 \ln \left (a x +1\right ) x^{3} a^{3}+1818 a \ln \left (a x -1\right ) x -3636 a^{3} \ln \left (a x -1\right ) x^{3}+909 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((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/64/a^9/(a*x+1)^2+11/64/a^9/(a*x+1)+47/128/a^9*ln(a*x+1)+ 
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))
 

Fricas [A] (verification not implemented)

Time = 0.07 (sec) , antiderivative size = 233, normalized size of antiderivative = 1.60 \[ \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 - 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 ) + 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 ) + 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((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 - 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) + 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) + 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.83 (sec) , antiderivative size = 158, normalized size of antiderivative = 1.08 \[ \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 {303 \log {\left (x - \frac {1}{a} \right )}}{128} - \frac {47 \log {\left (x + \frac {1}{a} \right )}}{128}}{a^{9} c^{4}}\right ) \] Input:

integrate((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) + (303*log(x - 1/a)/128 - 47*log(x + 1/a)/128) 
/(a**9*c**4))
 

Maxima [A] (verification not implemented)

Time = 0.06 (sec) , antiderivative size = 146, normalized size of antiderivative = 1.00 \[ \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 {47 \, \log \left (a x + 1\right )}{128 \, a c^{4}} - \frac {303 \, \log \left (a x - 1\right )}{128 \, a c^{4}} \] Input:

integrate((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 + 47/128*log(a*x + 1)/(a*c^4) - 303/128 
*log(a*x - 1)/(a*c^4)
 

Giac [A] (verification not implemented)

Time = 0.13 (sec) , antiderivative size = 97, normalized size of antiderivative = 0.66 \[ \int \frac {e^{2 \text {arctanh}(a x)}}{\left (c-\frac {c}{a^2 x^2}\right )^4} \, dx=-\frac {x}{c^{4}} + \frac {47 \, \log \left ({\left | a x + 1 \right |}\right )}{128 \, a c^{4}} - \frac {303 \, \log \left ({\left | a x - 1 \right |}\right )}{128 \, a c^{4}} + \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 x + 1\right )}^{2} {\left (a x - 1\right )}^{4} a c^{4}} \] Input:

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

Output:

-x/c^4 + 47/128*log(abs(a*x + 1))/(a*c^4) - 303/128*log(abs(a*x - 1))/(a*c 
^4) + 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*x + 1)^2*(a*x - 1)^4*a*c^4)
 

Mupad [B] (verification not implemented)

Time = 22.69 (sec) , antiderivative size = 144, 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 {47\,\ln \left (a\,x+1\right )}{128\,a\,c^4}-\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 {303\,\ln \left (a\,x-1\right )}{128\,a\,c^4}-\frac {x}{c^4} \] Input:

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

Output:

(47*log(a*x + 1))/(128*a*c^4) - ((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)/(a^2*c^4*x^2 - c^ 
4 - 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) 
 - (303*log(a*x - 1))/(128*a*c^4) - x/c^4
 

Reduce [B] (verification not implemented)

Time = 0.14 (sec) , antiderivative size = 277, normalized size of antiderivative = 1.90 \[ \int \frac {e^{2 \text {arctanh}(a x)}}{\left (c-\frac {c}{a^2 x^2}\right )^4} \, dx=\frac {-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 )+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 )-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((a*x+1)^2/(-a^2*x^2+1)/(c-c/a^2/x^2)^4,x)
 

Output:

( - 909*log(a*x - 1)*a**6*x**6 + 1818*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) + 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) - 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))