3.5.57 \(\int e^{\arcsin (a+b x)^2} x^3 \, dx\) [457]

3.5.57.1 Optimal result

Integrand size = 14, antiderivative size = 381 \[ \int e^{\arcsin (a+b x)^2} x^3 \, dx=\frac {e \sqrt {\pi } \text {erf}(1-i \arcsin (a+b x))}{16 b^4}+\frac {3 a^2 e \sqrt {\pi } \text {erf}(1-i \arcsin (a+b x))}{8 b^4}-\frac {e^4 \sqrt {\pi } \text {erf}(2-i \arcsin (a+b x))}{32 b^4}+\frac {e \sqrt {\pi } \text {erf}(1+i \arcsin (a+b x))}{16 b^4}+\frac {3 a^2 e \sqrt {\pi } \text {erf}(1+i \arcsin (a+b x))}{8 b^4}-\frac {e^4 \sqrt {\pi } \text {erf}(2+i \arcsin (a+b x))}{32 b^4}-\frac {3 a \sqrt [4]{e} \sqrt {\pi } \text {erfi}\left (\frac {1}{2} (-i+2 \arcsin (a+b x))\right )}{16 b^4}-\frac {a^3 \sqrt [4]{e} \sqrt {\pi } \text {erfi}\left (\frac {1}{2} (-i+2 \arcsin (a+b x))\right )}{4 b^4}-\frac {3 a \sqrt [4]{e} \sqrt {\pi } \text {erfi}\left (\frac {1}{2} (i+2 \arcsin (a+b x))\right )}{16 b^4}-\frac {a^3 \sqrt [4]{e} \sqrt {\pi } \text {erfi}\left (\frac {1}{2} (i+2 \arcsin (a+b x))\right )}{4 b^4}+\frac {3 a e^{9/4} \sqrt {\pi } \text {erfi}\left (\frac {1}{2} (-3 i+2 \arcsin (a+b x))\right )}{16 b^4}+\frac {3 a e^{9/4} \sqrt {\pi } \text {erfi}\left (\frac {1}{2} (3 i+2 \arcsin (a+b x))\right )}{16 b^4} \]

output

1/16*I*exp(1)*erfi(-I+arcsin(b*x+a))*Pi^(1/2)/b^4+3/8*I*a^2*exp(1)*erfi(-I 
+arcsin(b*x+a))*Pi^(1/2)/b^4-1/16*I*exp(1)*erfi(I+arcsin(b*x+a))*Pi^(1/2)/ 
b^4-3/8*I*a^2*exp(1)*erfi(I+arcsin(b*x+a))*Pi^(1/2)/b^4-1/32*I*exp(4)*erfi 
(-2*I+arcsin(b*x+a))*Pi^(1/2)/b^4+1/32*I*exp(4)*erfi(2*I+arcsin(b*x+a))*Pi 
^(1/2)/b^4-3/16*a*exp(1/4)*erfi(-1/2*I+arcsin(b*x+a))*Pi^(1/2)/b^4-1/4*a^3 
*exp(1/4)*erfi(-1/2*I+arcsin(b*x+a))*Pi^(1/2)/b^4-3/16*a*exp(1/4)*erfi(1/2 
*I+arcsin(b*x+a))*Pi^(1/2)/b^4-1/4*a^3*exp(1/4)*erfi(1/2*I+arcsin(b*x+a))* 
Pi^(1/2)/b^4+3/16*a*exp(9/4)*erfi(-3/2*I+arcsin(b*x+a))*Pi^(1/2)/b^4+3/16* 
a*exp(9/4)*erfi(3/2*I+arcsin(b*x+a))*Pi^(1/2)/b^4
 

3.5.57.2 Mathematica [A] (verified)

Time = 0.28 (sec) , antiderivative size = 221, normalized size of antiderivative = 0.58 \[ \int e^{\arcsin (a+b x)^2} x^3 \, dx=-\frac {\sqrt {\pi } \left (-2 \left (e+6 a^2 e\right ) \text {erf}(1-i \arcsin (a+b x))+e^4 \text {erf}(2-i \arcsin (a+b x))+\sqrt [4]{e} \left (-2 i a \left (3+4 a^2\right ) \text {erf}\left (\frac {1}{2}+i \arcsin (a+b x)\right )-2 \left (1+6 a^2\right ) e^{3/4} \text {erf}(1+i \arcsin (a+b x))+6 i a e^2 \text {erf}\left (\frac {3}{2}+i \arcsin (a+b x)\right )+e^{15/4} \text {erf}(2+i \arcsin (a+b x))+6 a \text {erfi}\left (\frac {1}{2} (i+2 \arcsin (a+b x))\right )+8 a^3 \text {erfi}\left (\frac {1}{2} (i+2 \arcsin (a+b x))\right )-6 a e^2 \text {erfi}\left (\frac {1}{2} (3 i+2 \arcsin (a+b x))\right )\right )\right )}{32 b^4} \]

input

Integrate[E^ArcSin[a + b*x]^2*x^3,x]
 

output

-1/32*(Sqrt[Pi]*(-2*(E + 6*a^2*E)*Erf[1 - I*ArcSin[a + b*x]] + E^4*Erf[2 - 
 I*ArcSin[a + b*x]] + E^(1/4)*((-2*I)*a*(3 + 4*a^2)*Erf[1/2 + I*ArcSin[a + 
 b*x]] - 2*(1 + 6*a^2)*E^(3/4)*Erf[1 + I*ArcSin[a + b*x]] + (6*I)*a*E^2*Er 
f[3/2 + I*ArcSin[a + b*x]] + E^(15/4)*Erf[2 + I*ArcSin[a + b*x]] + 6*a*Erf 
i[(I + 2*ArcSin[a + b*x])/2] + 8*a^3*Erfi[(I + 2*ArcSin[a + b*x])/2] - 6*a 
*E^2*Erfi[(3*I + 2*ArcSin[a + b*x])/2])))/b^4
 

3.5.57.3 Rubi [A] (verified)

Time = 0.83 (sec) , antiderivative size = 350, normalized size of antiderivative = 0.92, number of steps used = 7, number of rules used = 6, \(\frac {\text {number of rules}}{\text {integrand size}}\) = 0.429, Rules used = {5335, 25, 7292, 27, 7293, 2009}

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

\(\Big \downarrow \) 5335

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

\(\Big \downarrow \) 25

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

\(\Big \downarrow \) 7292

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

\(\Big \downarrow \) 27

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

\(\Big \downarrow \) 7293

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

\(\Big \downarrow \) 2009

\(\displaystyle -\frac {\frac {1}{4} \sqrt [4]{e} \sqrt {\pi } a^3 \text {erfi}\left (\frac {1}{2} (2 \arcsin (a+b x)-i)\right )+\frac {1}{4} \sqrt [4]{e} \sqrt {\pi } a^3 \text {erfi}\left (\frac {1}{2} (2 \arcsin (a+b x)+i)\right )-\frac {3}{8} e \sqrt {\pi } a^2 \text {erf}(1-i \arcsin (a+b x))-\frac {3}{8} e \sqrt {\pi } a^2 \text {erf}(1+i \arcsin (a+b x))-\frac {1}{16} e \sqrt {\pi } \text {erf}(1-i \arcsin (a+b x))+\frac {1}{32} e^4 \sqrt {\pi } \text {erf}(2-i \arcsin (a+b x))-\frac {1}{16} e \sqrt {\pi } \text {erf}(1+i \arcsin (a+b x))+\frac {1}{32} e^4 \sqrt {\pi } \text {erf}(2+i \arcsin (a+b x))+\frac {3}{16} \sqrt [4]{e} \sqrt {\pi } a \text {erfi}\left (\frac {1}{2} (2 \arcsin (a+b x)-i)\right )+\frac {3}{16} \sqrt [4]{e} \sqrt {\pi } a \text {erfi}\left (\frac {1}{2} (2 \arcsin (a+b x)+i)\right )-\frac {3}{16} e^{9/4} \sqrt {\pi } a \text {erfi}\left (\frac {1}{2} (2 \arcsin (a+b x)-3 i)\right )-\frac {3}{16} e^{9/4} \sqrt {\pi } a \text {erfi}\left (\frac {1}{2} (2 \arcsin (a+b x)+3 i)\right )}{b^4}\)

input

Int[E^ArcSin[a + b*x]^2*x^3,x]
 

output

-((-1/16*(E*Sqrt[Pi]*Erf[1 - I*ArcSin[a + b*x]]) - (3*a^2*E*Sqrt[Pi]*Erf[1 
 - I*ArcSin[a + b*x]])/8 + (E^4*Sqrt[Pi]*Erf[2 - I*ArcSin[a + b*x]])/32 - 
(E*Sqrt[Pi]*Erf[1 + I*ArcSin[a + b*x]])/16 - (3*a^2*E*Sqrt[Pi]*Erf[1 + I*A 
rcSin[a + b*x]])/8 + (E^4*Sqrt[Pi]*Erf[2 + I*ArcSin[a + b*x]])/32 + (3*a*E 
^(1/4)*Sqrt[Pi]*Erfi[(-I + 2*ArcSin[a + b*x])/2])/16 + (a^3*E^(1/4)*Sqrt[P 
i]*Erfi[(-I + 2*ArcSin[a + b*x])/2])/4 + (3*a*E^(1/4)*Sqrt[Pi]*Erfi[(I + 2 
*ArcSin[a + b*x])/2])/16 + (a^3*E^(1/4)*Sqrt[Pi]*Erfi[(I + 2*ArcSin[a + b* 
x])/2])/4 - (3*a*E^(9/4)*Sqrt[Pi]*Erfi[(-3*I + 2*ArcSin[a + b*x])/2])/16 - 
 (3*a*E^(9/4)*Sqrt[Pi]*Erfi[(3*I + 2*ArcSin[a + b*x])/2])/16)/b^4)
 

3.5.57.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 2009

Int[u_, x_Symbol] :> Simp[IntSum[u, x], x] /; SumQ[u]
 

rule 5335

Int[(u_.)*(f_)^(ArcSin[(a_.) + (b_.)*(x_)]^(n_.)*(c_.)), x_Symbol] :> Simp[ 
1/b   Subst[Int[(u /. x -> -a/b + Sin[x]/b)*f^(c*x^n)*Cos[x], x], x, ArcSin 
[a + b*x]], x] /; FreeQ[{a, b, c, f}, x] && IGtQ[n, 0]
 

rule 7292

Int[u_, x_Symbol] :> With[{v = NormalizeIntegrand[u, x]}, Int[v, x] /; v =! 
= u]
 

rule 7293

Int[u_, x_Symbol] :> With[{v = ExpandIntegrand[u, x]}, Int[v, x] /; SumQ[v] 
]
 

3.5.57.4 Maple [F]

\[\int {\mathrm e}^{\arcsin \left (b x +a \right )^{2}} x^{3}d x\]

input

int(exp(arcsin(b*x+a)^2)*x^3,x)
 

output

int(exp(arcsin(b*x+a)^2)*x^3,x)
 

3.5.57.5 Fricas [F]

\[ \int e^{\arcsin (a+b x)^2} x^3 \, dx=\int { x^{3} e^{\left (\arcsin \left (b x + a\right )^{2}\right )} \,d x } \]

input

integrate(exp(arcsin(b*x+a)^2)*x^3,x, algorithm="fricas")
 

output

integral(x^3*e^(arcsin(b*x + a)^2), x)
 

3.5.57.6 Sympy [F]

\[ \int e^{\arcsin (a+b x)^2} x^3 \, dx=\int x^{3} e^{\operatorname {asin}^{2}{\left (a + b x \right )}}\, dx \]

input

integrate(exp(asin(b*x+a)**2)*x**3,x)
 

output

Integral(x**3*exp(asin(a + b*x)**2), x)
 

3.5.57.7 Maxima [F]

\[ \int e^{\arcsin (a+b x)^2} x^3 \, dx=\int { x^{3} e^{\left (\arcsin \left (b x + a\right )^{2}\right )} \,d x } \]

input

integrate(exp(arcsin(b*x+a)^2)*x^3,x, algorithm="maxima")
 

output

integrate(x^3*e^(arcsin(b*x + a)^2), x)
 

3.5.57.8 Giac [F]

\[ \int e^{\arcsin (a+b x)^2} x^3 \, dx=\int { x^{3} e^{\left (\arcsin \left (b x + a\right )^{2}\right )} \,d x } \]

input

integrate(exp(arcsin(b*x+a)^2)*x^3,x, algorithm="giac")
 

output

integrate(x^3*e^(arcsin(b*x + a)^2), x)
 

3.5.57.9 Mupad [F(-1)]

Timed out. \[ \int e^{\arcsin (a+b x)^2} x^3 \, dx=\int x^3\,{\mathrm {e}}^{{\mathrm {asin}\left (a+b\,x\right )}^2} \,d x \]

input

int(x^3*exp(asin(a + b*x)^2),x)
 

output

int(x^3*exp(asin(a + b*x)^2), x)