\(\int e^{-3 \text {arctanh}(a x)} (e x)^m \sqrt {c+a c x} \, dx\) [474]

Optimal result

Integrand size = 24, antiderivative size = 54 \[ \int e^{-3 \text {arctanh}(a x)} (e x)^m \sqrt {c+a c x} \, dx=\frac {c (e x)^{1+m} \sqrt {1+a x} \operatorname {AppellF1}\left (1+m,-\frac {3}{2},1,2+m,a x,-a x\right )}{e (1+m) \sqrt {c+a c x}} \] Output:

c*(e*x)^(1+m)*(a*x+1)^(1/2)*AppellF1(1+m,-3/2,1,2+m,a*x,-a*x)/e/(1+m)/(a*c 
*x+c)^(1/2)
 

Mathematica [F]

\[ \int e^{-3 \text {arctanh}(a x)} (e x)^m \sqrt {c+a c x} \, dx=\int e^{-3 \text {arctanh}(a x)} (e x)^m \sqrt {c+a c x} \, dx \] Input:

Integrate[((e*x)^m*Sqrt[c + a*c*x])/E^(3*ArcTanh[a*x]),x]
 

Output:

Integrate[((e*x)^m*Sqrt[c + a*c*x])/E^(3*ArcTanh[a*x]), x]
 

Rubi [A] (verified)

Time = 0.48 (sec) , antiderivative size = 53, normalized size of antiderivative = 0.98, number of steps used = 3, number of rules used = 3, \(\frac {\text {number of rules}}{\text {integrand size}}\) = 0.125, Rules used = {6680, 37, 150}

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 e^{-3 \text {arctanh}(a x)} \sqrt {a c x+c} (e x)^m \, dx\)

\(\Big \downarrow \) 6680

\(\displaystyle \int \frac {(1-a x)^{3/2} \sqrt {a c x+c} (e x)^m}{(a x+1)^{3/2}}dx\)

\(\Big \downarrow \) 37

\(\displaystyle \frac {\sqrt {a c x+c} \int \frac {(e x)^m (1-a x)^{3/2}}{a x+1}dx}{\sqrt {a x+1}}\)

\(\Big \downarrow \) 150

\(\displaystyle \frac {\sqrt {a c x+c} (e x)^{m+1} \operatorname {AppellF1}\left (m+1,-\frac {3}{2},1,m+2,a x,-a x\right )}{e (m+1) \sqrt {a x+1}}\)

Input:

Int[((e*x)^m*Sqrt[c + a*c*x])/E^(3*ArcTanh[a*x]),x]
 

Output:

((e*x)^(1 + m)*Sqrt[c + a*c*x]*AppellF1[1 + m, -3/2, 1, 2 + m, a*x, -(a*x) 
])/(e*(1 + m)*Sqrt[1 + a*x])
 

Defintions of rubi rules used

rule 37

Int[(u_.)*((a_) + (b_.)*(x_))^(m_)*((c_) + (d_.)*(x_))^(n_), x_Symbol] :> S 
imp[(a + b*x)^m/(c + d*x)^m   Int[u*(c + d*x)^(m + n), x], x] /; FreeQ[{a, 
b, c, d, m, n}, x] && EqQ[b*c - a*d, 0] &&  !SimplerQ[a + b*x, c + d*x]
 

rule 150

Int[((b_.)*(x_))^(m_)*((c_) + (d_.)*(x_))^(n_)*((e_) + (f_.)*(x_))^(p_), x_ 
] :> Simp[c^n*e^p*((b*x)^(m + 1)/(b*(m + 1)))*AppellF1[m + 1, -n, -p, m + 2 
, (-d)*(x/c), (-f)*(x/e)], x] /; FreeQ[{b, c, d, e, f, m, n, p}, x] &&  !In 
tegerQ[m] &&  !IntegerQ[n] && GtQ[c, 0] && (IntegerQ[p] || GtQ[e, 0])
 

rule 6680

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

Maple [F]

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

Input:

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

Output:

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

Fricas [F]

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

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

Output:

integral(-sqrt(-a^2*x^2 + 1)*sqrt(a*c*x + c)*(a*x - 1)*(e*x)^m/(a^2*x^2 + 
2*a*x + 1), x)
 

Sympy [F(-1)]

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

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

Output:

Timed out
 

Maxima [F]

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

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

Output:

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

Giac [F(-2)]

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

integrate((e*x)^m*(a*c*x+c)^(1/2)/(a*x+1)^3*(-a^2*x^2+1)^(3/2),x, algorith 
m="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 \text {arctanh}(a x)} (e x)^m \sqrt {c+a c x} \, dx=\int \frac {{\left (e\,x\right )}^m\,{\left (1-a^2\,x^2\right )}^{3/2}\,\sqrt {c+a\,c\,x}}{{\left (a\,x+1\right )}^3} \,d x \] Input:

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

Output:

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

Reduce [F]

\[ \int e^{-3 \text {arctanh}(a x)} (e x)^m \sqrt {c+a c x} \, dx=\frac {2 e^{m} \sqrt {c}\, \left (-x^{m} \sqrt {-a x +1}\, a x -4 x^{m} \sqrt {-a x +1}\, m -5 x^{m} \sqrt {-a x +1}+8 \left (\int \frac {x^{m} \sqrt {-a x +1}\, x}{2 a^{2} m \,x^{2}+3 a^{2} x^{2}-2 m -3}d x \right ) a^{2} m^{3}+30 \left (\int \frac {x^{m} \sqrt {-a x +1}\, x}{2 a^{2} m \,x^{2}+3 a^{2} x^{2}-2 m -3}d x \right ) a^{2} m^{2}+39 \left (\int \frac {x^{m} \sqrt {-a x +1}\, x}{2 a^{2} m \,x^{2}+3 a^{2} x^{2}-2 m -3}d x \right ) a^{2} m +18 \left (\int \frac {x^{m} \sqrt {-a x +1}\, x}{2 a^{2} m \,x^{2}+3 a^{2} x^{2}-2 m -3}d x \right ) a^{2}-8 \left (\int \frac {x^{m} \sqrt {-a x +1}}{2 a^{2} m \,x^{3}+3 a^{2} x^{3}-2 m x -3 x}d x \right ) m^{3}-22 \left (\int \frac {x^{m} \sqrt {-a x +1}}{2 a^{2} m \,x^{3}+3 a^{2} x^{3}-2 m x -3 x}d x \right ) m^{2}-15 \left (\int \frac {x^{m} \sqrt {-a x +1}}{2 a^{2} m \,x^{3}+3 a^{2} x^{3}-2 m x -3 x}d x \right ) m \right )}{a \left (2 m +3\right )} \] Input:

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

Output:

(2*e**m*sqrt(c)*( - x**m*sqrt( - a*x + 1)*a*x - 4*x**m*sqrt( - a*x + 1)*m 
- 5*x**m*sqrt( - a*x + 1) + 8*int((x**m*sqrt( - a*x + 1)*x)/(2*a**2*m*x**2 
 + 3*a**2*x**2 - 2*m - 3),x)*a**2*m**3 + 30*int((x**m*sqrt( - a*x + 1)*x)/ 
(2*a**2*m*x**2 + 3*a**2*x**2 - 2*m - 3),x)*a**2*m**2 + 39*int((x**m*sqrt( 
- a*x + 1)*x)/(2*a**2*m*x**2 + 3*a**2*x**2 - 2*m - 3),x)*a**2*m + 18*int(( 
x**m*sqrt( - a*x + 1)*x)/(2*a**2*m*x**2 + 3*a**2*x**2 - 2*m - 3),x)*a**2 - 
 8*int((x**m*sqrt( - a*x + 1))/(2*a**2*m*x**3 + 3*a**2*x**3 - 2*m*x - 3*x) 
,x)*m**3 - 22*int((x**m*sqrt( - a*x + 1))/(2*a**2*m*x**3 + 3*a**2*x**3 - 2 
*m*x - 3*x),x)*m**2 - 15*int((x**m*sqrt( - a*x + 1))/(2*a**2*m*x**3 + 3*a* 
*2*x**3 - 2*m*x - 3*x),x)*m))/(a*(2*m + 3))