\(\int (a+b \text {arcsinh}(c+d x))^4 \, dx\) [150]

Optimal result

Integrand size = 12, antiderivative size = 115 \[ \int (a+b \text {arcsinh}(c+d x))^4 \, dx=24 b^4 x-\frac {24 b^3 \sqrt {1+(c+d x)^2} (a+b \text {arcsinh}(c+d x))}{d}+\frac {12 b^2 (c+d x) (a+b \text {arcsinh}(c+d x))^2}{d}-\frac {4 b \sqrt {1+(c+d x)^2} (a+b \text {arcsinh}(c+d x))^3}{d}+\frac {(c+d x) (a+b \text {arcsinh}(c+d x))^4}{d} \]

[Out]

24*b^4*x+12*b^2*(d*x+c)*(a+b*arcsinh(d*x+c))^2/d+(d*x+c)*(a+b*arcsinh(d*x+c))^4/d-24*b^3*(a+b*arcsinh(d*x+c))*
(1+(d*x+c)^2)^(1/2)/d-4*b*(a+b*arcsinh(d*x+c))^3*(1+(d*x+c)^2)^(1/2)/d

Rubi [A] (verified)

Time = 0.11 (sec) , antiderivative size = 115, normalized size of antiderivative = 1.00, number of steps used = 6, number of rules used = 4, \(\frac {\text {number of rules}}{\text {integrand size}}\) = 0.333, Rules used = {5858, 5772, 5798, 8} \[ \int (a+b \text {arcsinh}(c+d x))^4 \, dx=-\frac {24 b^3 \sqrt {(c+d x)^2+1} (a+b \text {arcsinh}(c+d x))}{d}+\frac {12 b^2 (c+d x) (a+b \text {arcsinh}(c+d x))^2}{d}-\frac {4 b \sqrt {(c+d x)^2+1} (a+b \text {arcsinh}(c+d x))^3}{d}+\frac {(c+d x) (a+b \text {arcsinh}(c+d x))^4}{d}+24 b^4 x \]

[In]

Int[(a + b*ArcSinh[c + d*x])^4,x]

[Out]

24*b^4*x - (24*b^3*Sqrt[1 + (c + d*x)^2]*(a + b*ArcSinh[c + d*x]))/d + (12*b^2*(c + d*x)*(a + b*ArcSinh[c + d*
x])^2)/d - (4*b*Sqrt[1 + (c + d*x)^2]*(a + b*ArcSinh[c + d*x])^3)/d + ((c + d*x)*(a + b*ArcSinh[c + d*x])^4)/d

Rule 8

Int[a_, x_Symbol] :> Simp[a*x, x] /; FreeQ[a, x]

Rule 5772

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 5798

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/(2*c*(p + 1)))*Simp[(d + e*x^2)^p/(1 + c^2*x^2)
^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 5858

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]

Rubi steps \begin{align*} \text {integral}& = \frac {\text {Subst}\left (\int (a+b \text {arcsinh}(x))^4 \, dx,x,c+d x\right )}{d} \\ & = \frac {(c+d x) (a+b \text {arcsinh}(c+d x))^4}{d}-\frac {(4 b) \text {Subst}\left (\int \frac {x (a+b \text {arcsinh}(x))^3}{\sqrt {1+x^2}} \, dx,x,c+d x\right )}{d} \\ & = -\frac {4 b \sqrt {1+(c+d x)^2} (a+b \text {arcsinh}(c+d x))^3}{d}+\frac {(c+d x) (a+b \text {arcsinh}(c+d x))^4}{d}+\frac {\left (12 b^2\right ) \text {Subst}\left (\int (a+b \text {arcsinh}(x))^2 \, dx,x,c+d x\right )}{d} \\ & = \frac {12 b^2 (c+d x) (a+b \text {arcsinh}(c+d x))^2}{d}-\frac {4 b \sqrt {1+(c+d x)^2} (a+b \text {arcsinh}(c+d x))^3}{d}+\frac {(c+d x) (a+b \text {arcsinh}(c+d x))^4}{d}-\frac {\left (24 b^3\right ) \text {Subst}\left (\int \frac {x (a+b \text {arcsinh}(x))}{\sqrt {1+x^2}} \, dx,x,c+d x\right )}{d} \\ & = -\frac {24 b^3 \sqrt {1+(c+d x)^2} (a+b \text {arcsinh}(c+d x))}{d}+\frac {12 b^2 (c+d x) (a+b \text {arcsinh}(c+d x))^2}{d}-\frac {4 b \sqrt {1+(c+d x)^2} (a+b \text {arcsinh}(c+d x))^3}{d}+\frac {(c+d x) (a+b \text {arcsinh}(c+d x))^4}{d}+\frac {\left (24 b^4\right ) \text {Subst}(\int 1 \, dx,x,c+d x)}{d} \\ & = 24 b^4 x-\frac {24 b^3 \sqrt {1+(c+d x)^2} (a+b \text {arcsinh}(c+d x))}{d}+\frac {12 b^2 (c+d x) (a+b \text {arcsinh}(c+d x))^2}{d}-\frac {4 b \sqrt {1+(c+d x)^2} (a+b \text {arcsinh}(c+d x))^3}{d}+\frac {(c+d x) (a+b \text {arcsinh}(c+d x))^4}{d} \\ \end{align*}

Mathematica [A] (verified)

Time = 0.23 (sec) , antiderivative size = 226, normalized size of antiderivative = 1.97 \[ \int (a+b \text {arcsinh}(c+d x))^4 \, dx=\frac {\left (a^4+12 a^2 b^2+24 b^4\right ) (c+d x)-4 a b \left (a^2+6 b^2\right ) \sqrt {1+(c+d x)^2}-4 b \left (-a^3 (c+d x)-6 a b^2 (c+d x)+3 a^2 b \sqrt {1+(c+d x)^2}+6 b^3 \sqrt {1+(c+d x)^2}\right ) \text {arcsinh}(c+d x)+6 b^2 \left (a^2 (c+d x)+2 b^2 (c+d x)-2 a b \sqrt {1+(c+d x)^2}\right ) \text {arcsinh}(c+d x)^2-4 b^3 \left (-a (c+d x)+b \sqrt {1+(c+d x)^2}\right ) \text {arcsinh}(c+d x)^3+b^4 (c+d x) \text {arcsinh}(c+d x)^4}{d} \]

[In]

Integrate[(a + b*ArcSinh[c + d*x])^4,x]

[Out]

((a^4 + 12*a^2*b^2 + 24*b^4)*(c + d*x) - 4*a*b*(a^2 + 6*b^2)*Sqrt[1 + (c + d*x)^2] - 4*b*(-(a^3*(c + d*x)) - 6
*a*b^2*(c + d*x) + 3*a^2*b*Sqrt[1 + (c + d*x)^2] + 6*b^3*Sqrt[1 + (c + d*x)^2])*ArcSinh[c + d*x] + 6*b^2*(a^2*
(c + d*x) + 2*b^2*(c + d*x) - 2*a*b*Sqrt[1 + (c + d*x)^2])*ArcSinh[c + d*x]^2 - 4*b^3*(-(a*(c + d*x)) + b*Sqrt
[1 + (c + d*x)^2])*ArcSinh[c + d*x]^3 + b^4*(c + d*x)*ArcSinh[c + d*x]^4)/d

Maple [B] (verified)

Leaf count of result is larger than twice the leaf count of optimal. \(244\) vs. \(2(111)=222\).

Time = 0.17 (sec) , antiderivative size = 245, normalized size of antiderivative = 2.13

method	result	size
derivativedivides	\(\frac {\left (d x +c \right ) a^{4}+b^{4} \left (\left (d x +c \right ) \operatorname {arcsinh}\left (d x +c \right )^{4}-4 \operatorname {arcsinh}\left (d x +c \right )^{3} \sqrt {1+\left (d x +c \right )^{2}}+12 \left (d x +c \right ) \operatorname {arcsinh}\left (d x +c \right )^{2}-24 \,\operatorname {arcsinh}\left (d x +c \right ) \sqrt {1+\left (d x +c \right )^{2}}+24 d x +24 c \right )+4 a \,b^{3} \left (\left (d x +c \right ) \operatorname {arcsinh}\left (d x +c \right )^{3}-3 \operatorname {arcsinh}\left (d x +c \right )^{2} \sqrt {1+\left (d x +c \right )^{2}}+6 \left (d x +c \right ) \operatorname {arcsinh}\left (d x +c \right )-6 \sqrt {1+\left (d x +c \right )^{2}}\right )+6 a^{2} b^{2} \left (\left (d x +c \right ) \operatorname {arcsinh}\left (d x +c \right )^{2}-2 \,\operatorname {arcsinh}\left (d x +c \right ) \sqrt {1+\left (d x +c \right )^{2}}+2 d x +2 c \right )+4 b \,a^{3} \left (\left (d x +c \right ) \operatorname {arcsinh}\left (d x +c \right )-\sqrt {1+\left (d x +c \right )^{2}}\right )}{d}\)	\(245\)
default	\(\frac {\left (d x +c \right ) a^{4}+b^{4} \left (\left (d x +c \right ) \operatorname {arcsinh}\left (d x +c \right )^{4}-4 \operatorname {arcsinh}\left (d x +c \right )^{3} \sqrt {1+\left (d x +c \right )^{2}}+12 \left (d x +c \right ) \operatorname {arcsinh}\left (d x +c \right )^{2}-24 \,\operatorname {arcsinh}\left (d x +c \right ) \sqrt {1+\left (d x +c \right )^{2}}+24 d x +24 c \right )+4 a \,b^{3} \left (\left (d x +c \right ) \operatorname {arcsinh}\left (d x +c \right )^{3}-3 \operatorname {arcsinh}\left (d x +c \right )^{2} \sqrt {1+\left (d x +c \right )^{2}}+6 \left (d x +c \right ) \operatorname {arcsinh}\left (d x +c \right )-6 \sqrt {1+\left (d x +c \right )^{2}}\right )+6 a^{2} b^{2} \left (\left (d x +c \right ) \operatorname {arcsinh}\left (d x +c \right )^{2}-2 \,\operatorname {arcsinh}\left (d x +c \right ) \sqrt {1+\left (d x +c \right )^{2}}+2 d x +2 c \right )+4 b \,a^{3} \left (\left (d x +c \right ) \operatorname {arcsinh}\left (d x +c \right )-\sqrt {1+\left (d x +c \right )^{2}}\right )}{d}\)	\(245\)
parts	\(a^{4} x +\frac {b^{4} \left (\left (d x +c \right ) \operatorname {arcsinh}\left (d x +c \right )^{4}-4 \operatorname {arcsinh}\left (d x +c \right )^{3} \sqrt {1+\left (d x +c \right )^{2}}+12 \left (d x +c \right ) \operatorname {arcsinh}\left (d x +c \right )^{2}-24 \,\operatorname {arcsinh}\left (d x +c \right ) \sqrt {1+\left (d x +c \right )^{2}}+24 d x +24 c \right )}{d}+\frac {4 a \,b^{3} \left (\left (d x +c \right ) \operatorname {arcsinh}\left (d x +c \right )^{3}-3 \operatorname {arcsinh}\left (d x +c \right )^{2} \sqrt {1+\left (d x +c \right )^{2}}+6 \left (d x +c \right ) \operatorname {arcsinh}\left (d x +c \right )-6 \sqrt {1+\left (d x +c \right )^{2}}\right )}{d}+\frac {6 a^{2} b^{2} \left (\left (d x +c \right ) \operatorname {arcsinh}\left (d x +c \right )^{2}-2 \,\operatorname {arcsinh}\left (d x +c \right ) \sqrt {1+\left (d x +c \right )^{2}}+2 d x +2 c \right )}{d}+\frac {4 b \,a^{3} \left (\left (d x +c \right ) \operatorname {arcsinh}\left (d x +c \right )-\sqrt {1+\left (d x +c \right )^{2}}\right )}{d}\)	\(249\)

[In]

int((a+b*arcsinh(d*x+c))^4,x,method=_RETURNVERBOSE)

[Out]

1/d*((d*x+c)*a^4+b^4*((d*x+c)*arcsinh(d*x+c)^4-4*arcsinh(d*x+c)^3*(1+(d*x+c)^2)^(1/2)+12*(d*x+c)*arcsinh(d*x+c
)^2-24*arcsinh(d*x+c)*(1+(d*x+c)^2)^(1/2)+24*d*x+24*c)+4*a*b^3*((d*x+c)*arcsinh(d*x+c)^3-3*arcsinh(d*x+c)^2*(1
+(d*x+c)^2)^(1/2)+6*(d*x+c)*arcsinh(d*x+c)-6*(1+(d*x+c)^2)^(1/2))+6*a^2*b^2*((d*x+c)*arcsinh(d*x+c)^2-2*arcsin
h(d*x+c)*(1+(d*x+c)^2)^(1/2)+2*d*x+2*c)+4*b*a^3*((d*x+c)*arcsinh(d*x+c)-(1+(d*x+c)^2)^(1/2)))

Fricas [B] (verification not implemented)

Leaf count of result is larger than twice the leaf count of optimal. 344 vs. \(2 (111) = 222\).

Time = 0.28 (sec) , antiderivative size = 344, normalized size of antiderivative = 2.99 \[ \int (a+b \text {arcsinh}(c+d x))^4 \, dx=\frac {{\left (b^{4} d x + b^{4} c\right )} \log \left (d x + c + \sqrt {d^{2} x^{2} + 2 \, c d x + c^{2} + 1}\right )^{4} + 4 \, {\left (a b^{3} d x + a b^{3} c - \sqrt {d^{2} x^{2} + 2 \, c d x + c^{2} + 1} b^{4}\right )} \log \left (d x + c + \sqrt {d^{2} x^{2} + 2 \, c d x + c^{2} + 1}\right )^{3} + {\left (a^{4} + 12 \, a^{2} b^{2} + 24 \, b^{4}\right )} d x - 6 \, {\left (2 \, \sqrt {d^{2} x^{2} + 2 \, c d x + c^{2} + 1} a b^{3} - {\left (a^{2} b^{2} + 2 \, b^{4}\right )} d x - {\left (a^{2} b^{2} + 2 \, b^{4}\right )} c\right )} \log \left (d x + c + \sqrt {d^{2} x^{2} + 2 \, c d x + c^{2} + 1}\right )^{2} + 4 \, {\left ({\left (a^{3} b + 6 \, a b^{3}\right )} d x + {\left (a^{3} b + 6 \, a b^{3}\right )} c - 3 \, {\left (a^{2} b^{2} + 2 \, b^{4}\right )} \sqrt {d^{2} x^{2} + 2 \, c d x + c^{2} + 1}\right )} \log \left (d x + c + \sqrt {d^{2} x^{2} + 2 \, c d x + c^{2} + 1}\right ) - 4 \, {\left (a^{3} b + 6 \, a b^{3}\right )} \sqrt {d^{2} x^{2} + 2 \, c d x + c^{2} + 1}}{d} \]

[In]

integrate((a+b*arcsinh(d*x+c))^4,x, algorithm="fricas")

[Out]

((b^4*d*x + b^4*c)*log(d*x + c + sqrt(d^2*x^2 + 2*c*d*x + c^2 + 1))^4 + 4*(a*b^3*d*x + a*b^3*c - sqrt(d^2*x^2
+ 2*c*d*x + c^2 + 1)*b^4)*log(d*x + c + sqrt(d^2*x^2 + 2*c*d*x + c^2 + 1))^3 + (a^4 + 12*a^2*b^2 + 24*b^4)*d*x
- 6*(2*sqrt(d^2*x^2 + 2*c*d*x + c^2 + 1)*a*b^3 - (a^2*b^2 + 2*b^4)*d*x - (a^2*b^2 + 2*b^4)*c)*log(d*x + c + s
qrt(d^2*x^2 + 2*c*d*x + c^2 + 1))^2 + 4*((a^3*b + 6*a*b^3)*d*x + (a^3*b + 6*a*b^3)*c - 3*(a^2*b^2 + 2*b^4)*sqr
t(d^2*x^2 + 2*c*d*x + c^2 + 1))*log(d*x + c + sqrt(d^2*x^2 + 2*c*d*x + c^2 + 1)) - 4*(a^3*b + 6*a*b^3)*sqrt(d^
2*x^2 + 2*c*d*x + c^2 + 1))/d

Sympy [B] (verification not implemented)

Leaf count of result is larger than twice the leaf count of optimal. 444 vs. \(2 (105) = 210\).

Time = 0.27 (sec) , antiderivative size = 444, normalized size of antiderivative = 3.86 \[ \int (a+b \text {arcsinh}(c+d x))^4 \, dx=\begin {cases} a^{4} x + \frac {4 a^{3} b c \operatorname {asinh}{\left (c + d x \right )}}{d} + 4 a^{3} b x \operatorname {asinh}{\left (c + d x \right )} - \frac {4 a^{3} b \sqrt {c^{2} + 2 c d x + d^{2} x^{2} + 1}}{d} + \frac {6 a^{2} b^{2} c \operatorname {asinh}^{2}{\left (c + d x \right )}}{d} + 6 a^{2} b^{2} x \operatorname {asinh}^{2}{\left (c + d x \right )} + 12 a^{2} b^{2} x - \frac {12 a^{2} b^{2} \sqrt {c^{2} + 2 c d x + d^{2} x^{2} + 1} \operatorname {asinh}{\left (c + d x \right )}}{d} + \frac {4 a b^{3} c \operatorname {asinh}^{3}{\left (c + d x \right )}}{d} + \frac {24 a b^{3} c \operatorname {asinh}{\left (c + d x \right )}}{d} + 4 a b^{3} x \operatorname {asinh}^{3}{\left (c + d x \right )} + 24 a b^{3} x \operatorname {asinh}{\left (c + d x \right )} - \frac {12 a b^{3} \sqrt {c^{2} + 2 c d x + d^{2} x^{2} + 1} \operatorname {asinh}^{2}{\left (c + d x \right )}}{d} - \frac {24 a b^{3} \sqrt {c^{2} + 2 c d x + d^{2} x^{2} + 1}}{d} + \frac {b^{4} c \operatorname {asinh}^{4}{\left (c + d x \right )}}{d} + \frac {12 b^{4} c \operatorname {asinh}^{2}{\left (c + d x \right )}}{d} + b^{4} x \operatorname {asinh}^{4}{\left (c + d x \right )} + 12 b^{4} x \operatorname {asinh}^{2}{\left (c + d x \right )} + 24 b^{4} x - \frac {4 b^{4} \sqrt {c^{2} + 2 c d x + d^{2} x^{2} + 1} \operatorname {asinh}^{3}{\left (c + d x \right )}}{d} - \frac {24 b^{4} \sqrt {c^{2} + 2 c d x + d^{2} x^{2} + 1} \operatorname {asinh}{\left (c + d x \right )}}{d} & \text {for}\: d \neq 0 \\x \left (a + b \operatorname {asinh}{\left (c \right )}\right )^{4} & \text {otherwise} \end {cases} \]

[In]

integrate((a+b*asinh(d*x+c))**4,x)

[Out]

Piecewise((a**4*x + 4*a**3*b*c*asinh(c + d*x)/d + 4*a**3*b*x*asinh(c + d*x) - 4*a**3*b*sqrt(c**2 + 2*c*d*x + d
**2*x**2 + 1)/d + 6*a**2*b**2*c*asinh(c + d*x)**2/d + 6*a**2*b**2*x*asinh(c + d*x)**2 + 12*a**2*b**2*x - 12*a*
*2*b**2*sqrt(c**2 + 2*c*d*x + d**2*x**2 + 1)*asinh(c + d*x)/d + 4*a*b**3*c*asinh(c + d*x)**3/d + 24*a*b**3*c*a
sinh(c + d*x)/d + 4*a*b**3*x*asinh(c + d*x)**3 + 24*a*b**3*x*asinh(c + d*x) - 12*a*b**3*sqrt(c**2 + 2*c*d*x +
d**2*x**2 + 1)*asinh(c + d*x)**2/d - 24*a*b**3*sqrt(c**2 + 2*c*d*x + d**2*x**2 + 1)/d + b**4*c*asinh(c + d*x)*
*4/d + 12*b**4*c*asinh(c + d*x)**2/d + b**4*x*asinh(c + d*x)**4 + 12*b**4*x*asinh(c + d*x)**2 + 24*b**4*x - 4*
b**4*sqrt(c**2 + 2*c*d*x + d**2*x**2 + 1)*asinh(c + d*x)**3/d - 24*b**4*sqrt(c**2 + 2*c*d*x + d**2*x**2 + 1)*a
sinh(c + d*x)/d, Ne(d, 0)), (x*(a + b*asinh(c))**4, True))

Maxima [F]

\[ \int (a+b \text {arcsinh}(c+d x))^4 \, dx=\int { {\left (b \operatorname {arsinh}\left (d x + c\right ) + a\right )}^{4} \,d x } \]

[In]

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

[Out]

b^4*x*log(d*x + c + sqrt(d^2*x^2 + 2*c*d*x + c^2 + 1))^4 + a^4*x + 4*((d*x + c)*arcsinh(d*x + c) - sqrt((d*x +
c)^2 + 1))*a^3*b/d + integrate(2*(2*((c^3 + c)*a*b^3 + (a*b^3*d^3 - b^4*d^3)*x^3 + (3*a*b^3*c*d^2 - 2*b^4*c*d
^2)*x^2 + ((3*c^2*d + d)*a*b^3 - (c^2*d + d)*b^4)*x + ((c^2 + 1)*a*b^3 + (a*b^3*d^2 - b^4*d^2)*x^2 + (2*a*b^3*
c*d - b^4*c*d)*x)*sqrt(d^2*x^2 + 2*c*d*x + c^2 + 1))*log(d*x + c + sqrt(d^2*x^2 + 2*c*d*x + c^2 + 1))^3 + 3*(a
^2*b^2*d^3*x^3 + 3*a^2*b^2*c*d^2*x^2 + (3*c^2*d + d)*a^2*b^2*x + (c^3 + c)*a^2*b^2 + (a^2*b^2*d^2*x^2 + 2*a^2*
b^2*c*d*x + (c^2 + 1)*a^2*b^2)*sqrt(d^2*x^2 + 2*c*d*x + c^2 + 1))*log(d*x + c + sqrt(d^2*x^2 + 2*c*d*x + c^2 +
1))^2)/(d^3*x^3 + 3*c*d^2*x^2 + c^3 + (3*c^2*d + d)*x + (d^2*x^2 + 2*c*d*x + c^2 + 1)^(3/2) + c), x)

Giac [F]

\[ \int (a+b \text {arcsinh}(c+d x))^4 \, dx=\int { {\left (b \operatorname {arsinh}\left (d x + c\right ) + a\right )}^{4} \,d x } \]

[In]

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

[Out]

integrate((b*arcsinh(d*x + c) + a)^4, x)

Mupad [F(-1)]

Timed out. \[ \int (a+b \text {arcsinh}(c+d x))^4 \, dx=\int {\left (a+b\,\mathrm {asinh}\left (c+d\,x\right )\right )}^4 \,d x \]

[In]

int((a + b*asinh(c + d*x))^4,x)

[Out]

int((a + b*asinh(c + d*x))^4, x)