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

Optimal result

Integrand size = 25, antiderivative size = 194 \[ \int \frac {e^{3 \text {arctanh}(a x)} \left (c-a^2 c x^2\right )^p}{x^3} \, dx=-\frac {\left (c-a^2 c x^2\right )^p}{2 x^2 \sqrt {1-a^2 x^2}}-\frac {3 a \left (c-a^2 c x^2\right )^p}{x \sqrt {1-a^2 x^2}}+a^3 (7-6 p) x \left (1-a^2 x^2\right )^{-p} \left (c-a^2 c x^2\right )^p \operatorname {Hypergeometric2F1}\left (\frac {1}{2},\frac {3}{2}-p,\frac {3}{2},a^2 x^2\right )+\frac {a^2 (9-2 p) \left (c-a^2 c x^2\right )^p \operatorname {Hypergeometric2F1}\left (1,-\frac {1}{2}+p,\frac {1}{2}+p,1-a^2 x^2\right )}{2 (1-2 p) \sqrt {1-a^2 x^2}} \] Output:

-1/2*(-a^2*c*x^2+c)^p/x^2/(-a^2*x^2+1)^(1/2)-3*a*(-a^2*c*x^2+c)^p/x/(-a^2* 
x^2+1)^(1/2)+a^3*(7-6*p)*x*(-a^2*c*x^2+c)^p*hypergeom([1/2, 3/2-p],[3/2],a 
^2*x^2)/((-a^2*x^2+1)^p)+1/2*a^2*(9-2*p)*(-a^2*c*x^2+c)^p*hypergeom([1, -1 
/2+p],[1/2+p],-a^2*x^2+1)/(1-2*p)/(-a^2*x^2+1)^(1/2)
 

Mathematica [A] (verified)

Time = 0.29 (sec) , antiderivative size = 154, normalized size of antiderivative = 0.79 \[ \int \frac {e^{3 \text {arctanh}(a x)} \left (c-a^2 c x^2\right )^p}{x^3} \, dx=a \left (1-a^2 x^2\right )^{-p} \left (c-a^2 c x^2\right )^p \left (-\frac {3 \operatorname {Hypergeometric2F1}\left (-\frac {1}{2},\frac {3}{2}-p,\frac {1}{2},a^2 x^2\right )}{x}+a \left (a x \operatorname {Hypergeometric2F1}\left (\frac {1}{2},\frac {3}{2}-p,\frac {3}{2},a^2 x^2\right )+\frac {\left (1-a^2 x^2\right )^{-\frac {1}{2}+p} \left (3 \operatorname {Hypergeometric2F1}\left (1,-\frac {1}{2}+p,\frac {1}{2}+p,1-a^2 x^2\right )+\operatorname {Hypergeometric2F1}\left (2,-\frac {1}{2}+p,\frac {1}{2}+p,1-a^2 x^2\right )\right )}{1-2 p}\right )\right ) \] Input:

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

Output:

(a*(c - a^2*c*x^2)^p*((-3*Hypergeometric2F1[-1/2, 3/2 - p, 1/2, a^2*x^2])/ 
x + a*(a*x*Hypergeometric2F1[1/2, 3/2 - p, 3/2, a^2*x^2] + ((1 - a^2*x^2)^ 
(-1/2 + p)*(3*Hypergeometric2F1[1, -1/2 + p, 1/2 + p, 1 - a^2*x^2] + Hyper 
geometric2F1[2, -1/2 + p, 1/2 + p, 1 - a^2*x^2]))/(1 - 2*p))))/(1 - a^2*x^ 
2)^p
 

Rubi [A] (verified)

Time = 0.58 (sec) , antiderivative size = 162, normalized size of antiderivative = 0.84, number of steps used = 9, number of rules used = 8, \(\frac {\text {number of rules}}{\text {integrand size}}\) = 0.320, Rules used = {6703, 6698, 543, 354, 87, 75, 359, 237}

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^{3 \text {arctanh}(a x)} \left (c-a^2 c x^2\right )^p}{x^3} \, dx\)

\(\Big \downarrow \) 6703

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

\(\Big \downarrow \) 6698

\(\displaystyle \left (1-a^2 x^2\right )^{-p} \left (c-a^2 c x^2\right )^p \int \frac {(a x+1)^3 \left (1-a^2 x^2\right )^{p-\frac {3}{2}}}{x^3}dx\)

\(\Big \downarrow \) 543

\(\displaystyle \left (1-a^2 x^2\right )^{-p} \left (c-a^2 c x^2\right )^p \left (\int \frac {\left (1-a^2 x^2\right )^{p-\frac {3}{2}} \left (3 a^2 x^2+1\right )}{x^3}dx+\int \frac {\left (1-a^2 x^2\right )^{p-\frac {3}{2}} \left (x^2 a^3+3 a\right )}{x^2}dx\right )\)

\(\Big \downarrow \) 354

\(\displaystyle \left (1-a^2 x^2\right )^{-p} \left (c-a^2 c x^2\right )^p \left (\frac {1}{2} \int \frac {\left (1-a^2 x^2\right )^{p-\frac {3}{2}} \left (3 a^2 x^2+1\right )}{x^4}dx^2+\int \frac {\left (1-a^2 x^2\right )^{p-\frac {3}{2}} \left (x^2 a^3+3 a\right )}{x^2}dx\right )\)

\(\Big \downarrow \) 87

\(\displaystyle \left (1-a^2 x^2\right )^{-p} \left (c-a^2 c x^2\right )^p \left (\frac {1}{2} \left (\frac {1}{2} a^2 (9-2 p) \int \frac {\left (1-a^2 x^2\right )^{p-\frac {3}{2}}}{x^2}dx^2-\frac {\left (1-a^2 x^2\right )^{p-\frac {1}{2}}}{x^2}\right )+\int \frac {\left (1-a^2 x^2\right )^{p-\frac {3}{2}} \left (x^2 a^3+3 a\right )}{x^2}dx\right )\)

\(\Big \downarrow \) 75

\(\displaystyle \left (1-a^2 x^2\right )^{-p} \left (c-a^2 c x^2\right )^p \left (\int \frac {\left (1-a^2 x^2\right )^{p-\frac {3}{2}} \left (x^2 a^3+3 a\right )}{x^2}dx+\frac {1}{2} \left (\frac {a^2 (9-2 p) \left (1-a^2 x^2\right )^{p-\frac {1}{2}} \operatorname {Hypergeometric2F1}\left (1,p-\frac {1}{2},p+\frac {1}{2},1-a^2 x^2\right )}{1-2 p}-\frac {\left (1-a^2 x^2\right )^{p-\frac {1}{2}}}{x^2}\right )\right )\)

\(\Big \downarrow \) 359

\(\displaystyle \left (1-a^2 x^2\right )^{-p} \left (c-a^2 c x^2\right )^p \left (a^3 (7-6 p) \int \left (1-a^2 x^2\right )^{p-\frac {3}{2}}dx+\frac {1}{2} \left (\frac {a^2 (9-2 p) \left (1-a^2 x^2\right )^{p-\frac {1}{2}} \operatorname {Hypergeometric2F1}\left (1,p-\frac {1}{2},p+\frac {1}{2},1-a^2 x^2\right )}{1-2 p}-\frac {\left (1-a^2 x^2\right )^{p-\frac {1}{2}}}{x^2}\right )-\frac {3 a \left (1-a^2 x^2\right )^{p-\frac {1}{2}}}{x}\right )\)

\(\Big \downarrow \) 237

\(\displaystyle \left (1-a^2 x^2\right )^{-p} \left (c-a^2 c x^2\right )^p \left (\frac {1}{2} \left (\frac {a^2 (9-2 p) \left (1-a^2 x^2\right )^{p-\frac {1}{2}} \operatorname {Hypergeometric2F1}\left (1,p-\frac {1}{2},p+\frac {1}{2},1-a^2 x^2\right )}{1-2 p}-\frac {\left (1-a^2 x^2\right )^{p-\frac {1}{2}}}{x^2}\right )-\frac {3 a \left (1-a^2 x^2\right )^{p-\frac {1}{2}}}{x}+a^3 (7-6 p) x \operatorname {Hypergeometric2F1}\left (\frac {1}{2},\frac {3}{2}-p,\frac {3}{2},a^2 x^2\right )\right )\)

Input:

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

Output:

((c - a^2*c*x^2)^p*((-3*a*(1 - a^2*x^2)^(-1/2 + p))/x + a^3*(7 - 6*p)*x*Hy 
pergeometric2F1[1/2, 3/2 - p, 3/2, a^2*x^2] + (-((1 - a^2*x^2)^(-1/2 + p)/ 
x^2) + (a^2*(9 - 2*p)*(1 - a^2*x^2)^(-1/2 + p)*Hypergeometric2F1[1, -1/2 + 
 p, 1/2 + p, 1 - a^2*x^2])/(1 - 2*p))/2))/(1 - a^2*x^2)^p
 

Defintions of rubi rules used

rule 75

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

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 237

Int[((a_) + (b_.)*(x_)^2)^(p_), x_Symbol] :> Simp[a^p*x*Hypergeometric2F1[- 
p, 1/2, 1/2 + 1, (-b)*(x^2/a)], x] /; FreeQ[{a, b, p}, x] &&  !IntegerQ[2*p 
] && GtQ[a, 0]
 

rule 354

Int[(x_)^(m_.)*((a_) + (b_.)*(x_)^2)^(p_.)*((c_) + (d_.)*(x_)^2)^(q_.), x_S 
ymbol] :> Simp[1/2   Subst[Int[x^((m - 1)/2)*(a + b*x)^p*(c + d*x)^q, x], x 
, x^2], x] /; FreeQ[{a, b, c, d, p, q}, x] && NeQ[b*c - a*d, 0] && IntegerQ 
[(m - 1)/2]
 

rule 359

Int[((e_.)*(x_))^(m_.)*((a_) + (b_.)*(x_)^2)^(p_.)*((c_) + (d_.)*(x_)^2), x 
_Symbol] :> Simp[c*(e*x)^(m + 1)*((a + b*x^2)^(p + 1)/(a*e*(m + 1))), x] + 
Simp[(a*d*(m + 1) - b*c*(m + 2*p + 3))/(a*e^2*(m + 1))   Int[(e*x)^(m + 2)* 
(a + b*x^2)^p, x], x] /; FreeQ[{a, b, c, d, e, p}, x] && NeQ[b*c - a*d, 0] 
&& LtQ[m, -1] &&  !ILtQ[p, -1]
 

rule 543

Int[(x_)^(m_.)*((c_) + (d_.)*(x_))^(n_)*((a_) + (b_.)*(x_)^2)^(p_), x_Symbo 
l] :> Module[{k}, Int[x^m*Sum[Binomial[n, 2*k]*c^(n - 2*k)*d^(2*k)*x^(2*k), 
 {k, 0, n/2}]*(a + b*x^2)^p, x] + Int[x^(m + 1)*Sum[Binomial[n, 2*k + 1]*c^ 
(n - 2*k - 1)*d^(2*k + 1)*x^(2*k), {k, 0, (n - 1)/2}]*(a + b*x^2)^p, x]] /; 
 FreeQ[{a, b, c, d, p}, x] && IGtQ[n, 1] && IntegerQ[m] &&  !IntegerQ[2*p] 
&&  !(EqQ[m, 1] && EqQ[b*c^2 + a*d^2, 0])
 

rule 6698

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

rule 6703

Int[E^(ArcTanh[(a_.)*(x_)]*(n_.))*(x_)^(m_.)*((c_) + (d_.)*(x_)^2)^(p_), x_ 
Symbol] :> Simp[c^IntPart[p]*((c + d*x^2)^FracPart[p]/(1 - a^2*x^2)^FracPar 
t[p])   Int[x^m*(1 - a^2*x^2)^p*E^(n*ArcTanh[a*x]), x], x] /; FreeQ[{a, c, 
d, m, n, p}, x] && EqQ[a^2*c + d, 0] &&  !(IntegerQ[p] || GtQ[c, 0]) &&  !I 
ntegerQ[n/2]
 

Maple [F]

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

Input:

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

Output:

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

Fricas [F]

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

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

Output:

integral(sqrt(-a^2*x^2 + 1)*(a*x + 1)*(-a^2*c*x^2 + c)^p/(a^2*x^5 - 2*a*x^ 
4 + x^3), x)
 

Sympy [F]

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

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

Output:

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

Maxima [F]

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

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

Output:

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

Giac [F]

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

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

Output:

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

Mupad [F(-1)]

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

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

Output:

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

Reduce [F]

\[ \int \frac {e^{3 \text {arctanh}(a x)} \left (c-a^2 c x^2\right )^p}{x^3} \, dx =\text {Too large to display} \] Input:

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

Output:

(4*( - a**2*c*x**2 + c)**p*sqrt( - a**2*x**2 + 1)*a**2*p**2*x**2 - 6*( - a 
**2*c*x**2 + c)**p*sqrt( - a**2*x**2 + 1)*a**2*p*x**2 + 2*( - a**2*c*x**2 
+ c)**p*sqrt( - a**2*x**2 + 1)*a**2*x**2 + 2*( - a**2*c*x**2 + c)**p*sqrt( 
 - a**2*x**2 + 1)*a*p*x - 2*( - a**2*c*x**2 + c)**p*sqrt( - a**2*x**2 + 1) 
*a*x + ( - a**2*c*x**2 + c)**p*sqrt( - a**2*x**2 + 1) - 4*int(( - a**2*c*x 
**2 + c)**p/(sqrt( - a**2*x**2 + 1)*a*x**3 - sqrt( - a**2*x**2 + 1)*x**2), 
x)*a**2*x**3 + 4*int(( - a**2*c*x**2 + c)**p/(sqrt( - a**2*x**2 + 1)*a*x** 
3 - sqrt( - a**2*x**2 + 1)*x**2),x)*a*x**2 - 2*int(( - a**2*c*x**2 + c)**p 
/(sqrt( - a**2*x**2 + 1)*a*x**2 - sqrt( - a**2*x**2 + 1)*x),x)*a**3*x**3 + 
 2*int(( - a**2*c*x**2 + c)**p/(sqrt( - a**2*x**2 + 1)*a*x**2 - sqrt( - a* 
*2*x**2 + 1)*x),x)*a**2*x**2 - 8*int((( - a**2*c*x**2 + c)**p*sqrt( - a**2 
*x**2 + 1)*x)/(a**3*x**3 - a**2*x**2 - a*x + 1),x)*a**5*p**3*x**3 + 12*int 
((( - a**2*c*x**2 + c)**p*sqrt( - a**2*x**2 + 1)*x)/(a**3*x**3 - a**2*x**2 
 - a*x + 1),x)*a**5*p**2*x**3 - 4*int((( - a**2*c*x**2 + c)**p*sqrt( - a** 
2*x**2 + 1)*x)/(a**3*x**3 - a**2*x**2 - a*x + 1),x)*a**5*p*x**3 + 8*int((( 
 - a**2*c*x**2 + c)**p*sqrt( - a**2*x**2 + 1)*x)/(a**3*x**3 - a**2*x**2 - 
a*x + 1),x)*a**4*p**3*x**2 - 12*int((( - a**2*c*x**2 + c)**p*sqrt( - a**2* 
x**2 + 1)*x)/(a**3*x**3 - a**2*x**2 - a*x + 1),x)*a**4*p**2*x**2 + 4*int(( 
( - a**2*c*x**2 + c)**p*sqrt( - a**2*x**2 + 1)*x)/(a**3*x**3 - a**2*x**2 - 
 a*x + 1),x)*a**4*p*x**2 - 2*int((( - a**2*c*x**2 + c)**p*sqrt( - a**2*...