3.253 \(\int (a+b \sin ^{-1}(c+d x))^{5/2} \, dx\)

Optimal. Leaf size=204 \[ -\frac{15 \sqrt{\frac{\pi }{2}} b^{5/2} \sin \left (\frac{a}{b}\right ) \text{FresnelC}\left (\frac{\sqrt{\frac{2}{\pi }} \sqrt{a+b \sin ^{-1}(c+d x)}}{\sqrt{b}}\right )}{4 d}+\frac{15 \sqrt{\frac{\pi }{2}} b^{5/2} \cos \left (\frac{a}{b}\right ) S\left (\frac{\sqrt{\frac{2}{\pi }} \sqrt{a+b \sin ^{-1}(c+d x)}}{\sqrt{b}}\right )}{4 d}-\frac{15 b^2 (c+d x) \sqrt{a+b \sin ^{-1}(c+d x)}}{4 d}+\frac{5 b \sqrt{1-(c+d x)^2} \left (a+b \sin ^{-1}(c+d x)\right )^{3/2}}{2 d}+\frac{(c+d x) \left (a+b \sin ^{-1}(c+d x)\right )^{5/2}}{d} \]

[Out]

(-15*b^2*(c + d*x)*Sqrt[a + b*ArcSin[c + d*x]])/(4*d) + (5*b*Sqrt[1 - (c + d*x)^2]*(a + b*ArcSin[c + d*x])^(3/
2))/(2*d) + ((c + d*x)*(a + b*ArcSin[c + d*x])^(5/2))/d + (15*b^(5/2)*Sqrt[Pi/2]*Cos[a/b]*FresnelS[(Sqrt[2/Pi]
*Sqrt[a + b*ArcSin[c + d*x]])/Sqrt[b]])/(4*d) - (15*b^(5/2)*Sqrt[Pi/2]*FresnelC[(Sqrt[2/Pi]*Sqrt[a + b*ArcSin[
c + d*x]])/Sqrt[b]]*Sin[a/b])/(4*d)

________________________________________________________________________________________

Rubi [A]  time = 0.415817, antiderivative size = 204, normalized size of antiderivative = 1., number of steps used = 10, number of rules used = 9, integrand size = 14, \(\frac{\text{number of rules}}{\text{integrand size}}\) = 0.643, Rules used = {4803, 4619, 4677, 4723, 3306, 3305, 3351, 3304, 3352} \[ -\frac{15 \sqrt{\frac{\pi }{2}} b^{5/2} \sin \left (\frac{a}{b}\right ) \text{FresnelC}\left (\frac{\sqrt{\frac{2}{\pi }} \sqrt{a+b \sin ^{-1}(c+d x)}}{\sqrt{b}}\right )}{4 d}+\frac{15 \sqrt{\frac{\pi }{2}} b^{5/2} \cos \left (\frac{a}{b}\right ) S\left (\frac{\sqrt{\frac{2}{\pi }} \sqrt{a+b \sin ^{-1}(c+d x)}}{\sqrt{b}}\right )}{4 d}-\frac{15 b^2 (c+d x) \sqrt{a+b \sin ^{-1}(c+d x)}}{4 d}+\frac{5 b \sqrt{1-(c+d x)^2} \left (a+b \sin ^{-1}(c+d x)\right )^{3/2}}{2 d}+\frac{(c+d x) \left (a+b \sin ^{-1}(c+d x)\right )^{5/2}}{d} \]

Antiderivative was successfully verified.

[In]

Int[(a + b*ArcSin[c + d*x])^(5/2),x]

[Out]

(-15*b^2*(c + d*x)*Sqrt[a + b*ArcSin[c + d*x]])/(4*d) + (5*b*Sqrt[1 - (c + d*x)^2]*(a + b*ArcSin[c + d*x])^(3/
2))/(2*d) + ((c + d*x)*(a + b*ArcSin[c + d*x])^(5/2))/d + (15*b^(5/2)*Sqrt[Pi/2]*Cos[a/b]*FresnelS[(Sqrt[2/Pi]
*Sqrt[a + b*ArcSin[c + d*x]])/Sqrt[b]])/(4*d) - (15*b^(5/2)*Sqrt[Pi/2]*FresnelC[(Sqrt[2/Pi]*Sqrt[a + b*ArcSin[
c + d*x]])/Sqrt[b]]*Sin[a/b])/(4*d)

Rule 4803

Int[((a_.) + ArcSin[(c_) + (d_.)*(x_)]*(b_.))^(n_.), x_Symbol] :> Dist[1/d, Subst[Int[(a + b*ArcSin[x])^n, x],
x, c + d*x], x] /; FreeQ[{a, b, c, d, n}, x]

Rule 4619

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

Rule 4677

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

Rule 4723

Int[((a_.) + ArcSin[(c_.)*(x_)]*(b_.))^(n_.)*(x_)^(m_.)*((d_) + (e_.)*(x_)^2)^(p_.), x_Symbol] :> Dist[d^p/c^(
m + 1), Subst[Int[(a + b*x)^n*Sin[x]^m*Cos[x]^(2*p + 1), x], x, ArcSin[c*x]], x] /; FreeQ[{a, b, c, d, e, n},
x] && EqQ[c^2*d + e, 0] && IntegerQ[2*p] && GtQ[p, -1] && IGtQ[m, 0] && (IntegerQ[p] || GtQ[d, 0])

Rule 3306

Int[sin[(e_.) + (f_.)*(x_)]/Sqrt[(c_.) + (d_.)*(x_)], x_Symbol] :> Dist[Cos[(d*e - c*f)/d], Int[Sin[(c*f)/d +
f*x]/Sqrt[c + d*x], x], x] + Dist[Sin[(d*e - c*f)/d], Int[Cos[(c*f)/d + f*x]/Sqrt[c + d*x], x], x] /; FreeQ[{c
, d, e, f}, x] && ComplexFreeQ[f] && NeQ[d*e - c*f, 0]

Rule 3305

Int[sin[(e_.) + (f_.)*(x_)]/Sqrt[(c_.) + (d_.)*(x_)], x_Symbol] :> Dist[2/d, Subst[Int[Sin[(f*x^2)/d], x], x,
Sqrt[c + d*x]], x] /; FreeQ[{c, d, e, f}, x] && ComplexFreeQ[f] && EqQ[d*e - c*f, 0]

Rule 3351

Int[Sin[(d_.)*((e_.) + (f_.)*(x_))^2], x_Symbol] :> Simp[(Sqrt[Pi/2]*FresnelS[Sqrt[2/Pi]*Rt[d, 2]*(e + f*x)])/
(f*Rt[d, 2]), x] /; FreeQ[{d, e, f}, x]

Rule 3304

Int[sin[Pi/2 + (e_.) + (f_.)*(x_)]/Sqrt[(c_.) + (d_.)*(x_)], x_Symbol] :> Dist[2/d, Subst[Int[Cos[(f*x^2)/d],
x], x, Sqrt[c + d*x]], x] /; FreeQ[{c, d, e, f}, x] && ComplexFreeQ[f] && EqQ[d*e - c*f, 0]

Rule 3352

Int[Cos[(d_.)*((e_.) + (f_.)*(x_))^2], x_Symbol] :> Simp[(Sqrt[Pi/2]*FresnelC[Sqrt[2/Pi]*Rt[d, 2]*(e + f*x)])/
(f*Rt[d, 2]), x] /; FreeQ[{d, e, f}, x]

Rubi steps

\begin{align*} \int \left (a+b \sin ^{-1}(c+d x)\right )^{5/2} \, dx &=\frac{\operatorname{Subst}\left (\int \left (a+b \sin ^{-1}(x)\right )^{5/2} \, dx,x,c+d x\right )}{d}\\ &=\frac{(c+d x) \left (a+b \sin ^{-1}(c+d x)\right )^{5/2}}{d}-\frac{(5 b) \operatorname{Subst}\left (\int \frac{x \left (a+b \sin ^{-1}(x)\right )^{3/2}}{\sqrt{1-x^2}} \, dx,x,c+d x\right )}{2 d}\\ &=\frac{5 b \sqrt{1-(c+d x)^2} \left (a+b \sin ^{-1}(c+d x)\right )^{3/2}}{2 d}+\frac{(c+d x) \left (a+b \sin ^{-1}(c+d x)\right )^{5/2}}{d}-\frac{\left (15 b^2\right ) \operatorname{Subst}\left (\int \sqrt{a+b \sin ^{-1}(x)} \, dx,x,c+d x\right )}{4 d}\\ &=-\frac{15 b^2 (c+d x) \sqrt{a+b \sin ^{-1}(c+d x)}}{4 d}+\frac{5 b \sqrt{1-(c+d x)^2} \left (a+b \sin ^{-1}(c+d x)\right )^{3/2}}{2 d}+\frac{(c+d x) \left (a+b \sin ^{-1}(c+d x)\right )^{5/2}}{d}+\frac{\left (15 b^3\right ) \operatorname{Subst}\left (\int \frac{x}{\sqrt{1-x^2} \sqrt{a+b \sin ^{-1}(x)}} \, dx,x,c+d x\right )}{8 d}\\ &=-\frac{15 b^2 (c+d x) \sqrt{a+b \sin ^{-1}(c+d x)}}{4 d}+\frac{5 b \sqrt{1-(c+d x)^2} \left (a+b \sin ^{-1}(c+d x)\right )^{3/2}}{2 d}+\frac{(c+d x) \left (a+b \sin ^{-1}(c+d x)\right )^{5/2}}{d}+\frac{\left (15 b^3\right ) \operatorname{Subst}\left (\int \frac{\sin (x)}{\sqrt{a+b x}} \, dx,x,\sin ^{-1}(c+d x)\right )}{8 d}\\ &=-\frac{15 b^2 (c+d x) \sqrt{a+b \sin ^{-1}(c+d x)}}{4 d}+\frac{5 b \sqrt{1-(c+d x)^2} \left (a+b \sin ^{-1}(c+d x)\right )^{3/2}}{2 d}+\frac{(c+d x) \left (a+b \sin ^{-1}(c+d x)\right )^{5/2}}{d}+\frac{\left (15 b^3 \cos \left (\frac{a}{b}\right )\right ) \operatorname{Subst}\left (\int \frac{\sin \left (\frac{a}{b}+x\right )}{\sqrt{a+b x}} \, dx,x,\sin ^{-1}(c+d x)\right )}{8 d}-\frac{\left (15 b^3 \sin \left (\frac{a}{b}\right )\right ) \operatorname{Subst}\left (\int \frac{\cos \left (\frac{a}{b}+x\right )}{\sqrt{a+b x}} \, dx,x,\sin ^{-1}(c+d x)\right )}{8 d}\\ &=-\frac{15 b^2 (c+d x) \sqrt{a+b \sin ^{-1}(c+d x)}}{4 d}+\frac{5 b \sqrt{1-(c+d x)^2} \left (a+b \sin ^{-1}(c+d x)\right )^{3/2}}{2 d}+\frac{(c+d x) \left (a+b \sin ^{-1}(c+d x)\right )^{5/2}}{d}+\frac{\left (15 b^2 \cos \left (\frac{a}{b}\right )\right ) \operatorname{Subst}\left (\int \sin \left (\frac{x^2}{b}\right ) \, dx,x,\sqrt{a+b \sin ^{-1}(c+d x)}\right )}{4 d}-\frac{\left (15 b^2 \sin \left (\frac{a}{b}\right )\right ) \operatorname{Subst}\left (\int \cos \left (\frac{x^2}{b}\right ) \, dx,x,\sqrt{a+b \sin ^{-1}(c+d x)}\right )}{4 d}\\ &=-\frac{15 b^2 (c+d x) \sqrt{a+b \sin ^{-1}(c+d x)}}{4 d}+\frac{5 b \sqrt{1-(c+d x)^2} \left (a+b \sin ^{-1}(c+d x)\right )^{3/2}}{2 d}+\frac{(c+d x) \left (a+b \sin ^{-1}(c+d x)\right )^{5/2}}{d}+\frac{15 b^{5/2} \sqrt{\frac{\pi }{2}} \cos \left (\frac{a}{b}\right ) S\left (\frac{\sqrt{\frac{2}{\pi }} \sqrt{a+b \sin ^{-1}(c+d x)}}{\sqrt{b}}\right )}{4 d}-\frac{15 b^{5/2} \sqrt{\frac{\pi }{2}} C\left (\frac{\sqrt{\frac{2}{\pi }} \sqrt{a+b \sin ^{-1}(c+d x)}}{\sqrt{b}}\right ) \sin \left (\frac{a}{b}\right )}{4 d}\\ \end{align*}

Mathematica [C]  time = 1.57182, size = 432, normalized size = 2.12 \[ \frac{e^{-\frac{i a}{b}} \left (2 b \left (2 a^2 \sqrt{-\frac{i \left (a+b \sin ^{-1}(c+d x)\right )}{b}} \text{Gamma}\left (\frac{3}{2},-\frac{i \left (a+b \sin ^{-1}(c+d x)\right )}{b}\right )+2 a^2 e^{\frac{2 i a}{b}} \sqrt{\frac{i \left (a+b \sin ^{-1}(c+d x)\right )}{b}} \text{Gamma}\left (\frac{3}{2},\frac{i \left (a+b \sin ^{-1}(c+d x)\right )}{b}\right )+e^{\frac{i a}{b}} \left (a+b \sin ^{-1}(c+d x)\right ) \left (2 \sin ^{-1}(c+d x) \left (4 a (c+d x)+5 b \sqrt{-c^2-2 c d x-d^2 x^2+1}\right )+10 a \sqrt{-c^2-2 c d x-d^2 x^2+1}-15 b (c+d x)+4 b (c+d x) \sin ^{-1}(c+d x)^2\right )\right )+\frac{i \sqrt{\frac{\pi }{2}} \left (4 a^2+15 b^2\right ) \left (-1+e^{\frac{2 i a}{b}}\right ) \sqrt{a+b \sin ^{-1}(c+d x)} \text{FresnelC}\left (\sqrt{\frac{2}{\pi }} \sqrt{\frac{1}{b}} \sqrt{a+b \sin ^{-1}(c+d x)}\right )}{\sqrt{\frac{1}{b}}}+\frac{\sqrt{\frac{\pi }{2}} \left (4 a^2+15 b^2\right ) \left (1+e^{\frac{2 i a}{b}}\right ) \sqrt{a+b \sin ^{-1}(c+d x)} S\left (\sqrt{\frac{1}{b}} \sqrt{\frac{2}{\pi }} \sqrt{a+b \sin ^{-1}(c+d x)}\right )}{\sqrt{\frac{1}{b}}}\right )}{8 d \sqrt{a+b \sin ^{-1}(c+d x)}} \]

Warning: Unable to verify antiderivative.

[In]

Integrate[(a + b*ArcSin[c + d*x])^(5/2),x]

[Out]

((I*(4*a^2 + 15*b^2)*(-1 + E^(((2*I)*a)/b))*Sqrt[Pi/2]*Sqrt[a + b*ArcSin[c + d*x]]*FresnelC[Sqrt[b^(-1)]*Sqrt[
2/Pi]*Sqrt[a + b*ArcSin[c + d*x]]])/Sqrt[b^(-1)] + ((4*a^2 + 15*b^2)*(1 + E^(((2*I)*a)/b))*Sqrt[Pi/2]*Sqrt[a +
b*ArcSin[c + d*x]]*FresnelS[Sqrt[b^(-1)]*Sqrt[2/Pi]*Sqrt[a + b*ArcSin[c + d*x]]])/Sqrt[b^(-1)] + 2*b*(E^((I*a
)/b)*(a + b*ArcSin[c + d*x])*(-15*b*(c + d*x) + 10*a*Sqrt[1 - c^2 - 2*c*d*x - d^2*x^2] + 2*(4*a*(c + d*x) + 5*
b*Sqrt[1 - c^2 - 2*c*d*x - d^2*x^2])*ArcSin[c + d*x] + 4*b*(c + d*x)*ArcSin[c + d*x]^2) + 2*a^2*Sqrt[((-I)*(a
+ b*ArcSin[c + d*x]))/b]*Gamma[3/2, ((-I)*(a + b*ArcSin[c + d*x]))/b] + 2*a^2*E^(((2*I)*a)/b)*Sqrt[(I*(a + b*A
rcSin[c + d*x]))/b]*Gamma[3/2, (I*(a + b*ArcSin[c + d*x]))/b]))/(8*d*E^((I*a)/b)*Sqrt[a + b*ArcSin[c + d*x]])

________________________________________________________________________________________

Maple [B]  time = 0., size = 433, normalized size = 2.1 \begin{align*}{\frac{1}{8\,d} \left ( 15\,\sqrt{{b}^{-1}}\sqrt{\pi }\sqrt{2}\sqrt{a+b\arcsin \left ( dx+c \right ) }\cos \left ({\frac{a}{b}} \right ){\it FresnelS} \left ({\frac{\sqrt{2}\sqrt{a+b\arcsin \left ( dx+c \right ) }}{\sqrt{{b}^{-1}}\sqrt{\pi }b}} \right ){b}^{3}-15\,\sqrt{{b}^{-1}}\sqrt{\pi }\sqrt{2}\sqrt{a+b\arcsin \left ( dx+c \right ) }\sin \left ({\frac{a}{b}} \right ){\it FresnelC} \left ({\frac{\sqrt{2}\sqrt{a+b\arcsin \left ( dx+c \right ) }}{\sqrt{{b}^{-1}}\sqrt{\pi }b}} \right ){b}^{3}+8\, \left ( \arcsin \left ( dx+c \right ) \right ) ^{3}\sin \left ({\frac{a+b\arcsin \left ( dx+c \right ) }{b}}-{\frac{a}{b}} \right ){b}^{3}+24\, \left ( \arcsin \left ( dx+c \right ) \right ) ^{2}\sin \left ({\frac{a+b\arcsin \left ( dx+c \right ) }{b}}-{\frac{a}{b}} \right ) a{b}^{2}+20\, \left ( \arcsin \left ( dx+c \right ) \right ) ^{2}\cos \left ({\frac{a+b\arcsin \left ( dx+c \right ) }{b}}-{\frac{a}{b}} \right ){b}^{3}+24\,\arcsin \left ( dx+c \right ) \sin \left ({\frac{a+b\arcsin \left ( dx+c \right ) }{b}}-{\frac{a}{b}} \right ){a}^{2}b-30\,\arcsin \left ( dx+c \right ) \sin \left ({\frac{a+b\arcsin \left ( dx+c \right ) }{b}}-{\frac{a}{b}} \right ){b}^{3}+40\,\arcsin \left ( dx+c \right ) \cos \left ({\frac{a+b\arcsin \left ( dx+c \right ) }{b}}-{\frac{a}{b}} \right ) a{b}^{2}+8\,\sin \left ({\frac{a+b\arcsin \left ( dx+c \right ) }{b}}-{\frac{a}{b}} \right ){a}^{3}-30\,\sin \left ({\frac{a+b\arcsin \left ( dx+c \right ) }{b}}-{\frac{a}{b}} \right ) a{b}^{2}+20\,\cos \left ({\frac{a+b\arcsin \left ( dx+c \right ) }{b}}-{\frac{a}{b}} \right ){a}^{2}b \right ){\frac{1}{\sqrt{a+b\arcsin \left ( dx+c \right ) }}}} \end{align*}

Verification of antiderivative is not currently implemented for this CAS.

[In]

int((a+b*arcsin(d*x+c))^(5/2),x)

[Out]

1/8/d/(a+b*arcsin(d*x+c))^(1/2)*(15*(1/b)^(1/2)*Pi^(1/2)*2^(1/2)*(a+b*arcsin(d*x+c))^(1/2)*cos(a/b)*FresnelS(2
^(1/2)/Pi^(1/2)/(1/b)^(1/2)*(a+b*arcsin(d*x+c))^(1/2)/b)*b^3-15*(1/b)^(1/2)*Pi^(1/2)*2^(1/2)*(a+b*arcsin(d*x+c
))^(1/2)*sin(a/b)*FresnelC(2^(1/2)/Pi^(1/2)/(1/b)^(1/2)*(a+b*arcsin(d*x+c))^(1/2)/b)*b^3+8*arcsin(d*x+c)^3*sin
((a+b*arcsin(d*x+c))/b-a/b)*b^3+24*arcsin(d*x+c)^2*sin((a+b*arcsin(d*x+c))/b-a/b)*a*b^2+20*arcsin(d*x+c)^2*cos
((a+b*arcsin(d*x+c))/b-a/b)*b^3+24*arcsin(d*x+c)*sin((a+b*arcsin(d*x+c))/b-a/b)*a^2*b-30*arcsin(d*x+c)*sin((a+
b*arcsin(d*x+c))/b-a/b)*b^3+40*arcsin(d*x+c)*cos((a+b*arcsin(d*x+c))/b-a/b)*a*b^2+8*sin((a+b*arcsin(d*x+c))/b-
a/b)*a^3-30*sin((a+b*arcsin(d*x+c))/b-a/b)*a*b^2+20*cos((a+b*arcsin(d*x+c))/b-a/b)*a^2*b)

________________________________________________________________________________________

Maxima [F]  time = 0., size = 0, normalized size = 0. \begin{align*} \int{\left (b \arcsin \left (d x + c\right ) + a\right )}^{\frac{5}{2}}\,{d x} \end{align*}

Verification of antiderivative is not currently implemented for this CAS.

[In]

integrate((a+b*arcsin(d*x+c))^(5/2),x, algorithm="maxima")

[Out]

integrate((b*arcsin(d*x + c) + a)^(5/2), x)

________________________________________________________________________________________

Fricas [F(-2)]  time = 0., size = 0, normalized size = 0. \begin{align*} \text{Exception raised: UnboundLocalError} \end{align*}

Verification of antiderivative is not currently implemented for this CAS.

[In]

integrate((a+b*arcsin(d*x+c))^(5/2),x, algorithm="fricas")

[Out]

Exception raised: UnboundLocalError

________________________________________________________________________________________

Sympy [F(-1)]  time = 0., size = 0, normalized size = 0. \begin{align*} \text{Timed out} \end{align*}

Verification of antiderivative is not currently implemented for this CAS.

[In]

integrate((a+b*asin(d*x+c))**(5/2),x)

[Out]

Timed out

________________________________________________________________________________________

Giac [B]  time = 2.47833, size = 1766, normalized size = 8.66 \begin{align*} \text{result too large to display} \end{align*}

Verification of antiderivative is not currently implemented for this CAS.

[In]

integrate((a+b*arcsin(d*x+c))^(5/2),x, algorithm="giac")

[Out]

-1/2*sqrt(2)*sqrt(pi)*a^2*b^3*i*erf(-1/2*sqrt(2)*sqrt(b*arcsin(d*x + c) + a)*i/sqrt(abs(b)) - 1/2*sqrt(2)*sqrt
(b*arcsin(d*x + c) + a)*sqrt(abs(b))/b)*e^(a*i/b)/((b^3*i/sqrt(abs(b)) + b^2*sqrt(abs(b)))*d) - 1/2*sqrt(2)*sq
rt(pi)*a^2*b^3*i*erf(1/2*sqrt(2)*sqrt(b*arcsin(d*x + c) + a)*i/sqrt(abs(b)) - 1/2*sqrt(2)*sqrt(b*arcsin(d*x +
c) + a)*sqrt(abs(b))/b)*e^(-a*i/b)/((b^3*i/sqrt(abs(b)) - b^2*sqrt(abs(b)))*d) + 3/4*sqrt(2)*sqrt(pi)*a*b^4*er
f(-1/2*sqrt(2)*sqrt(b*arcsin(d*x + c) + a)*i/sqrt(abs(b)) - 1/2*sqrt(2)*sqrt(b*arcsin(d*x + c) + a)*sqrt(abs(b
))/b)*e^(a*i/b)/((b^3*i/sqrt(abs(b)) + b^2*sqrt(abs(b)))*d) + 1/2*sqrt(2)*sqrt(pi)*a^2*b^2*i*erf(-1/2*sqrt(2)*
sqrt(b*arcsin(d*x + c) + a)*i/sqrt(abs(b)) - 1/2*sqrt(2)*sqrt(b*arcsin(d*x + c) + a)*sqrt(abs(b))/b)*e^(a*i/b)
/((b^2*i/sqrt(abs(b)) + b*sqrt(abs(b)))*d) - 15/16*sqrt(2)*sqrt(pi)*b^4*i*erf(-1/2*sqrt(2)*sqrt(b*arcsin(d*x +
c) + a)*i/sqrt(abs(b)) - 1/2*sqrt(2)*sqrt(b*arcsin(d*x + c) + a)*sqrt(abs(b))/b)*e^(a*i/b)/((b^2*i/sqrt(abs(b
)) + b*sqrt(abs(b)))*d) - 3/4*sqrt(2)*sqrt(pi)*a*b^4*erf(1/2*sqrt(2)*sqrt(b*arcsin(d*x + c) + a)*i/sqrt(abs(b)
) - 1/2*sqrt(2)*sqrt(b*arcsin(d*x + c) + a)*sqrt(abs(b))/b)*e^(-a*i/b)/((b^3*i/sqrt(abs(b)) - b^2*sqrt(abs(b))
)*d) + 1/2*sqrt(2)*sqrt(pi)*a^2*b^2*i*erf(1/2*sqrt(2)*sqrt(b*arcsin(d*x + c) + a)*i/sqrt(abs(b)) - 1/2*sqrt(2)
*sqrt(b*arcsin(d*x + c) + a)*sqrt(abs(b))/b)*e^(-a*i/b)/((b^2*i/sqrt(abs(b)) - b*sqrt(abs(b)))*d) - 15/16*sqrt
(2)*sqrt(pi)*b^4*i*erf(1/2*sqrt(2)*sqrt(b*arcsin(d*x + c) + a)*i/sqrt(abs(b)) - 1/2*sqrt(2)*sqrt(b*arcsin(d*x
+ c) + a)*sqrt(abs(b))/b)*e^(-a*i/b)/((b^2*i/sqrt(abs(b)) - b*sqrt(abs(b)))*d) - 1/2*sqrt(b*arcsin(d*x + c) +
a)*b^2*i*arcsin(d*x + c)^2*e^(i*arcsin(d*x + c))/d + 1/2*sqrt(b*arcsin(d*x + c) + a)*b^2*i*arcsin(d*x + c)^2*e
^(-i*arcsin(d*x + c))/d - 3/4*sqrt(2)*sqrt(pi)*a*b^3*erf(-1/2*sqrt(2)*sqrt(b*arcsin(d*x + c) + a)*i/sqrt(abs(b
)) - 1/2*sqrt(2)*sqrt(b*arcsin(d*x + c) + a)*sqrt(abs(b))/b)*e^(a*i/b)/((b^2*i/sqrt(abs(b)) + b*sqrt(abs(b)))*
d) + 3/4*sqrt(2)*sqrt(pi)*a*b^3*erf(1/2*sqrt(2)*sqrt(b*arcsin(d*x + c) + a)*i/sqrt(abs(b)) - 1/2*sqrt(2)*sqrt(
b*arcsin(d*x + c) + a)*sqrt(abs(b))/b)*e^(-a*i/b)/((b^2*i/sqrt(abs(b)) - b*sqrt(abs(b)))*d) - sqrt(b*arcsin(d*
x + c) + a)*a*b*i*arcsin(d*x + c)*e^(i*arcsin(d*x + c))/d + sqrt(b*arcsin(d*x + c) + a)*a*b*i*arcsin(d*x + c)*
e^(-i*arcsin(d*x + c))/d - 1/2*sqrt(b*arcsin(d*x + c) + a)*a^2*i*e^(i*arcsin(d*x + c))/d + 15/8*sqrt(b*arcsin(
d*x + c) + a)*b^2*i*e^(i*arcsin(d*x + c))/d + 5/4*sqrt(b*arcsin(d*x + c) + a)*b^2*arcsin(d*x + c)*e^(i*arcsin(
d*x + c))/d + 1/2*sqrt(b*arcsin(d*x + c) + a)*a^2*i*e^(-i*arcsin(d*x + c))/d - 15/8*sqrt(b*arcsin(d*x + c) + a
)*b^2*i*e^(-i*arcsin(d*x + c))/d + 5/4*sqrt(b*arcsin(d*x + c) + a)*b^2*arcsin(d*x + c)*e^(-i*arcsin(d*x + c))/
d + 5/4*sqrt(b*arcsin(d*x + c) + a)*a*b*e^(i*arcsin(d*x + c))/d + 5/4*sqrt(b*arcsin(d*x + c) + a)*a*b*e^(-i*ar
csin(d*x + c))/d