3.202 \(\int (a+b \sinh ^{-1}(c+d x))^{7/2} \, dx\)

Optimal. Leaf size=216 \[ \frac{105 \sqrt{\pi } b^{7/2} e^{a/b} \text{Erf}\left (\frac{\sqrt{a+b \sinh ^{-1}(c+d x)}}{\sqrt{b}}\right )}{32 d}+\frac{105 \sqrt{\pi } b^{7/2} e^{-\frac{a}{b}} \text{Erfi}\left (\frac{\sqrt{a+b \sinh ^{-1}(c+d x)}}{\sqrt{b}}\right )}{32 d}-\frac{105 b^3 \sqrt{(c+d x)^2+1} \sqrt{a+b \sinh ^{-1}(c+d x)}}{8 d}+\frac{35 b^2 (c+d x) \left (a+b \sinh ^{-1}(c+d x)\right )^{3/2}}{4 d}-\frac{7 b \sqrt{(c+d x)^2+1} \left (a+b \sinh ^{-1}(c+d x)\right )^{5/2}}{2 d}+\frac{(c+d x) \left (a+b \sinh ^{-1}(c+d x)\right )^{7/2}}{d} \]

[Out]

(-105*b^3*Sqrt[1 + (c + d*x)^2]*Sqrt[a + b*ArcSinh[c + d*x]])/(8*d) + (35*b^2*(c + d*x)*(a + b*ArcSinh[c + d*x
])^(3/2))/(4*d) - (7*b*Sqrt[1 + (c + d*x)^2]*(a + b*ArcSinh[c + d*x])^(5/2))/(2*d) + ((c + d*x)*(a + b*ArcSinh
[c + d*x])^(7/2))/d + (105*b^(7/2)*E^(a/b)*Sqrt[Pi]*Erf[Sqrt[a + b*ArcSinh[c + d*x]]/Sqrt[b]])/(32*d) + (105*b
^(7/2)*Sqrt[Pi]*Erfi[Sqrt[a + b*ArcSinh[c + d*x]]/Sqrt[b]])/(32*d*E^(a/b))

________________________________________________________________________________________

Rubi [A]  time = 0.418628, antiderivative size = 216, normalized size of antiderivative = 1., number of steps used = 11, number of rules used = 8, integrand size = 14, \(\frac{\text{number of rules}}{\text{integrand size}}\) = 0.571, Rules used = {5863, 5653, 5717, 5657, 3307, 2180, 2205, 2204} \[ \frac{105 \sqrt{\pi } b^{7/2} e^{a/b} \text{Erf}\left (\frac{\sqrt{a+b \sinh ^{-1}(c+d x)}}{\sqrt{b}}\right )}{32 d}+\frac{105 \sqrt{\pi } b^{7/2} e^{-\frac{a}{b}} \text{Erfi}\left (\frac{\sqrt{a+b \sinh ^{-1}(c+d x)}}{\sqrt{b}}\right )}{32 d}-\frac{105 b^3 \sqrt{(c+d x)^2+1} \sqrt{a+b \sinh ^{-1}(c+d x)}}{8 d}+\frac{35 b^2 (c+d x) \left (a+b \sinh ^{-1}(c+d x)\right )^{3/2}}{4 d}-\frac{7 b \sqrt{(c+d x)^2+1} \left (a+b \sinh ^{-1}(c+d x)\right )^{5/2}}{2 d}+\frac{(c+d x) \left (a+b \sinh ^{-1}(c+d x)\right )^{7/2}}{d} \]

Antiderivative was successfully verified.

[In]

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

[Out]

(-105*b^3*Sqrt[1 + (c + d*x)^2]*Sqrt[a + b*ArcSinh[c + d*x]])/(8*d) + (35*b^2*(c + d*x)*(a + b*ArcSinh[c + d*x
])^(3/2))/(4*d) - (7*b*Sqrt[1 + (c + d*x)^2]*(a + b*ArcSinh[c + d*x])^(5/2))/(2*d) + ((c + d*x)*(a + b*ArcSinh
[c + d*x])^(7/2))/d + (105*b^(7/2)*E^(a/b)*Sqrt[Pi]*Erf[Sqrt[a + b*ArcSinh[c + d*x]]/Sqrt[b]])/(32*d) + (105*b
^(7/2)*Sqrt[Pi]*Erfi[Sqrt[a + b*ArcSinh[c + d*x]]/Sqrt[b]])/(32*d*E^(a/b))

Rule 5863

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

Rule 5653

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

Rule 5717

Int[((a_.) + ArcSinh[(c_.)*(x_)]*(b_.))^(n_.)*(x_)*((d_) + (e_.)*(x_)^2)^(p_.), x_Symbol] :> Simp[((d + e*x^2)
^(p + 1)*(a + b*ArcSinh[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*ArcSinh[c*x])^(n - 1), x], x] /; FreeQ[{a,
b, c, d, e, p}, x] && EqQ[e, c^2*d] && GtQ[n, 0] && NeQ[p, -1]

Rule 5657

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

Rule 3307

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

Rule 2180

Int[(F_)^((g_.)*((e_.) + (f_.)*(x_)))/Sqrt[(c_.) + (d_.)*(x_)], x_Symbol] :> Dist[2/d, Subst[Int[F^(g*(e - (c*
f)/d) + (f*g*x^2)/d), x], x, Sqrt[c + d*x]], x] /; FreeQ[{F, c, d, e, f, g}, x] &&  !$UseGamma === True

Rule 2205

Int[(F_)^((a_.) + (b_.)*((c_.) + (d_.)*(x_))^2), x_Symbol] :> Simp[(F^a*Sqrt[Pi]*Erf[(c + d*x)*Rt[-(b*Log[F]),
2]])/(2*d*Rt[-(b*Log[F]), 2]), x] /; FreeQ[{F, a, b, c, d}, x] && NegQ[b]

Rule 2204

Int[(F_)^((a_.) + (b_.)*((c_.) + (d_.)*(x_))^2), x_Symbol] :> Simp[(F^a*Sqrt[Pi]*Erfi[(c + d*x)*Rt[b*Log[F], 2
]])/(2*d*Rt[b*Log[F], 2]), x] /; FreeQ[{F, a, b, c, d}, x] && PosQ[b]

Rubi steps

\begin{align*} \int \left (a+b \sinh ^{-1}(c+d x)\right )^{7/2} \, dx &=\frac{\operatorname{Subst}\left (\int \left (a+b \sinh ^{-1}(x)\right )^{7/2} \, dx,x,c+d x\right )}{d}\\ &=\frac{(c+d x) \left (a+b \sinh ^{-1}(c+d x)\right )^{7/2}}{d}-\frac{(7 b) \operatorname{Subst}\left (\int \frac{x \left (a+b \sinh ^{-1}(x)\right )^{5/2}}{\sqrt{1+x^2}} \, dx,x,c+d x\right )}{2 d}\\ &=-\frac{7 b \sqrt{1+(c+d x)^2} \left (a+b \sinh ^{-1}(c+d x)\right )^{5/2}}{2 d}+\frac{(c+d x) \left (a+b \sinh ^{-1}(c+d x)\right )^{7/2}}{d}+\frac{\left (35 b^2\right ) \operatorname{Subst}\left (\int \left (a+b \sinh ^{-1}(x)\right )^{3/2} \, dx,x,c+d x\right )}{4 d}\\ &=\frac{35 b^2 (c+d x) \left (a+b \sinh ^{-1}(c+d x)\right )^{3/2}}{4 d}-\frac{7 b \sqrt{1+(c+d x)^2} \left (a+b \sinh ^{-1}(c+d x)\right )^{5/2}}{2 d}+\frac{(c+d x) \left (a+b \sinh ^{-1}(c+d x)\right )^{7/2}}{d}-\frac{\left (105 b^3\right ) \operatorname{Subst}\left (\int \frac{x \sqrt{a+b \sinh ^{-1}(x)}}{\sqrt{1+x^2}} \, dx,x,c+d x\right )}{8 d}\\ &=-\frac{105 b^3 \sqrt{1+(c+d x)^2} \sqrt{a+b \sinh ^{-1}(c+d x)}}{8 d}+\frac{35 b^2 (c+d x) \left (a+b \sinh ^{-1}(c+d x)\right )^{3/2}}{4 d}-\frac{7 b \sqrt{1+(c+d x)^2} \left (a+b \sinh ^{-1}(c+d x)\right )^{5/2}}{2 d}+\frac{(c+d x) \left (a+b \sinh ^{-1}(c+d x)\right )^{7/2}}{d}+\frac{\left (105 b^4\right ) \operatorname{Subst}\left (\int \frac{1}{\sqrt{a+b \sinh ^{-1}(x)}} \, dx,x,c+d x\right )}{16 d}\\ &=-\frac{105 b^3 \sqrt{1+(c+d x)^2} \sqrt{a+b \sinh ^{-1}(c+d x)}}{8 d}+\frac{35 b^2 (c+d x) \left (a+b \sinh ^{-1}(c+d x)\right )^{3/2}}{4 d}-\frac{7 b \sqrt{1+(c+d x)^2} \left (a+b \sinh ^{-1}(c+d x)\right )^{5/2}}{2 d}+\frac{(c+d x) \left (a+b \sinh ^{-1}(c+d x)\right )^{7/2}}{d}+\frac{\left (105 b^3\right ) \operatorname{Subst}\left (\int \frac{\cosh \left (\frac{a}{b}-\frac{x}{b}\right )}{\sqrt{x}} \, dx,x,a+b \sinh ^{-1}(c+d x)\right )}{16 d}\\ &=-\frac{105 b^3 \sqrt{1+(c+d x)^2} \sqrt{a+b \sinh ^{-1}(c+d x)}}{8 d}+\frac{35 b^2 (c+d x) \left (a+b \sinh ^{-1}(c+d x)\right )^{3/2}}{4 d}-\frac{7 b \sqrt{1+(c+d x)^2} \left (a+b \sinh ^{-1}(c+d x)\right )^{5/2}}{2 d}+\frac{(c+d x) \left (a+b \sinh ^{-1}(c+d x)\right )^{7/2}}{d}+\frac{\left (105 b^3\right ) \operatorname{Subst}\left (\int \frac{e^{-i \left (\frac{i a}{b}-\frac{i x}{b}\right )}}{\sqrt{x}} \, dx,x,a+b \sinh ^{-1}(c+d x)\right )}{32 d}+\frac{\left (105 b^3\right ) \operatorname{Subst}\left (\int \frac{e^{i \left (\frac{i a}{b}-\frac{i x}{b}\right )}}{\sqrt{x}} \, dx,x,a+b \sinh ^{-1}(c+d x)\right )}{32 d}\\ &=-\frac{105 b^3 \sqrt{1+(c+d x)^2} \sqrt{a+b \sinh ^{-1}(c+d x)}}{8 d}+\frac{35 b^2 (c+d x) \left (a+b \sinh ^{-1}(c+d x)\right )^{3/2}}{4 d}-\frac{7 b \sqrt{1+(c+d x)^2} \left (a+b \sinh ^{-1}(c+d x)\right )^{5/2}}{2 d}+\frac{(c+d x) \left (a+b \sinh ^{-1}(c+d x)\right )^{7/2}}{d}+\frac{\left (105 b^3\right ) \operatorname{Subst}\left (\int e^{\frac{a}{b}-\frac{x^2}{b}} \, dx,x,\sqrt{a+b \sinh ^{-1}(c+d x)}\right )}{16 d}+\frac{\left (105 b^3\right ) \operatorname{Subst}\left (\int e^{-\frac{a}{b}+\frac{x^2}{b}} \, dx,x,\sqrt{a+b \sinh ^{-1}(c+d x)}\right )}{16 d}\\ &=-\frac{105 b^3 \sqrt{1+(c+d x)^2} \sqrt{a+b \sinh ^{-1}(c+d x)}}{8 d}+\frac{35 b^2 (c+d x) \left (a+b \sinh ^{-1}(c+d x)\right )^{3/2}}{4 d}-\frac{7 b \sqrt{1+(c+d x)^2} \left (a+b \sinh ^{-1}(c+d x)\right )^{5/2}}{2 d}+\frac{(c+d x) \left (a+b \sinh ^{-1}(c+d x)\right )^{7/2}}{d}+\frac{105 b^{7/2} e^{a/b} \sqrt{\pi } \text{erf}\left (\frac{\sqrt{a+b \sinh ^{-1}(c+d x)}}{\sqrt{b}}\right )}{32 d}+\frac{105 b^{7/2} e^{-\frac{a}{b}} \sqrt{\pi } \text{erfi}\left (\frac{\sqrt{a+b \sinh ^{-1}(c+d x)}}{\sqrt{b}}\right )}{32 d}\\ \end{align*}

Mathematica [B]  time = 4.53561, size = 698, normalized size = 3.23 \[ \frac{16 a^3 e^{-\frac{a}{b}} \sqrt{a+b \sinh ^{-1}(c+d x)} \left (\frac{\text{Gamma}\left (\frac{3}{2},-\frac{a+b \sinh ^{-1}(c+d x)}{b}\right )}{\sqrt{-\frac{a+b \sinh ^{-1}(c+d x)}{b}}}-\frac{e^{\frac{2 a}{b}} \text{Gamma}\left (\frac{3}{2},\frac{a}{b}+\sinh ^{-1}(c+d x)\right )}{\sqrt{\frac{a}{b}+\sinh ^{-1}(c+d x)}}\right )+6 a \sqrt{b} \left (\sqrt{\pi } \left (4 a^2-12 a b+15 b^2\right ) \left (\sinh \left (\frac{a}{b}\right )+\cosh \left (\frac{a}{b}\right )\right ) \text{Erf}\left (\frac{\sqrt{a+b \sinh ^{-1}(c+d x)}}{\sqrt{b}}\right )+\sqrt{\pi } \left (4 a^2+12 a b+15 b^2\right ) \left (\sinh \left (\frac{a}{b}\right )-\cosh \left (\frac{a}{b}\right )\right ) \text{Erfi}\left (\frac{\sqrt{a+b \sinh ^{-1}(c+d x)}}{\sqrt{b}}\right )+4 \sqrt{b} \sqrt{a+b \sinh ^{-1}(c+d x)} \left (2 \sqrt{(c+d x)^2+1} \left (a-5 b \sinh ^{-1}(c+d x)\right )+b (c+d x) \left (4 \sinh ^{-1}(c+d x)^2+15\right )\right )\right )+\sqrt{b} \left (\sqrt{\pi } \left (36 a^2 b-8 a^3-90 a b^2+105 b^3\right ) \left (\sinh \left (\frac{a}{b}\right )+\cosh \left (\frac{a}{b}\right )\right ) \text{Erf}\left (\frac{\sqrt{a+b \sinh ^{-1}(c+d x)}}{\sqrt{b}}\right )+\sqrt{\pi } \left (36 a^2 b+8 a^3+90 a b^2+105 b^3\right ) \left (\cosh \left (\frac{a}{b}\right )-\sinh \left (\frac{a}{b}\right )\right ) \text{Erfi}\left (\frac{\sqrt{a+b \sinh ^{-1}(c+d x)}}{\sqrt{b}}\right )+4 \sqrt{b} \sqrt{a+b \sinh ^{-1}(c+d x)} \left (\sqrt{(c+d x)^2+1} \left (-4 a^2+4 a b \sinh ^{-1}(c+d x)-7 b^2 \left (4 \sinh ^{-1}(c+d x)^2+15\right )\right )+2 b (c+d x) \left (-10 a+4 b \sinh ^{-1}(c+d x)^3+35 b \sinh ^{-1}(c+d x)\right )\right )\right )+12 a^2 \sqrt{b} \left (\sqrt{\pi } (3 b-2 a) \left (\sinh \left (\frac{a}{b}\right )+\cosh \left (\frac{a}{b}\right )\right ) \text{Erf}\left (\frac{\sqrt{a+b \sinh ^{-1}(c+d x)}}{\sqrt{b}}\right )+\sqrt{\pi } (2 a+3 b) \left (\cosh \left (\frac{a}{b}\right )-\sinh \left (\frac{a}{b}\right )\right ) \text{Erfi}\left (\frac{\sqrt{a+b \sinh ^{-1}(c+d x)}}{\sqrt{b}}\right )+4 \sqrt{b} \left (2 (c+d x) \sinh ^{-1}(c+d x)-3 \sqrt{(c+d x)^2+1}\right ) \sqrt{a+b \sinh ^{-1}(c+d x)}\right )}{32 d} \]

Warning: Unable to verify antiderivative.

[In]

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

[Out]

((16*a^3*Sqrt[a + b*ArcSinh[c + d*x]]*(-((E^((2*a)/b)*Gamma[3/2, a/b + ArcSinh[c + d*x]])/Sqrt[a/b + ArcSinh[c
+ d*x]]) + Gamma[3/2, -((a + b*ArcSinh[c + d*x])/b)]/Sqrt[-((a + b*ArcSinh[c + d*x])/b)]))/E^(a/b) + 12*a^2*S
qrt[b]*(4*Sqrt[b]*Sqrt[a + b*ArcSinh[c + d*x]]*(-3*Sqrt[1 + (c + d*x)^2] + 2*(c + d*x)*ArcSinh[c + d*x]) + (2*
a + 3*b)*Sqrt[Pi]*Erfi[Sqrt[a + b*ArcSinh[c + d*x]]/Sqrt[b]]*(Cosh[a/b] - Sinh[a/b]) + (-2*a + 3*b)*Sqrt[Pi]*E
rf[Sqrt[a + b*ArcSinh[c + d*x]]/Sqrt[b]]*(Cosh[a/b] + Sinh[a/b])) + 6*a*Sqrt[b]*(4*Sqrt[b]*Sqrt[a + b*ArcSinh[
c + d*x]]*(2*Sqrt[1 + (c + d*x)^2]*(a - 5*b*ArcSinh[c + d*x]) + b*(c + d*x)*(15 + 4*ArcSinh[c + d*x]^2)) + (4*
a^2 + 12*a*b + 15*b^2)*Sqrt[Pi]*Erfi[Sqrt[a + b*ArcSinh[c + d*x]]/Sqrt[b]]*(-Cosh[a/b] + Sinh[a/b]) + (4*a^2 -
12*a*b + 15*b^2)*Sqrt[Pi]*Erf[Sqrt[a + b*ArcSinh[c + d*x]]/Sqrt[b]]*(Cosh[a/b] + Sinh[a/b])) + Sqrt[b]*(4*Sqr
t[b]*Sqrt[a + b*ArcSinh[c + d*x]]*(2*b*(c + d*x)*(-10*a + 35*b*ArcSinh[c + d*x] + 4*b*ArcSinh[c + d*x]^3) + Sq
rt[1 + (c + d*x)^2]*(-4*a^2 + 4*a*b*ArcSinh[c + d*x] - 7*b^2*(15 + 4*ArcSinh[c + d*x]^2))) + (8*a^3 + 36*a^2*b
+ 90*a*b^2 + 105*b^3)*Sqrt[Pi]*Erfi[Sqrt[a + b*ArcSinh[c + d*x]]/Sqrt[b]]*(Cosh[a/b] - Sinh[a/b]) + (-8*a^3 +
36*a^2*b - 90*a*b^2 + 105*b^3)*Sqrt[Pi]*Erf[Sqrt[a + b*ArcSinh[c + d*x]]/Sqrt[b]]*(Cosh[a/b] + Sinh[a/b])))/(
32*d)

________________________________________________________________________________________

Maple [F]  time = 0.082, size = 0, normalized size = 0. \begin{align*} \int \left ( a+b{\it Arcsinh} \left ( dx+c \right ) \right ) ^{{\frac{7}{2}}}\, dx \end{align*}

Verification of antiderivative is not currently implemented for this CAS.

[In]

int((a+b*arcsinh(d*x+c))^(7/2),x)

[Out]

int((a+b*arcsinh(d*x+c))^(7/2),x)

________________________________________________________________________________________

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

Verification of antiderivative is not currently implemented for this CAS.

[In]

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

[Out]

integrate((b*arcsinh(d*x + c) + a)^(7/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*arcsinh(d*x+c))^(7/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*asinh(d*x+c))**(7/2),x)

[Out]

Timed out

________________________________________________________________________________________

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

Verification of antiderivative is not currently implemented for this CAS.

[In]

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

[Out]

integrate((b*arcsinh(d*x + c) + a)^(7/2), x)