\(\int \frac {\arcsin (a x)^4}{x^3} \, dx\) [37]

Optimal result

Integrand size = 10, antiderivative size = 119 \[ \int \frac {\arcsin (a x)^4}{x^3} \, dx=-2 i a^2 \arcsin (a x)^3-\frac {2 a \sqrt {1-a^2 x^2} \arcsin (a x)^3}{x}-\frac {\arcsin (a x)^4}{2 x^2}+6 a^2 \arcsin (a x)^2 \log \left (1-e^{2 i \arcsin (a x)}\right )-6 i a^2 \arcsin (a x) \operatorname {PolyLog}\left (2,e^{2 i \arcsin (a x)}\right )+3 a^2 \operatorname {PolyLog}\left (3,e^{2 i \arcsin (a x)}\right ) \] Output:

-2*I*a^2*arcsin(a*x)^3-2*a*(-a^2*x^2+1)^(1/2)*arcsin(a*x)^3/x-1/2*arcsin(a 
*x)^4/x^2+6*a^2*arcsin(a*x)^2*ln(1-(I*a*x+(-a^2*x^2+1)^(1/2))^2)-6*I*a^2*a 
rcsin(a*x)*polylog(2,(I*a*x+(-a^2*x^2+1)^(1/2))^2)+3*a^2*polylog(3,(I*a*x+ 
(-a^2*x^2+1)^(1/2))^2)
 

Mathematica [A] (verified)

Time = 0.21 (sec) , antiderivative size = 124, normalized size of antiderivative = 1.04 \[ \int \frac {\arcsin (a x)^4}{x^3} \, dx=-\frac {\arcsin (a x)^4}{2 x^2}+\frac {1}{4} a^2 \left (-i \pi ^3+8 i \arcsin (a x)^3-\frac {8 \sqrt {1-a^2 x^2} \arcsin (a x)^3}{a x}+24 \arcsin (a x)^2 \log \left (1-e^{-2 i \arcsin (a x)}\right )+24 i \arcsin (a x) \operatorname {PolyLog}\left (2,e^{-2 i \arcsin (a x)}\right )+12 \operatorname {PolyLog}\left (3,e^{-2 i \arcsin (a x)}\right )\right ) \] Input:

Integrate[ArcSin[a*x]^4/x^3,x]
 

Output:

-1/2*ArcSin[a*x]^4/x^2 + (a^2*((-I)*Pi^3 + (8*I)*ArcSin[a*x]^3 - (8*Sqrt[1 
 - a^2*x^2]*ArcSin[a*x]^3)/(a*x) + 24*ArcSin[a*x]^2*Log[1 - E^((-2*I)*ArcS 
in[a*x])] + (24*I)*ArcSin[a*x]*PolyLog[2, E^((-2*I)*ArcSin[a*x])] + 12*Pol 
yLog[3, E^((-2*I)*ArcSin[a*x])]))/4
 

Rubi [A] (verified)

Time = 0.88 (sec) , antiderivative size = 134, normalized size of antiderivative = 1.13, number of steps used = 12, number of rules used = 11, \(\frac {\text {number of rules}}{\text {integrand size}}\) = 1.100, Rules used = {5138, 5186, 5136, 3042, 25, 4200, 25, 2620, 3011, 2720, 7143}

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 {\arcsin (a x)^4}{x^3} \, dx\)

\(\Big \downarrow \) 5138

\(\displaystyle 2 a \int \frac {\arcsin (a x)^3}{x^2 \sqrt {1-a^2 x^2}}dx-\frac {\arcsin (a x)^4}{2 x^2}\)

\(\Big \downarrow \) 5186

\(\displaystyle 2 a \left (3 a \int \frac {\arcsin (a x)^2}{x}dx-\frac {\sqrt {1-a^2 x^2} \arcsin (a x)^3}{x}\right )-\frac {\arcsin (a x)^4}{2 x^2}\)

\(\Big \downarrow \) 5136

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

\(\Big \downarrow \) 3042

\(\displaystyle 2 a \left (3 a \int -\arcsin (a x)^2 \tan \left (\arcsin (a x)+\frac {\pi }{2}\right )d\arcsin (a x)-\frac {\sqrt {1-a^2 x^2} \arcsin (a x)^3}{x}\right )-\frac {\arcsin (a x)^4}{2 x^2}\)

\(\Big \downarrow \) 25

\(\displaystyle 2 a \left (-3 a \int \arcsin (a x)^2 \tan \left (\arcsin (a x)+\frac {\pi }{2}\right )d\arcsin (a x)-\frac {\sqrt {1-a^2 x^2} \arcsin (a x)^3}{x}\right )-\frac {\arcsin (a x)^4}{2 x^2}\)

\(\Big \downarrow \) 4200

\(\displaystyle -\frac {\arcsin (a x)^4}{2 x^2}+2 a \left (-\frac {\sqrt {1-a^2 x^2} \arcsin (a x)^3}{x}+3 a \left (2 i \int -\frac {e^{2 i \arcsin (a x)} \arcsin (a x)^2}{1-e^{2 i \arcsin (a x)}}d\arcsin (a x)-\frac {1}{3} i \arcsin (a x)^3\right )\right )\)

\(\Big \downarrow \) 25

\(\displaystyle -\frac {\arcsin (a x)^4}{2 x^2}+2 a \left (-\frac {\sqrt {1-a^2 x^2} \arcsin (a x)^3}{x}+3 a \left (-2 i \int \frac {e^{2 i \arcsin (a x)} \arcsin (a x)^2}{1-e^{2 i \arcsin (a x)}}d\arcsin (a x)-\frac {1}{3} i \arcsin (a x)^3\right )\right )\)

\(\Big \downarrow \) 2620

\(\displaystyle -\frac {\arcsin (a x)^4}{2 x^2}+2 a \left (-\frac {\sqrt {1-a^2 x^2} \arcsin (a x)^3}{x}+3 a \left (-2 i \left (\frac {1}{2} i \arcsin (a x)^2 \log \left (1-e^{2 i \arcsin (a x)}\right )-i \int \arcsin (a x) \log \left (1-e^{2 i \arcsin (a x)}\right )d\arcsin (a x)\right )-\frac {1}{3} i \arcsin (a x)^3\right )\right )\)

\(\Big \downarrow \) 3011

\(\displaystyle -\frac {\arcsin (a x)^4}{2 x^2}+2 a \left (-\frac {\sqrt {1-a^2 x^2} \arcsin (a x)^3}{x}+3 a \left (-2 i \left (\frac {1}{2} i \arcsin (a x)^2 \log \left (1-e^{2 i \arcsin (a x)}\right )-i \left (\frac {1}{2} i \arcsin (a x) \operatorname {PolyLog}\left (2,e^{2 i \arcsin (a x)}\right )-\frac {1}{2} i \int \operatorname {PolyLog}\left (2,e^{2 i \arcsin (a x)}\right )d\arcsin (a x)\right )\right )-\frac {1}{3} i \arcsin (a x)^3\right )\right )\)

\(\Big \downarrow \) 2720

\(\displaystyle -\frac {\arcsin (a x)^4}{2 x^2}+2 a \left (-\frac {\sqrt {1-a^2 x^2} \arcsin (a x)^3}{x}+3 a \left (-2 i \left (\frac {1}{2} i \arcsin (a x)^2 \log \left (1-e^{2 i \arcsin (a x)}\right )-i \left (\frac {1}{2} i \arcsin (a x) \operatorname {PolyLog}\left (2,e^{2 i \arcsin (a x)}\right )-\frac {1}{4} \int e^{-2 i \arcsin (a x)} \operatorname {PolyLog}\left (2,e^{2 i \arcsin (a x)}\right )de^{2 i \arcsin (a x)}\right )\right )-\frac {1}{3} i \arcsin (a x)^3\right )\right )\)

\(\Big \downarrow \) 7143

\(\displaystyle -\frac {\arcsin (a x)^4}{2 x^2}+2 a \left (-\frac {\sqrt {1-a^2 x^2} \arcsin (a x)^3}{x}+3 a \left (-2 i \left (\frac {1}{2} i \arcsin (a x)^2 \log \left (1-e^{2 i \arcsin (a x)}\right )-i \left (\frac {1}{2} i \arcsin (a x) \operatorname {PolyLog}\left (2,e^{2 i \arcsin (a x)}\right )-\frac {1}{4} \operatorname {PolyLog}\left (3,e^{2 i \arcsin (a x)}\right )\right )\right )-\frac {1}{3} i \arcsin (a x)^3\right )\right )\)

Input:

Int[ArcSin[a*x]^4/x^3,x]
 

Output:

-1/2*ArcSin[a*x]^4/x^2 + 2*a*(-((Sqrt[1 - a^2*x^2]*ArcSin[a*x]^3)/x) + 3*a 
*((-1/3*I)*ArcSin[a*x]^3 - (2*I)*((I/2)*ArcSin[a*x]^2*Log[1 - E^((2*I)*Arc 
Sin[a*x])] - I*((I/2)*ArcSin[a*x]*PolyLog[2, E^((2*I)*ArcSin[a*x])] - Poly 
Log[3, E^((2*I)*ArcSin[a*x])]/4))))
 

Defintions of rubi rules used

rule 25

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

rule 2620

Int[(((F_)^((g_.)*((e_.) + (f_.)*(x_))))^(n_.)*((c_.) + (d_.)*(x_))^(m_.))/ 
((a_) + (b_.)*((F_)^((g_.)*((e_.) + (f_.)*(x_))))^(n_.)), x_Symbol] :> Simp 
[((c + d*x)^m/(b*f*g*n*Log[F]))*Log[1 + b*((F^(g*(e + f*x)))^n/a)], x] - Si 
mp[d*(m/(b*f*g*n*Log[F]))   Int[(c + d*x)^(m - 1)*Log[1 + b*((F^(g*(e + f*x 
)))^n/a)], x], x] /; FreeQ[{F, a, b, c, d, e, f, g, n}, x] && IGtQ[m, 0]
 

rule 2720

Int[u_, x_Symbol] :> With[{v = FunctionOfExponential[u, x]}, Simp[v/D[v, x] 
   Subst[Int[FunctionOfExponentialFunction[u, x]/x, x], x, v], x]] /; Funct 
ionOfExponentialQ[u, x] &&  !MatchQ[u, (w_)*((a_.)*(v_)^(n_))^(m_) /; FreeQ 
[{a, m, n}, x] && IntegerQ[m*n]] &&  !MatchQ[u, E^((c_.)*((a_.) + (b_.)*x)) 
*(F_)[v_] /; FreeQ[{a, b, c}, x] && InverseFunctionQ[F[x]]]
 

rule 3011

Int[Log[1 + (e_.)*((F_)^((c_.)*((a_.) + (b_.)*(x_))))^(n_.)]*((f_.) + (g_.) 
*(x_))^(m_.), x_Symbol] :> Simp[(-(f + g*x)^m)*(PolyLog[2, (-e)*(F^(c*(a + 
b*x)))^n]/(b*c*n*Log[F])), x] + Simp[g*(m/(b*c*n*Log[F]))   Int[(f + g*x)^( 
m - 1)*PolyLog[2, (-e)*(F^(c*(a + b*x)))^n], x], x] /; FreeQ[{F, a, b, c, e 
, f, g, n}, x] && GtQ[m, 0]
 

rule 3042

Int[u_, x_Symbol] :> Int[DeactivateTrig[u, x], x] /; FunctionOfTrigOfLinear 
Q[u, x]
 

rule 4200

Int[((c_.) + (d_.)*(x_))^(m_.)*tan[(e_.) + Pi*(k_.) + (f_.)*(x_)], x_Symbol 
] :> Simp[I*((c + d*x)^(m + 1)/(d*(m + 1))), x] - Simp[2*I   Int[(c + d*x)^ 
m*E^(2*I*k*Pi)*(E^(2*I*(e + f*x))/(1 + E^(2*I*k*Pi)*E^(2*I*(e + f*x)))), x] 
, x] /; FreeQ[{c, d, e, f}, x] && IntegerQ[4*k] && IGtQ[m, 0]
 

rule 5136

Int[((a_.) + ArcSin[(c_.)*(x_)]*(b_.))^(n_.)/(x_), x_Symbol] :> Subst[Int[( 
a + b*x)^n*Cot[x], x], x, ArcSin[c*x]] /; FreeQ[{a, b, c}, x] && IGtQ[n, 0]
 

rule 5138

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

rule 5186

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

rule 7143

Int[PolyLog[n_, (c_.)*((a_.) + (b_.)*(x_))^(p_.)]/((d_.) + (e_.)*(x_)), x_S 
ymbol] :> Simp[PolyLog[n + 1, c*(a + b*x)^p]/(e*p), x] /; FreeQ[{a, b, c, d 
, e, n, p}, x] && EqQ[b*d, a*e]
 

Maple [A] (verified)

Time = 0.29 (sec) , antiderivative size = 219, normalized size of antiderivative = 1.84

method	result	size
derivativedivides	\(a^{2} \left (-\frac {\arcsin \left (a x \right )^{3} \left (-4 i a^{2} x^{2}+4 \sqrt {-a^{2} x^{2}+1}\, x a +\arcsin \left (a x \right )\right )}{2 a^{2} x^{2}}-4 i \arcsin \left (a x \right )^{3}+6 \arcsin \left (a x \right )^{2} \ln \left (1+i a x +\sqrt {-a^{2} x^{2}+1}\right )-12 i \arcsin \left (a x \right ) \operatorname {polylog}\left (2, -i a x -\sqrt {-a^{2} x^{2}+1}\right )+12 \operatorname {polylog}\left (3, -i a x -\sqrt {-a^{2} x^{2}+1}\right )+6 \arcsin \left (a x \right )^{2} \ln \left (1-i a x -\sqrt {-a^{2} x^{2}+1}\right )-12 i \arcsin \left (a x \right ) \operatorname {polylog}\left (2, i a x +\sqrt {-a^{2} x^{2}+1}\right )+12 \operatorname {polylog}\left (3, i a x +\sqrt {-a^{2} x^{2}+1}\right )\right )\)	\(219\)
default	\(a^{2} \left (-\frac {\arcsin \left (a x \right )^{3} \left (-4 i a^{2} x^{2}+4 \sqrt {-a^{2} x^{2}+1}\, x a +\arcsin \left (a x \right )\right )}{2 a^{2} x^{2}}-4 i \arcsin \left (a x \right )^{3}+6 \arcsin \left (a x \right )^{2} \ln \left (1+i a x +\sqrt {-a^{2} x^{2}+1}\right )-12 i \arcsin \left (a x \right ) \operatorname {polylog}\left (2, -i a x -\sqrt {-a^{2} x^{2}+1}\right )+12 \operatorname {polylog}\left (3, -i a x -\sqrt {-a^{2} x^{2}+1}\right )+6 \arcsin \left (a x \right )^{2} \ln \left (1-i a x -\sqrt {-a^{2} x^{2}+1}\right )-12 i \arcsin \left (a x \right ) \operatorname {polylog}\left (2, i a x +\sqrt {-a^{2} x^{2}+1}\right )+12 \operatorname {polylog}\left (3, i a x +\sqrt {-a^{2} x^{2}+1}\right )\right )\)	\(219\)

Input:

int(arcsin(a*x)^4/x^3,x,method=_RETURNVERBOSE)
 

Output:

a^2*(-1/2*arcsin(a*x)^3*(-4*I*a^2*x^2+4*(-a^2*x^2+1)^(1/2)*x*a+arcsin(a*x) 
)/a^2/x^2-4*I*arcsin(a*x)^3+6*arcsin(a*x)^2*ln(1+I*a*x+(-a^2*x^2+1)^(1/2)) 
-12*I*arcsin(a*x)*polylog(2,-I*a*x-(-a^2*x^2+1)^(1/2))+12*polylog(3,-I*a*x 
-(-a^2*x^2+1)^(1/2))+6*arcsin(a*x)^2*ln(1-I*a*x-(-a^2*x^2+1)^(1/2))-12*I*a 
rcsin(a*x)*polylog(2,I*a*x+(-a^2*x^2+1)^(1/2))+12*polylog(3,I*a*x+(-a^2*x^ 
2+1)^(1/2)))
 

Fricas [F]

\[ \int \frac {\arcsin (a x)^4}{x^3} \, dx=\int { \frac {\arcsin \left (a x\right )^{4}}{x^{3}} \,d x } \] Input:

integrate(arcsin(a*x)^4/x^3,x, algorithm="fricas")
 

Output:

integral(arcsin(a*x)^4/x^3, x)
 

Sympy [F]

\[ \int \frac {\arcsin (a x)^4}{x^3} \, dx=\int \frac {\operatorname {asin}^{4}{\left (a x \right )}}{x^{3}}\, dx \] Input:

integrate(asin(a*x)**4/x**3,x)
 

Output:

Integral(asin(a*x)**4/x**3, x)
                                                                                    
                                                                                    
 

Maxima [F]

\[ \int \frac {\arcsin (a x)^4}{x^3} \, dx=\int { \frac {\arcsin \left (a x\right )^{4}}{x^{3}} \,d x } \] Input:

integrate(arcsin(a*x)^4/x^3,x, algorithm="maxima")
 

Output:

-1/2*(arctan2(a*x, sqrt(a*x + 1)*sqrt(-a*x + 1))^4 + 4*a*x^2*integrate(sqr 
t(a*x + 1)*sqrt(-a*x + 1)*arctan2(a*x, sqrt(a*x + 1)*sqrt(-a*x + 1))^3/(a^ 
2*x^4 - x^2), x))/x^2
 

Giac [F]

\[ \int \frac {\arcsin (a x)^4}{x^3} \, dx=\int { \frac {\arcsin \left (a x\right )^{4}}{x^{3}} \,d x } \] Input:

integrate(arcsin(a*x)^4/x^3,x, algorithm="giac")
 

Output:

integrate(arcsin(a*x)^4/x^3, x)
 

Mupad [F(-1)]

Timed out. \[ \int \frac {\arcsin (a x)^4}{x^3} \, dx=\int \frac {{\mathrm {asin}\left (a\,x\right )}^4}{x^3} \,d x \] Input:

int(asin(a*x)^4/x^3,x)
 

Output:

int(asin(a*x)^4/x^3, x)
 

Reduce [F]

\[ \int \frac {\arcsin (a x)^4}{x^3} \, dx=\int \frac {\mathit {asin} \left (a x \right )^{4}}{x^{3}}d x \] Input:

int(asin(a*x)^4/x^3,x)
 

Output:

int(asin(a*x)**4/x**3,x)