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

Optimal result

Integrand size = 25, antiderivative size = 161 \[ \int e^{-3 \coth ^{-1}(a x)} \sqrt {c-\frac {c}{a x}} x \, dx=-\frac {47 \sqrt {c-\frac {c}{a x}}}{4 a^2 \sqrt {1-\frac {1}{a^2 x^2}}}-\frac {11 \sqrt {c-\frac {c}{a x}} x}{4 a \sqrt {1-\frac {1}{a^2 x^2}}}+\frac {\left (c-\frac {c}{a x}\right )^{3/2} x^2}{2 c \sqrt {1-\frac {1}{a^2 x^2}}}+\frac {47 \sqrt {c} \text {arctanh}\left (\frac {\sqrt {c} \sqrt {1-\frac {1}{a^2 x^2}}}{\sqrt {c-\frac {c}{a x}}}\right )}{4 a^2} \] Output:

-47/4*(c-c/a/x)^(1/2)/a^2/(1-1/a^2/x^2)^(1/2)-11/4*(c-c/a/x)^(1/2)*x/a/(1- 
1/a^2/x^2)^(1/2)+1/2*(c-c/a/x)^(3/2)*x^2/c/(1-1/a^2/x^2)^(1/2)+47/4*c^(1/2 
)*arctanh(c^(1/2)*(1-1/a^2/x^2)^(1/2)/(c-c/a/x)^(1/2))/a^2
 

Mathematica [A] (verified)

Time = 2.08 (sec) , antiderivative size = 151, normalized size of antiderivative = 0.94 \[ \int e^{-3 \coth ^{-1}(a x)} \sqrt {c-\frac {c}{a x}} x \, dx=\frac {\sqrt {1-\frac {1}{a^2 x^2}} \sqrt {c-\frac {c}{a x}} x^2 \left (-47-13 a x+2 a^2 x^2\right )}{-4+4 a^2 x^2}-\frac {47 \sqrt {c} \log (1-a x)}{8 a^2}+\frac {47 \sqrt {c} \log \left (2 a^2 \sqrt {c} \sqrt {1-\frac {1}{a^2 x^2}} \sqrt {c-\frac {c}{a x}} x^2+c \left (-1-a x+2 a^2 x^2\right )\right )}{8 a^2} \] Input:

Integrate[(Sqrt[c - c/(a*x)]*x)/E^(3*ArcCoth[a*x]),x]
 

Output:

(Sqrt[1 - 1/(a^2*x^2)]*Sqrt[c - c/(a*x)]*x^2*(-47 - 13*a*x + 2*a^2*x^2))/( 
-4 + 4*a^2*x^2) - (47*Sqrt[c]*Log[1 - a*x])/(8*a^2) + (47*Sqrt[c]*Log[2*a^ 
2*Sqrt[c]*Sqrt[1 - 1/(a^2*x^2)]*Sqrt[c - c/(a*x)]*x^2 + c*(-1 - a*x + 2*a^ 
2*x^2)])/(8*a^2)
 

Rubi [A] (verified)

Time = 0.64 (sec) , antiderivative size = 116, normalized size of antiderivative = 0.72, number of steps used = 10, number of rules used = 9, \(\frac {\text {number of rules}}{\text {integrand size}}\) = 0.360, Rules used = {6733, 585, 27, 100, 27, 87, 61, 73, 221}

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

\(\Big \downarrow \) 6733

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

\(\Big \downarrow \) 585

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

\(\Big \downarrow \) 27

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

\(\Big \downarrow \) 100

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

\(\Big \downarrow \) 27

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

\(\Big \downarrow \) 87

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

\(\Big \downarrow \) 61

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

\(\Big \downarrow \) 73

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

\(\Big \downarrow \) 221

\(\displaystyle -\frac {\left (\frac {1}{4} \left (\frac {47}{2} \left (\frac {2}{\sqrt {\frac {1}{a x}+1}}-2 \text {arctanh}\left (\sqrt {\frac {1}{a x}+1}\right )\right )+\frac {13 a x}{\sqrt {\frac {1}{a x}+1}}\right )-\frac {a^2 x^2}{2 \sqrt {\frac {1}{a x}+1}}\right ) \sqrt {c-\frac {c}{a x}}}{a^2 \sqrt {1-\frac {1}{a x}}}\)

Input:

Int[(Sqrt[c - c/(a*x)]*x)/E^(3*ArcCoth[a*x]),x]
 

Output:

-((Sqrt[c - c/(a*x)]*(-1/2*(a^2*x^2)/Sqrt[1 + 1/(a*x)] + ((13*a*x)/Sqrt[1 
+ 1/(a*x)] + (47*(2/Sqrt[1 + 1/(a*x)] - 2*ArcTanh[Sqrt[1 + 1/(a*x)]]))/2)/ 
4))/(a^2*Sqrt[1 - 1/(a*x)]))
 

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 61

Int[((a_.) + (b_.)*(x_))^(m_)*((c_.) + (d_.)*(x_))^(n_), x_Symbol] :> Simp[ 
(a + b*x)^(m + 1)*((c + d*x)^(n + 1)/((b*c - a*d)*(m + 1))), x] - Simp[d*(( 
m + n + 2)/((b*c - a*d)*(m + 1)))   Int[(a + b*x)^(m + 1)*(c + d*x)^n, x], 
x] /; FreeQ[{a, b, c, d, n}, x] && LtQ[m, -1] &&  !(LtQ[n, -1] && (EqQ[a, 0 
] || (NeQ[c, 0] && LtQ[m - n, 0] && IntegerQ[n]))) && IntLinearQ[a, b, c, d 
, m, n, x]
 

rule 73

Int[((a_.) + (b_.)*(x_))^(m_)*((c_.) + (d_.)*(x_))^(n_), x_Symbol] :> With[ 
{p = Denominator[m]}, Simp[p/b   Subst[Int[x^(p*(m + 1) - 1)*(c - a*(d/b) + 
 d*(x^p/b))^n, x], x, (a + b*x)^(1/p)], x]] /; FreeQ[{a, b, c, d}, x] && Lt 
Q[-1, m, 0] && LeQ[-1, n, 0] && LeQ[Denominator[n], Denominator[m]] && IntL 
inearQ[a, b, c, d, m, n, x]
 

rule 87

Int[((a_.) + (b_.)*(x_))*((c_.) + (d_.)*(x_))^(n_.)*((e_.) + (f_.)*(x_))^(p 
_.), x_] :> Simp[(-(b*e - a*f))*(c + d*x)^(n + 1)*((e + f*x)^(p + 1)/(f*(p 
+ 1)*(c*f - d*e))), x] - Simp[(a*d*f*(n + p + 2) - b*(d*e*(n + 1) + c*f*(p 
+ 1)))/(f*(p + 1)*(c*f - d*e))   Int[(c + d*x)^n*(e + f*x)^(p + 1), x], x] 
/; FreeQ[{a, b, c, d, e, f, n}, x] && LtQ[p, -1] && ( !LtQ[n, -1] || Intege 
rQ[p] ||  !(IntegerQ[n] ||  !(EqQ[e, 0] ||  !(EqQ[c, 0] || LtQ[p, n]))))
 

rule 100

Int[((a_.) + (b_.)*(x_))^2*((c_.) + (d_.)*(x_))^(n_)*((e_.) + (f_.)*(x_))^( 
p_), x_] :> Simp[(b*c - a*d)^2*(c + d*x)^(n + 1)*((e + f*x)^(p + 1)/(d^2*(d 
*e - c*f)*(n + 1))), x] - Simp[1/(d^2*(d*e - c*f)*(n + 1))   Int[(c + d*x)^ 
(n + 1)*(e + f*x)^p*Simp[a^2*d^2*f*(n + p + 2) + b^2*c*(d*e*(n + 1) + c*f*( 
p + 1)) - 2*a*b*d*(d*e*(n + 1) + c*f*(p + 1)) - b^2*d*(d*e - c*f)*(n + 1)*x 
, x], x], x] /; FreeQ[{a, b, c, d, e, f, n, p}, x] && (LtQ[n, -1] || (EqQ[n 
 + p + 3, 0] && NeQ[n, -1] && (SumSimplerQ[n, 1] ||  !SumSimplerQ[p, 1])))
 

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 585

Int[((e_.)*(x_))^(m_.)*((c_) + (d_.)*(x_))^(n_.)*((a_) + (b_.)*(x_)^2)^(p_) 
, x_Symbol] :> Simp[a^p*c^IntPart[n]*((c + d*x)^FracPart[n]/(1 + d*(x/c))^F 
racPart[n])   Int[(e*x)^m*(1 - d*(x/c))^p*(1 + d*(x/c))^(n + p), x], x] /; 
FreeQ[{a, b, c, d, e, m, n, p}, x] && EqQ[b*c^2 + a*d^2, 0] && GtQ[a, 0]
 

rule 6733

Int[E^(ArcCoth[(a_.)*(x_)]*(n_.))*((c_) + (d_.)/(x_))^(p_.)*(x_)^(m_.), x_S 
ymbol] :> Simp[-c^n   Subst[Int[(c + d*x)^(p - n)*((1 - x^2/a^2)^(n/2)/x^(m 
 + 2)), x], x, 1/x], x] /; FreeQ[{a, c, d, p}, x] && EqQ[c + a*d, 0] && Int 
egerQ[(n - 1)/2] && IntegerQ[m] && IntegerQ[2*p]
 

Maple [A] (verified)

Time = 0.10 (sec) , antiderivative size = 163, normalized size of antiderivative = 1.01

method	result	size
default	\(\frac {\left (\frac {a x -1}{a x +1}\right )^{\frac {3}{2}} \left (a x +1\right ) \sqrt {\frac {c \left (a x -1\right )}{a x}}\, x \left (4 a^{\frac {5}{2}} x^{2} \sqrt {x \left (a x +1\right )}-26 a^{\frac {3}{2}} x \sqrt {x \left (a x +1\right )}+47 \ln \left (\frac {2 \sqrt {x \left (a x +1\right )}\, \sqrt {a}+2 a x +1}{2 \sqrt {a}}\right ) a x -94 \sqrt {x \left (a x +1\right )}\, \sqrt {a}+47 \ln \left (\frac {2 \sqrt {x \left (a x +1\right )}\, \sqrt {a}+2 a x +1}{2 \sqrt {a}}\right )\right )}{8 \left (a x -1\right )^{2} a^{\frac {3}{2}} \sqrt {x \left (a x +1\right )}}\)	\(163\)
risch	\(\frac {\left (2 a x -15\right ) \left (a x +1\right ) x \sqrt {\frac {a x -1}{a x +1}}\, \sqrt {\frac {c \left (a x -1\right )}{a x}}}{4 a \left (a x -1\right )}+\frac {\left (\frac {47 \ln \left (\frac {\frac {1}{2} a c +a^{2} c x}{\sqrt {a^{2} c}}+\sqrt {a^{2} c \,x^{2}+a c x}\right )}{8 a \sqrt {a^{2} c}}-\frac {8 \sqrt {\left (x +\frac {1}{a}\right )^{2} a^{2} c -\left (x +\frac {1}{a}\right ) a c}}{a^{3} c \left (x +\frac {1}{a}\right )}\right ) \sqrt {\frac {a x -1}{a x +1}}\, \sqrt {\frac {c \left (a x -1\right )}{a x}}\, \sqrt {\left (a x +1\right ) a c x}}{a x -1}\)	\(193\)

Input:

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

Output:

1/8*((a*x-1)/(a*x+1))^(3/2)/(a*x-1)^2*(a*x+1)*(c*(a*x-1)/a/x)^(1/2)*x*(4*a 
^(5/2)*x^2*(x*(a*x+1))^(1/2)-26*a^(3/2)*x*(x*(a*x+1))^(1/2)+47*ln(1/2*(2*( 
x*(a*x+1))^(1/2)*a^(1/2)+2*a*x+1)/a^(1/2))*a*x-94*(x*(a*x+1))^(1/2)*a^(1/2 
)+47*ln(1/2*(2*(x*(a*x+1))^(1/2)*a^(1/2)+2*a*x+1)/a^(1/2)))/a^(3/2)/(x*(a* 
x+1))^(1/2)
 

Fricas [A] (verification not implemented)

Time = 0.16 (sec) , antiderivative size = 321, normalized size of antiderivative = 1.99 \[ \int e^{-3 \coth ^{-1}(a x)} \sqrt {c-\frac {c}{a x}} x \, dx=\left [\frac {47 \, {\left (a x - 1\right )} \sqrt {c} \log \left (-\frac {8 \, a^{3} c x^{3} - 7 \, a c x + 4 \, {\left (2 \, a^{3} x^{3} + 3 \, a^{2} x^{2} + a x\right )} \sqrt {c} \sqrt {\frac {a x - 1}{a x + 1}} \sqrt {\frac {a c x - c}{a x}} - c}{a x - 1}\right ) + 4 \, {\left (2 \, a^{3} x^{3} - 13 \, a^{2} x^{2} - 47 \, a x\right )} \sqrt {\frac {a x - 1}{a x + 1}} \sqrt {\frac {a c x - c}{a x}}}{16 \, {\left (a^{3} x - a^{2}\right )}}, -\frac {47 \, {\left (a x - 1\right )} \sqrt {-c} \arctan \left (\frac {2 \, {\left (a^{2} x^{2} + a x\right )} \sqrt {-c} \sqrt {\frac {a x - 1}{a x + 1}} \sqrt {\frac {a c x - c}{a x}}}{2 \, a^{2} c x^{2} - a c x - c}\right ) - 2 \, {\left (2 \, a^{3} x^{3} - 13 \, a^{2} x^{2} - 47 \, a x\right )} \sqrt {\frac {a x - 1}{a x + 1}} \sqrt {\frac {a c x - c}{a x}}}{8 \, {\left (a^{3} x - a^{2}\right )}}\right ] \] Input:

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

Output:

[1/16*(47*(a*x - 1)*sqrt(c)*log(-(8*a^3*c*x^3 - 7*a*c*x + 4*(2*a^3*x^3 + 3 
*a^2*x^2 + a*x)*sqrt(c)*sqrt((a*x - 1)/(a*x + 1))*sqrt((a*c*x - c)/(a*x)) 
- c)/(a*x - 1)) + 4*(2*a^3*x^3 - 13*a^2*x^2 - 47*a*x)*sqrt((a*x - 1)/(a*x 
+ 1))*sqrt((a*c*x - c)/(a*x)))/(a^3*x - a^2), -1/8*(47*(a*x - 1)*sqrt(-c)* 
arctan(2*(a^2*x^2 + a*x)*sqrt(-c)*sqrt((a*x - 1)/(a*x + 1))*sqrt((a*c*x - 
c)/(a*x))/(2*a^2*c*x^2 - a*c*x - c)) - 2*(2*a^3*x^3 - 13*a^2*x^2 - 47*a*x) 
*sqrt((a*x - 1)/(a*x + 1))*sqrt((a*c*x - c)/(a*x)))/(a^3*x - a^2)]
 

Sympy [F(-1)]

Timed out. \[ \int e^{-3 \coth ^{-1}(a x)} \sqrt {c-\frac {c}{a x}} x \, dx=\text {Timed out} \] Input:

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

Output:

Timed out
 

Maxima [F]

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

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

Output:

integrate(sqrt(c - c/(a*x))*x*((a*x - 1)/(a*x + 1))^(3/2), x)
 

Giac [F(-2)]

Exception generated. \[ \int e^{-3 \coth ^{-1}(a x)} \sqrt {c-\frac {c}{a x}} x \, dx=\text {Exception raised: TypeError} \] Input:

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

Output:

Exception raised: TypeError >> an error occurred running a Giac command:IN 
PUT:sage2:=int(sage0,sageVARx):;OUTPUT:sym2poly/r2sym(const gen & e,const 
index_m & i,const vecteur & l) Error: Bad Argument Value
 

Mupad [F(-1)]

Timed out. \[ \int e^{-3 \coth ^{-1}(a x)} \sqrt {c-\frac {c}{a x}} x \, dx=\int x\,\sqrt {c-\frac {c}{a\,x}}\,{\left (\frac {a\,x-1}{a\,x+1}\right )}^{3/2} \,d x \] Input:

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

Output:

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

Reduce [B] (verification not implemented)

Time = 0.15 (sec) , antiderivative size = 71, normalized size of antiderivative = 0.44 \[ \int e^{-3 \coth ^{-1}(a x)} \sqrt {c-\frac {c}{a x}} x \, dx=\frac {\sqrt {c}\, \left (47 \sqrt {a x +1}\, \mathrm {log}\left (\sqrt {a x +1}+\sqrt {x}\, \sqrt {a}\right )-36 \sqrt {a x +1}+2 \sqrt {x}\, \sqrt {a}\, a^{2} x^{2}-13 \sqrt {x}\, \sqrt {a}\, a x -47 \sqrt {x}\, \sqrt {a}\right )}{4 \sqrt {a x +1}\, a^{2}} \] Input:

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

Output:

(sqrt(c)*(47*sqrt(a*x + 1)*log(sqrt(a*x + 1) + sqrt(x)*sqrt(a)) - 36*sqrt( 
a*x + 1) + 2*sqrt(x)*sqrt(a)*a**2*x**2 - 13*sqrt(x)*sqrt(a)*a*x - 47*sqrt( 
x)*sqrt(a)))/(4*sqrt(a*x + 1)*a**2)