3.9.24 \(\int e^{-3 \coth ^{-1}(a x)} (c-\frac {c}{a^2 x^2})^2 \, dx\) [824]

3.9.24.1 Optimal result
3.9.24.2 Mathematica [A] (verified)
3.9.24.3 Rubi [A] (verified)
3.9.24.4 Maple [A] (verified)
3.9.24.5 Fricas [A] (verification not implemented)
3.9.24.6 Sympy [F]
3.9.24.7 Maxima [A] (verification not implemented)
3.9.24.8 Giac [A] (verification not implemented)
3.9.24.9 Mupad [B] (verification not implemented)

3.9.24.1 Optimal result

Integrand size = 22, antiderivative size = 195 \[ \int e^{-3 \coth ^{-1}(a x)} \left (c-\frac {c}{a^2 x^2}\right )^2 \, dx=\frac {5 c^2 \sqrt {1-\frac {1}{a x}} \sqrt {1+\frac {1}{a x}}}{2 a}+\frac {11 c^2 \left (1-\frac {1}{a x}\right )^{3/2} \sqrt {1+\frac {1}{a x}}}{6 a}+\frac {4 c^2 \left (1-\frac {1}{a x}\right )^{5/2} \sqrt {1+\frac {1}{a x}}}{3 a}+c^2 \left (1-\frac {1}{a x}\right )^{7/2} \sqrt {1+\frac {1}{a x}} x-\frac {c^2 \csc ^{-1}(a x)}{2 a}-\frac {3 c^2 \text {arctanh}\left (\sqrt {1-\frac {1}{a x}} \sqrt {1+\frac {1}{a x}}\right )}{a} \]

output
-1/2*c^2*arccsc(a*x)/a-3*c^2*arctanh((1-1/a/x)^(1/2)*(1+1/a/x)^(1/2))/a+11 
/6*c^2*(1-1/a/x)^(3/2)*(1+1/a/x)^(1/2)/a+4/3*c^2*(1-1/a/x)^(5/2)*(1+1/a/x) 
^(1/2)/a+c^2*(1-1/a/x)^(7/2)*x*(1+1/a/x)^(1/2)+5/2*c^2*(1-1/a/x)^(1/2)*(1+ 
1/a/x)^(1/2)/a
 
3.9.24.2 Mathematica [A] (verified)

Time = 0.20 (sec) , antiderivative size = 94, normalized size of antiderivative = 0.48 \[ \int e^{-3 \coth ^{-1}(a x)} \left (c-\frac {c}{a^2 x^2}\right )^2 \, dx=\frac {c^2 \left (\sqrt {1-\frac {1}{a^2 x^2}} \left (2-9 a x+16 a^2 x^2+6 a^3 x^3\right )-3 a^2 x^2 \arcsin \left (\frac {1}{a x}\right )-18 a^2 x^2 \log \left (\left (1+\sqrt {1-\frac {1}{a^2 x^2}}\right ) x\right )\right )}{6 a^3 x^2} \]

input
Integrate[(c - c/(a^2*x^2))^2/E^(3*ArcCoth[a*x]),x]
 
output
(c^2*(Sqrt[1 - 1/(a^2*x^2)]*(2 - 9*a*x + 16*a^2*x^2 + 6*a^3*x^3) - 3*a^2*x 
^2*ArcSin[1/(a*x)] - 18*a^2*x^2*Log[(1 + Sqrt[1 - 1/(a^2*x^2)])*x]))/(6*a^ 
3*x^2)
 
3.9.24.3 Rubi [A] (verified)

Time = 0.39 (sec) , antiderivative size = 189, normalized size of antiderivative = 0.97, number of steps used = 16, number of rules used = 15, \(\frac {\text {number of rules}}{\text {integrand size}}\) = 0.682, Rules used = {6748, 108, 25, 27, 171, 27, 171, 27, 171, 27, 175, 39, 103, 221, 223}

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 (c-\frac {c}{a^2 x^2}\right )^2 e^{-3 \coth ^{-1}(a x)} \, dx\)

\(\Big \downarrow \) 6748

\(\displaystyle -c^2 \int \left (1-\frac {1}{a x}\right )^{7/2} \sqrt {1+\frac {1}{a x}} x^2d\frac {1}{x}\)

\(\Big \downarrow \) 108

\(\displaystyle -c^2 \left (\int -\frac {\left (3 a+\frac {4}{x}\right ) \left (1-\frac {1}{a x}\right )^{5/2} x}{a^2 \sqrt {1+\frac {1}{a x}}}d\frac {1}{x}-x \left (1-\frac {1}{a x}\right )^{7/2} \sqrt {\frac {1}{a x}+1}\right )\)

\(\Big \downarrow \) 25

\(\displaystyle -c^2 \left (x \left (-\sqrt {\frac {1}{a x}+1}\right ) \left (1-\frac {1}{a x}\right )^{7/2}-\int \frac {\left (3 a+\frac {4}{x}\right ) \left (1-\frac {1}{a x}\right )^{5/2} x}{a^2 \sqrt {1+\frac {1}{a x}}}d\frac {1}{x}\right )\)

\(\Big \downarrow \) 27

\(\displaystyle -c^2 \left (x \left (-\sqrt {\frac {1}{a x}+1}\right ) \left (1-\frac {1}{a x}\right )^{7/2}-\frac {\int \frac {\left (3 a+\frac {4}{x}\right ) \left (1-\frac {1}{a x}\right )^{5/2} x}{\sqrt {1+\frac {1}{a x}}}d\frac {1}{x}}{a^2}\right )\)

\(\Big \downarrow \) 171

\(\displaystyle -c^2 \left (x \left (-\sqrt {\frac {1}{a x}+1}\right ) \left (1-\frac {1}{a x}\right )^{7/2}-\frac {\frac {1}{3} a \int \frac {\left (9 a+\frac {11}{x}\right ) \left (1-\frac {1}{a x}\right )^{3/2} x}{a \sqrt {1+\frac {1}{a x}}}d\frac {1}{x}+\frac {4}{3} a \sqrt {\frac {1}{a x}+1} \left (1-\frac {1}{a x}\right )^{5/2}}{a^2}\right )\)

\(\Big \downarrow \) 27

\(\displaystyle -c^2 \left (x \left (-\sqrt {\frac {1}{a x}+1}\right ) \left (1-\frac {1}{a x}\right )^{7/2}-\frac {\frac {1}{3} \int \frac {\left (9 a+\frac {11}{x}\right ) \left (1-\frac {1}{a x}\right )^{3/2} x}{\sqrt {1+\frac {1}{a x}}}d\frac {1}{x}+\frac {4}{3} a \sqrt {\frac {1}{a x}+1} \left (1-\frac {1}{a x}\right )^{5/2}}{a^2}\right )\)

\(\Big \downarrow \) 171

\(\displaystyle -c^2 \left (x \left (-\sqrt {\frac {1}{a x}+1}\right ) \left (1-\frac {1}{a x}\right )^{7/2}-\frac {\frac {1}{3} \left (\frac {1}{2} a \int \frac {3 \left (6 a+\frac {5}{x}\right ) \sqrt {1-\frac {1}{a x}} x}{a \sqrt {1+\frac {1}{a x}}}d\frac {1}{x}+\frac {11}{2} a \sqrt {\frac {1}{a x}+1} \left (1-\frac {1}{a x}\right )^{3/2}\right )+\frac {4}{3} a \sqrt {\frac {1}{a x}+1} \left (1-\frac {1}{a x}\right )^{5/2}}{a^2}\right )\)

\(\Big \downarrow \) 27

\(\displaystyle -c^2 \left (x \left (-\sqrt {\frac {1}{a x}+1}\right ) \left (1-\frac {1}{a x}\right )^{7/2}-\frac {\frac {1}{3} \left (\frac {3}{2} \int \frac {\left (6 a+\frac {5}{x}\right ) \sqrt {1-\frac {1}{a x}} x}{\sqrt {1+\frac {1}{a x}}}d\frac {1}{x}+\frac {11}{2} a \sqrt {\frac {1}{a x}+1} \left (1-\frac {1}{a x}\right )^{3/2}\right )+\frac {4}{3} a \sqrt {\frac {1}{a x}+1} \left (1-\frac {1}{a x}\right )^{5/2}}{a^2}\right )\)

\(\Big \downarrow \) 171

\(\displaystyle -c^2 \left (x \left (-\sqrt {\frac {1}{a x}+1}\right ) \left (1-\frac {1}{a x}\right )^{7/2}-\frac {\frac {1}{3} \left (\frac {3}{2} \left (a \int \frac {\left (6 a-\frac {1}{x}\right ) x}{a \sqrt {1-\frac {1}{a x}} \sqrt {1+\frac {1}{a x}}}d\frac {1}{x}+5 a \sqrt {1-\frac {1}{a x}} \sqrt {\frac {1}{a x}+1}\right )+\frac {11}{2} a \sqrt {\frac {1}{a x}+1} \left (1-\frac {1}{a x}\right )^{3/2}\right )+\frac {4}{3} a \sqrt {\frac {1}{a x}+1} \left (1-\frac {1}{a x}\right )^{5/2}}{a^2}\right )\)

\(\Big \downarrow \) 27

\(\displaystyle -c^2 \left (x \left (-\sqrt {\frac {1}{a x}+1}\right ) \left (1-\frac {1}{a x}\right )^{7/2}-\frac {\frac {1}{3} \left (\frac {3}{2} \left (\int \frac {\left (6 a-\frac {1}{x}\right ) x}{\sqrt {1-\frac {1}{a x}} \sqrt {1+\frac {1}{a x}}}d\frac {1}{x}+5 a \sqrt {1-\frac {1}{a x}} \sqrt {\frac {1}{a x}+1}\right )+\frac {11}{2} a \sqrt {\frac {1}{a x}+1} \left (1-\frac {1}{a x}\right )^{3/2}\right )+\frac {4}{3} a \sqrt {\frac {1}{a x}+1} \left (1-\frac {1}{a x}\right )^{5/2}}{a^2}\right )\)

\(\Big \downarrow \) 175

\(\displaystyle -c^2 \left (x \left (-\sqrt {\frac {1}{a x}+1}\right ) \left (1-\frac {1}{a x}\right )^{7/2}-\frac {\frac {1}{3} \left (\frac {3}{2} \left (6 a \int \frac {x}{\sqrt {1-\frac {1}{a x}} \sqrt {1+\frac {1}{a x}}}d\frac {1}{x}-\int \frac {1}{\sqrt {1-\frac {1}{a x}} \sqrt {1+\frac {1}{a x}}}d\frac {1}{x}+5 a \sqrt {1-\frac {1}{a x}} \sqrt {\frac {1}{a x}+1}\right )+\frac {11}{2} a \sqrt {\frac {1}{a x}+1} \left (1-\frac {1}{a x}\right )^{3/2}\right )+\frac {4}{3} a \sqrt {\frac {1}{a x}+1} \left (1-\frac {1}{a x}\right )^{5/2}}{a^2}\right )\)

\(\Big \downarrow \) 39

\(\displaystyle -c^2 \left (x \left (-\sqrt {\frac {1}{a x}+1}\right ) \left (1-\frac {1}{a x}\right )^{7/2}-\frac {\frac {1}{3} \left (\frac {3}{2} \left (-\int \frac {1}{\sqrt {1-\frac {1}{a^2 x^2}}}d\frac {1}{x}+6 a \int \frac {x}{\sqrt {1-\frac {1}{a x}} \sqrt {1+\frac {1}{a x}}}d\frac {1}{x}+5 a \sqrt {1-\frac {1}{a x}} \sqrt {\frac {1}{a x}+1}\right )+\frac {11}{2} a \sqrt {\frac {1}{a x}+1} \left (1-\frac {1}{a x}\right )^{3/2}\right )+\frac {4}{3} a \sqrt {\frac {1}{a x}+1} \left (1-\frac {1}{a x}\right )^{5/2}}{a^2}\right )\)

\(\Big \downarrow \) 103

\(\displaystyle -c^2 \left (x \left (-\sqrt {\frac {1}{a x}+1}\right ) \left (1-\frac {1}{a x}\right )^{7/2}-\frac {\frac {1}{3} \left (\frac {3}{2} \left (-\int \frac {1}{\sqrt {1-\frac {1}{a^2 x^2}}}d\frac {1}{x}-6 \int \frac {1}{\frac {1}{a}-\frac {1}{a x^2}}d\left (\sqrt {1-\frac {1}{a x}} \sqrt {1+\frac {1}{a x}}\right )+5 a \sqrt {1-\frac {1}{a x}} \sqrt {\frac {1}{a x}+1}\right )+\frac {11}{2} a \sqrt {\frac {1}{a x}+1} \left (1-\frac {1}{a x}\right )^{3/2}\right )+\frac {4}{3} a \sqrt {\frac {1}{a x}+1} \left (1-\frac {1}{a x}\right )^{5/2}}{a^2}\right )\)

\(\Big \downarrow \) 221

\(\displaystyle -c^2 \left (x \left (-\sqrt {\frac {1}{a x}+1}\right ) \left (1-\frac {1}{a x}\right )^{7/2}-\frac {\frac {1}{3} \left (\frac {3}{2} \left (-\int \frac {1}{\sqrt {1-\frac {1}{a^2 x^2}}}d\frac {1}{x}-6 a \text {arctanh}\left (\sqrt {1-\frac {1}{a x}} \sqrt {\frac {1}{a x}+1}\right )+5 a \sqrt {1-\frac {1}{a x}} \sqrt {\frac {1}{a x}+1}\right )+\frac {11}{2} a \sqrt {\frac {1}{a x}+1} \left (1-\frac {1}{a x}\right )^{3/2}\right )+\frac {4}{3} a \sqrt {\frac {1}{a x}+1} \left (1-\frac {1}{a x}\right )^{5/2}}{a^2}\right )\)

\(\Big \downarrow \) 223

\(\displaystyle -c^2 \left (x \left (-\sqrt {\frac {1}{a x}+1}\right ) \left (1-\frac {1}{a x}\right )^{7/2}-\frac {\frac {1}{3} \left (\frac {3}{2} \left (a \left (-\arcsin \left (\frac {1}{a x}\right )\right )-6 a \text {arctanh}\left (\sqrt {1-\frac {1}{a x}} \sqrt {\frac {1}{a x}+1}\right )+5 a \sqrt {1-\frac {1}{a x}} \sqrt {\frac {1}{a x}+1}\right )+\frac {11}{2} a \sqrt {\frac {1}{a x}+1} \left (1-\frac {1}{a x}\right )^{3/2}\right )+\frac {4}{3} a \sqrt {\frac {1}{a x}+1} \left (1-\frac {1}{a x}\right )^{5/2}}{a^2}\right )\)

input
Int[(c - c/(a^2*x^2))^2/E^(3*ArcCoth[a*x]),x]
 
output
-(c^2*(-((1 - 1/(a*x))^(7/2)*Sqrt[1 + 1/(a*x)]*x) - ((4*a*(1 - 1/(a*x))^(5 
/2)*Sqrt[1 + 1/(a*x)])/3 + ((11*a*(1 - 1/(a*x))^(3/2)*Sqrt[1 + 1/(a*x)])/2 
 + (3*(5*a*Sqrt[1 - 1/(a*x)]*Sqrt[1 + 1/(a*x)] - a*ArcSin[1/(a*x)] - 6*a*A 
rcTanh[Sqrt[1 - 1/(a*x)]*Sqrt[1 + 1/(a*x)]]))/2)/3)/a^2))
 

3.9.24.3.1 Defintions of rubi rules used

rule 25
Int[-(Fx_), x_Symbol] :> Simp[Identity[-1]   Int[Fx, x], x]
 

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 39
Int[((a_) + (b_.)*(x_))^(m_.)*((c_) + (d_.)*(x_))^(m_.), x_Symbol] :> Int[( 
a*c + b*d*x^2)^m, x] /; FreeQ[{a, b, c, d, m}, x] && EqQ[b*c + a*d, 0] && ( 
IntegerQ[m] || (GtQ[a, 0] && GtQ[c, 0]))
 

rule 103
Int[1/(Sqrt[(a_.) + (b_.)*(x_)]*Sqrt[(c_.) + (d_.)*(x_)]*((e_.) + (f_.)*(x_ 
))), x_] :> Simp[b*f   Subst[Int[1/(d*(b*e - a*f)^2 + b*f^2*x^2), x], x, Sq 
rt[a + b*x]*Sqrt[c + d*x]], x] /; FreeQ[{a, b, c, d, e, f}, x] && EqQ[2*b*d 
*e - f*(b*c + a*d), 0]
 

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

rule 171
Int[((a_.) + (b_.)*(x_))^(m_)*((c_.) + (d_.)*(x_))^(n_)*((e_.) + (f_.)*(x_) 
)^(p_)*((g_.) + (h_.)*(x_)), x_] :> Simp[h*(a + b*x)^m*(c + d*x)^(n + 1)*(( 
e + f*x)^(p + 1)/(d*f*(m + n + p + 2))), x] + Simp[1/(d*f*(m + n + p + 2)) 
  Int[(a + b*x)^(m - 1)*(c + d*x)^n*(e + f*x)^p*Simp[a*d*f*g*(m + n + p + 2 
) - h*(b*c*e*m + a*(d*e*(n + 1) + c*f*(p + 1))) + (b*d*f*g*(m + n + p + 2) 
+ h*(a*d*f*m - b*(d*e*(m + n + 1) + c*f*(m + p + 1))))*x, x], x], x] /; Fre 
eQ[{a, b, c, d, e, f, g, h, n, p}, x] && GtQ[m, 0] && NeQ[m + n + p + 2, 0] 
 && IntegersQ[2*m, 2*n, 2*p]
 

rule 175
Int[(((c_.) + (d_.)*(x_))^(n_)*((e_.) + (f_.)*(x_))^(p_)*((g_.) + (h_.)*(x_ 
)))/((a_.) + (b_.)*(x_)), x_] :> Simp[h/b   Int[(c + d*x)^n*(e + f*x)^p, x] 
, x] + Simp[(b*g - a*h)/b   Int[(c + d*x)^n*((e + f*x)^p/(a + b*x)), x], x] 
 /; FreeQ[{a, b, c, d, e, f, g, h, n, p}, x]
 

rule 221
Int[((a_) + (b_.)*(x_)^2)^(-1), x_Symbol] :> Simp[(Rt[-a/b, 2]/a)*ArcTanh[x 
/Rt[-a/b, 2]], x] /; FreeQ[{a, b}, x] && NegQ[a/b]
 

rule 223
Int[1/Sqrt[(a_) + (b_.)*(x_)^2], x_Symbol] :> Simp[ArcSin[Rt[-b, 2]*(x/Sqrt 
[a])]/Rt[-b, 2], x] /; FreeQ[{a, b}, x] && GtQ[a, 0] && NegQ[b]
 

rule 6748
Int[E^(ArcCoth[(a_.)*(x_)]*(n_.))*((c_) + (d_.)/(x_)^2)^(p_.), x_Symbol] :> 
 Simp[-c^p   Subst[Int[(1 - x/a)^(p - n/2)*((1 + x/a)^(p + n/2)/x^2), x], x 
, 1/x], x] /; FreeQ[{a, c, d, n, p}, x] && EqQ[c + a^2*d, 0] &&  !IntegerQ[ 
n/2] && (IntegerQ[p] || GtQ[c, 0]) &&  !IntegersQ[2*p, p + n/2]
 
3.9.24.4 Maple [A] (verified)

Time = 0.09 (sec) , antiderivative size = 157, normalized size of antiderivative = 0.81

method result size
risch \(\frac {\left (a x +1\right ) \left (16 a^{2} x^{2}-9 a x +2\right ) c^{2} \sqrt {\frac {a x -1}{a x +1}}}{6 x^{3} a^{4}}+\frac {\left (-\frac {3 a^{4} \ln \left (\frac {a^{2} x}{\sqrt {a^{2}}}+\sqrt {a^{2} x^{2}-1}\right )}{\sqrt {a^{2}}}-\frac {a^{3} \arctan \left (\frac {1}{\sqrt {a^{2} x^{2}-1}}\right )}{2}+a^{3} \sqrt {\left (a x -1\right ) \left (a x +1\right )}\right ) c^{2} \sqrt {\frac {a x -1}{a x +1}}\, \sqrt {\left (a x -1\right ) \left (a x +1\right )}}{a^{4} \left (a x -1\right )}\) \(157\)
default \(-\frac {\left (\frac {a x -1}{a x +1}\right )^{\frac {3}{2}} \left (a x +1\right )^{2} c^{2} \left (-18 \sqrt {a^{2} x^{2}-1}\, \sqrt {a^{2}}\, a^{4} x^{4}+18 \left (a^{2} x^{2}-1\right )^{\frac {3}{2}} \sqrt {a^{2}}\, a^{2} x^{2}+3 \sqrt {a^{2} x^{2}-1}\, \sqrt {a^{2}}\, a^{3} x^{3}+3 a^{3} x^{3} \sqrt {a^{2}}\, \arctan \left (\frac {1}{\sqrt {a^{2} x^{2}-1}}\right )+18 \ln \left (\frac {a^{2} x +\sqrt {a^{2} x^{2}-1}\, \sqrt {a^{2}}}{\sqrt {a^{2}}}\right ) a^{4} x^{3}-9 \left (a^{2} x^{2}-1\right )^{\frac {3}{2}} \sqrt {a^{2}}\, a x +2 \left (a^{2} x^{2}-1\right )^{\frac {3}{2}} \sqrt {a^{2}}\right )}{6 \left (a x -1\right ) \sqrt {\left (a x -1\right ) \left (a x +1\right )}\, a^{4} x^{3} \sqrt {a^{2}}}\) \(233\)

input
int((c-c/a^2/x^2)^2*((a*x-1)/(a*x+1))^(3/2),x,method=_RETURNVERBOSE)
 
output
1/6*(a*x+1)*(16*a^2*x^2-9*a*x+2)/x^3*c^2/a^4*((a*x-1)/(a*x+1))^(1/2)+(-3*a 
^4*ln(a^2*x/(a^2)^(1/2)+(a^2*x^2-1)^(1/2))/(a^2)^(1/2)-1/2*a^3*arctan(1/(a 
^2*x^2-1)^(1/2))+a^3*((a*x-1)*(a*x+1))^(1/2))*c^2/a^4*((a*x-1)/(a*x+1))^(1 
/2)/(a*x-1)*((a*x-1)*(a*x+1))^(1/2)
 
3.9.24.5 Fricas [A] (verification not implemented)

Time = 0.27 (sec) , antiderivative size = 156, normalized size of antiderivative = 0.80 \[ \int e^{-3 \coth ^{-1}(a x)} \left (c-\frac {c}{a^2 x^2}\right )^2 \, dx=\frac {6 \, a^{3} c^{2} x^{3} \arctan \left (\sqrt {\frac {a x - 1}{a x + 1}}\right ) - 18 \, a^{3} c^{2} x^{3} \log \left (\sqrt {\frac {a x - 1}{a x + 1}} + 1\right ) + 18 \, a^{3} c^{2} x^{3} \log \left (\sqrt {\frac {a x - 1}{a x + 1}} - 1\right ) + {\left (6 \, a^{4} c^{2} x^{4} + 22 \, a^{3} c^{2} x^{3} + 7 \, a^{2} c^{2} x^{2} - 7 \, a c^{2} x + 2 \, c^{2}\right )} \sqrt {\frac {a x - 1}{a x + 1}}}{6 \, a^{4} x^{3}} \]

input
integrate((c-c/a^2/x^2)^2*((a*x-1)/(a*x+1))^(3/2),x, algorithm="fricas")
 
output
1/6*(6*a^3*c^2*x^3*arctan(sqrt((a*x - 1)/(a*x + 1))) - 18*a^3*c^2*x^3*log( 
sqrt((a*x - 1)/(a*x + 1)) + 1) + 18*a^3*c^2*x^3*log(sqrt((a*x - 1)/(a*x + 
1)) - 1) + (6*a^4*c^2*x^4 + 22*a^3*c^2*x^3 + 7*a^2*c^2*x^2 - 7*a*c^2*x + 2 
*c^2)*sqrt((a*x - 1)/(a*x + 1)))/(a^4*x^3)
 
3.9.24.6 Sympy [F]

\[ \int e^{-3 \coth ^{-1}(a x)} \left (c-\frac {c}{a^2 x^2}\right )^2 \, dx=\frac {c^{2} \left (\int \left (- \frac {\sqrt {\frac {a x}{a x + 1} - \frac {1}{a x + 1}}}{a x^{5} + x^{4}}\right )\, dx + \int \frac {a \sqrt {\frac {a x}{a x + 1} - \frac {1}{a x + 1}}}{a x^{4} + x^{3}}\, dx + \int \frac {2 a^{2} \sqrt {\frac {a x}{a x + 1} - \frac {1}{a x + 1}}}{a x^{3} + x^{2}}\, dx + \int \left (- \frac {2 a^{3} \sqrt {\frac {a x}{a x + 1} - \frac {1}{a x + 1}}}{a x^{2} + x}\right )\, dx + \int \left (- \frac {a^{4} \sqrt {\frac {a x}{a x + 1} - \frac {1}{a x + 1}}}{a x + 1}\right )\, dx + \int \frac {a^{5} x \sqrt {\frac {a x}{a x + 1} - \frac {1}{a x + 1}}}{a x + 1}\, dx\right )}{a^{4}} \]

input
integrate((c-c/a**2/x**2)**2*((a*x-1)/(a*x+1))**(3/2),x)
 
output
c**2*(Integral(-sqrt(a*x/(a*x + 1) - 1/(a*x + 1))/(a*x**5 + x**4), x) + In 
tegral(a*sqrt(a*x/(a*x + 1) - 1/(a*x + 1))/(a*x**4 + x**3), x) + Integral( 
2*a**2*sqrt(a*x/(a*x + 1) - 1/(a*x + 1))/(a*x**3 + x**2), x) + Integral(-2 
*a**3*sqrt(a*x/(a*x + 1) - 1/(a*x + 1))/(a*x**2 + x), x) + Integral(-a**4* 
sqrt(a*x/(a*x + 1) - 1/(a*x + 1))/(a*x + 1), x) + Integral(a**5*x*sqrt(a*x 
/(a*x + 1) - 1/(a*x + 1))/(a*x + 1), x))/a**4
 
3.9.24.7 Maxima [A] (verification not implemented)

Time = 0.29 (sec) , antiderivative size = 224, normalized size of antiderivative = 1.15 \[ \int e^{-3 \coth ^{-1}(a x)} \left (c-\frac {c}{a^2 x^2}\right )^2 \, dx=\frac {1}{3} \, a {\left (\frac {3 \, c^{2} \arctan \left (\sqrt {\frac {a x - 1}{a x + 1}}\right )}{a^{2}} - \frac {9 \, c^{2} \log \left (\sqrt {\frac {a x - 1}{a x + 1}} + 1\right )}{a^{2}} + \frac {9 \, c^{2} \log \left (\sqrt {\frac {a x - 1}{a x + 1}} - 1\right )}{a^{2}} - \frac {21 \, c^{2} \left (\frac {a x - 1}{a x + 1}\right )^{\frac {7}{2}} - 17 \, c^{2} \left (\frac {a x - 1}{a x + 1}\right )^{\frac {5}{2}} - 37 \, c^{2} \left (\frac {a x - 1}{a x + 1}\right )^{\frac {3}{2}} - 15 \, c^{2} \sqrt {\frac {a x - 1}{a x + 1}}}{\frac {2 \, {\left (a x - 1\right )} a^{2}}{a x + 1} - \frac {2 \, {\left (a x - 1\right )}^{3} a^{2}}{{\left (a x + 1\right )}^{3}} - \frac {{\left (a x - 1\right )}^{4} a^{2}}{{\left (a x + 1\right )}^{4}} + a^{2}}\right )} \]

input
integrate((c-c/a^2/x^2)^2*((a*x-1)/(a*x+1))^(3/2),x, algorithm="maxima")
 
output
1/3*a*(3*c^2*arctan(sqrt((a*x - 1)/(a*x + 1)))/a^2 - 9*c^2*log(sqrt((a*x - 
 1)/(a*x + 1)) + 1)/a^2 + 9*c^2*log(sqrt((a*x - 1)/(a*x + 1)) - 1)/a^2 - ( 
21*c^2*((a*x - 1)/(a*x + 1))^(7/2) - 17*c^2*((a*x - 1)/(a*x + 1))^(5/2) - 
37*c^2*((a*x - 1)/(a*x + 1))^(3/2) - 15*c^2*sqrt((a*x - 1)/(a*x + 1)))/(2* 
(a*x - 1)*a^2/(a*x + 1) - 2*(a*x - 1)^3*a^2/(a*x + 1)^3 - (a*x - 1)^4*a^2/ 
(a*x + 1)^4 + a^2))
 
3.9.24.8 Giac [A] (verification not implemented)

Time = 0.29 (sec) , antiderivative size = 264, normalized size of antiderivative = 1.35 \[ \int e^{-3 \coth ^{-1}(a x)} \left (c-\frac {c}{a^2 x^2}\right )^2 \, dx=\frac {c^{2} \arctan \left (-x {\left | a \right |} + \sqrt {a^{2} x^{2} - 1}\right ) \mathrm {sgn}\left (a x + 1\right )}{a} + \frac {3 \, c^{2} \log \left ({\left | -x {\left | a \right |} + \sqrt {a^{2} x^{2} - 1} \right |}\right ) \mathrm {sgn}\left (a x + 1\right )}{{\left | a \right |}} + \frac {\sqrt {a^{2} x^{2} - 1} c^{2} \mathrm {sgn}\left (a x + 1\right )}{a} + \frac {9 \, {\left (x {\left | a \right |} - \sqrt {a^{2} x^{2} - 1}\right )}^{5} c^{2} {\left | a \right |} \mathrm {sgn}\left (a x + 1\right ) + 12 \, {\left (x {\left | a \right |} - \sqrt {a^{2} x^{2} - 1}\right )}^{4} a c^{2} \mathrm {sgn}\left (a x + 1\right ) + 36 \, {\left (x {\left | a \right |} - \sqrt {a^{2} x^{2} - 1}\right )}^{2} a c^{2} \mathrm {sgn}\left (a x + 1\right ) - 9 \, {\left (x {\left | a \right |} - \sqrt {a^{2} x^{2} - 1}\right )} c^{2} {\left | a \right |} \mathrm {sgn}\left (a x + 1\right ) + 16 \, a c^{2} \mathrm {sgn}\left (a x + 1\right )}{3 \, {\left ({\left (x {\left | a \right |} - \sqrt {a^{2} x^{2} - 1}\right )}^{2} + 1\right )}^{3} a {\left | a \right |}} \]

input
integrate((c-c/a^2/x^2)^2*((a*x-1)/(a*x+1))^(3/2),x, algorithm="giac")
 
output
c^2*arctan(-x*abs(a) + sqrt(a^2*x^2 - 1))*sgn(a*x + 1)/a + 3*c^2*log(abs(- 
x*abs(a) + sqrt(a^2*x^2 - 1)))*sgn(a*x + 1)/abs(a) + sqrt(a^2*x^2 - 1)*c^2 
*sgn(a*x + 1)/a + 1/3*(9*(x*abs(a) - sqrt(a^2*x^2 - 1))^5*c^2*abs(a)*sgn(a 
*x + 1) + 12*(x*abs(a) - sqrt(a^2*x^2 - 1))^4*a*c^2*sgn(a*x + 1) + 36*(x*a 
bs(a) - sqrt(a^2*x^2 - 1))^2*a*c^2*sgn(a*x + 1) - 9*(x*abs(a) - sqrt(a^2*x 
^2 - 1))*c^2*abs(a)*sgn(a*x + 1) + 16*a*c^2*sgn(a*x + 1))/(((x*abs(a) - sq 
rt(a^2*x^2 - 1))^2 + 1)^3*a*abs(a))
 
3.9.24.9 Mupad [B] (verification not implemented)

Time = 0.09 (sec) , antiderivative size = 183, normalized size of antiderivative = 0.94 \[ \int e^{-3 \coth ^{-1}(a x)} \left (c-\frac {c}{a^2 x^2}\right )^2 \, dx=\frac {5\,c^2\,\sqrt {\frac {a\,x-1}{a\,x+1}}+\frac {37\,c^2\,{\left (\frac {a\,x-1}{a\,x+1}\right )}^{3/2}}{3}+\frac {17\,c^2\,{\left (\frac {a\,x-1}{a\,x+1}\right )}^{5/2}}{3}-7\,c^2\,{\left (\frac {a\,x-1}{a\,x+1}\right )}^{7/2}}{a+\frac {2\,a\,\left (a\,x-1\right )}{a\,x+1}-\frac {2\,a\,{\left (a\,x-1\right )}^3}{{\left (a\,x+1\right )}^3}-\frac {a\,{\left (a\,x-1\right )}^4}{{\left (a\,x+1\right )}^4}}+\frac {c^2\,\mathrm {atan}\left (\sqrt {\frac {a\,x-1}{a\,x+1}}\right )}{a}-\frac {6\,c^2\,\mathrm {atanh}\left (\sqrt {\frac {a\,x-1}{a\,x+1}}\right )}{a} \]

input
int((c - c/(a^2*x^2))^2*((a*x - 1)/(a*x + 1))^(3/2),x)
 
output
(5*c^2*((a*x - 1)/(a*x + 1))^(1/2) + (37*c^2*((a*x - 1)/(a*x + 1))^(3/2))/ 
3 + (17*c^2*((a*x - 1)/(a*x + 1))^(5/2))/3 - 7*c^2*((a*x - 1)/(a*x + 1))^( 
7/2))/(a + (2*a*(a*x - 1))/(a*x + 1) - (2*a*(a*x - 1)^3)/(a*x + 1)^3 - (a* 
(a*x - 1)^4)/(a*x + 1)^4) + (c^2*atan(((a*x - 1)/(a*x + 1))^(1/2)))/a - (6 
*c^2*atanh(((a*x - 1)/(a*x + 1))^(1/2)))/a