\(\int (c e+d e x) (a+b \text {arcsinh}(c+d x))^3 \, dx\) [140]

Optimal result

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

[Out]

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

Rubi [A] (verified)

Time = 0.15 (sec) , antiderivative size = 161, normalized size of antiderivative = 1.00, number of steps used = 8, number of rules used = 7, \(\frac {\text {number of rules}}{\text {integrand size}}\) = 0.333, Rules used = {5859, 12, 5776, 5812, 5783, 327, 221} \[ \int (c e+d e x) (a+b \text {arcsinh}(c+d x))^3 \, dx=\frac {3 b^2 e (c+d x)^2 (a+b \text {arcsinh}(c+d x))}{4 d}-\frac {3 b e (c+d x) \sqrt {(c+d x)^2+1} (a+b \text {arcsinh}(c+d x))^2}{4 d}+\frac {e (c+d x)^2 (a+b \text {arcsinh}(c+d x))^3}{2 d}+\frac {e (a+b \text {arcsinh}(c+d x))^3}{4 d}+\frac {3 b^3 e \text {arcsinh}(c+d x)}{8 d}-\frac {3 b^3 e (c+d x) \sqrt {(c+d x)^2+1}}{8 d} \]

[In]

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

[Out]

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

Rule 12

Int[(a_)*(u_), x_Symbol] :> Dist[a, Int[u, x], x] /; FreeQ[a, x] &&  !MatchQ[u, (b_)*(v_) /; FreeQ[b, x]]

Rule 221

Int[1/Sqrt[(a_) + (b_.)*(x_)^2], x_Symbol] :> Simp[ArcSinh[Rt[b, 2]*(x/Sqrt[a])]/Rt[b, 2], x] /; FreeQ[{a, b},
x] && GtQ[a, 0] && PosQ[b]

Rule 327

Int[((c_.)*(x_))^(m_)*((a_) + (b_.)*(x_)^(n_))^(p_), x_Symbol] :> Simp[c^(n - 1)*(c*x)^(m - n + 1)*((a + b*x^n
)^(p + 1)/(b*(m + n*p + 1))), x] - Dist[a*c^n*((m - n + 1)/(b*(m + n*p + 1))), Int[(c*x)^(m - n)*(a + b*x^n)^p
, x], x] /; FreeQ[{a, b, c, p}, x] && IGtQ[n, 0] && GtQ[m, n - 1] && NeQ[m + n*p + 1, 0] && IntBinomialQ[a, b,
c, n, m, p, x]

Rule 5776

Int[((a_.) + ArcSinh[(c_.)*(x_)]*(b_.))^(n_.)*((d_.)*(x_))^(m_.), x_Symbol] :> Simp[(d*x)^(m + 1)*((a + b*ArcS
inh[c*x])^n/(d*(m + 1))), x] - Dist[b*c*(n/(d*(m + 1))), Int[(d*x)^(m + 1)*((a + b*ArcSinh[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 5783

Int[((a_.) + ArcSinh[(c_.)*(x_)]*(b_.))^(n_.)/Sqrt[(d_) + (e_.)*(x_)^2], x_Symbol] :> Simp[(1/(b*c*(n + 1)))*S
imp[Sqrt[1 + c^2*x^2]/Sqrt[d + e*x^2]]*(a + b*ArcSinh[c*x])^(n + 1), x] /; FreeQ[{a, b, c, d, e, n}, x] && EqQ
[e, c^2*d] && NeQ[n, -1]

Rule 5812

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

Rule 5859

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

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

Mathematica [A] (verified)

Time = 0.27 (sec) , antiderivative size = 200, normalized size of antiderivative = 1.24 \[ \int (c e+d e x) (a+b \text {arcsinh}(c+d x))^3 \, dx=\frac {e \left (2 a \left (2 a^2+3 b^2\right ) (c+d x)^2-3 b \left (2 a^2+b^2\right ) (c+d x) \sqrt {1+(c+d x)^2}+3 b \left (2 a^2+b^2\right ) \text {arcsinh}(c+d x)-6 b (c+d x) \left (-2 a^2 (c+d x)-b^2 (c+d x)+2 a b \sqrt {1+(c+d x)^2}\right ) \text {arcsinh}(c+d x)+6 b^2 \left (a+2 a (c+d x)^2-b (c+d x) \sqrt {1+(c+d x)^2}\right ) \text {arcsinh}(c+d x)^2+2 b^3 \left (1+2 (c+d x)^2\right ) \text {arcsinh}(c+d x)^3\right )}{8 d} \]

[In]

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

[Out]

(e*(2*a*(2*a^2 + 3*b^2)*(c + d*x)^2 - 3*b*(2*a^2 + b^2)*(c + d*x)*Sqrt[1 + (c + d*x)^2] + 3*b*(2*a^2 + b^2)*Ar
cSinh[c + d*x] - 6*b*(c + d*x)*(-2*a^2*(c + d*x) - b^2*(c + d*x) + 2*a*b*Sqrt[1 + (c + d*x)^2])*ArcSinh[c + d*
x] + 6*b^2*(a + 2*a*(c + d*x)^2 - b*(c + d*x)*Sqrt[1 + (c + d*x)^2])*ArcSinh[c + d*x]^2 + 2*b^3*(1 + 2*(c + d*
x)^2)*ArcSinh[c + d*x]^3))/(8*d)

Maple [A] (verified)

Time = 0.08 (sec) , antiderivative size = 243, normalized size of antiderivative = 1.51

method	result	size
derivativedivides	\(\frac {\frac {e \,a^{3} \left (d x +c \right )^{2}}{2}+e \,b^{3} \left (\frac {\operatorname {arcsinh}\left (d x +c \right )^{3} \left (1+\left (d x +c \right )^{2}\right )}{2}-\frac {3 \operatorname {arcsinh}\left (d x +c \right )^{2} \sqrt {1+\left (d x +c \right )^{2}}\, \left (d x +c \right )}{4}-\frac {\operatorname {arcsinh}\left (d x +c \right )^{3}}{4}+\frac {3 \left (1+\left (d x +c \right )^{2}\right ) \operatorname {arcsinh}\left (d x +c \right )}{4}-\frac {3 \left (d x +c \right ) \sqrt {1+\left (d x +c \right )^{2}}}{8}-\frac {3 \,\operatorname {arcsinh}\left (d x +c \right )}{8}\right )+3 e a \,b^{2} \left (\frac {\operatorname {arcsinh}\left (d x +c \right )^{2} \left (1+\left (d x +c \right )^{2}\right )}{2}-\frac {\sqrt {1+\left (d x +c \right )^{2}}\, \left (d x +c \right ) \operatorname {arcsinh}\left (d x +c \right )}{2}-\frac {\operatorname {arcsinh}\left (d x +c \right )^{2}}{4}+\frac {\left (d x +c \right )^{2}}{4}+\frac {1}{4}\right )+3 e \,a^{2} b \left (\frac {\left (d x +c \right )^{2} \operatorname {arcsinh}\left (d x +c \right )}{2}-\frac {\left (d x +c \right ) \sqrt {1+\left (d x +c \right )^{2}}}{4}+\frac {\operatorname {arcsinh}\left (d x +c \right )}{4}\right )}{d}\)	\(243\)
default	\(\frac {\frac {e \,a^{3} \left (d x +c \right )^{2}}{2}+e \,b^{3} \left (\frac {\operatorname {arcsinh}\left (d x +c \right )^{3} \left (1+\left (d x +c \right )^{2}\right )}{2}-\frac {3 \operatorname {arcsinh}\left (d x +c \right )^{2} \sqrt {1+\left (d x +c \right )^{2}}\, \left (d x +c \right )}{4}-\frac {\operatorname {arcsinh}\left (d x +c \right )^{3}}{4}+\frac {3 \left (1+\left (d x +c \right )^{2}\right ) \operatorname {arcsinh}\left (d x +c \right )}{4}-\frac {3 \left (d x +c \right ) \sqrt {1+\left (d x +c \right )^{2}}}{8}-\frac {3 \,\operatorname {arcsinh}\left (d x +c \right )}{8}\right )+3 e a \,b^{2} \left (\frac {\operatorname {arcsinh}\left (d x +c \right )^{2} \left (1+\left (d x +c \right )^{2}\right )}{2}-\frac {\sqrt {1+\left (d x +c \right )^{2}}\, \left (d x +c \right ) \operatorname {arcsinh}\left (d x +c \right )}{2}-\frac {\operatorname {arcsinh}\left (d x +c \right )^{2}}{4}+\frac {\left (d x +c \right )^{2}}{4}+\frac {1}{4}\right )+3 e \,a^{2} b \left (\frac {\left (d x +c \right )^{2} \operatorname {arcsinh}\left (d x +c \right )}{2}-\frac {\left (d x +c \right ) \sqrt {1+\left (d x +c \right )^{2}}}{4}+\frac {\operatorname {arcsinh}\left (d x +c \right )}{4}\right )}{d}\)	\(243\)
parts	\(e \,a^{3} \left (\frac {1}{2} d \,x^{2}+c x \right )+\frac {e \,b^{3} \left (\frac {\operatorname {arcsinh}\left (d x +c \right )^{3} \left (1+\left (d x +c \right )^{2}\right )}{2}-\frac {3 \operatorname {arcsinh}\left (d x +c \right )^{2} \sqrt {1+\left (d x +c \right )^{2}}\, \left (d x +c \right )}{4}-\frac {\operatorname {arcsinh}\left (d x +c \right )^{3}}{4}+\frac {3 \left (1+\left (d x +c \right )^{2}\right ) \operatorname {arcsinh}\left (d x +c \right )}{4}-\frac {3 \left (d x +c \right ) \sqrt {1+\left (d x +c \right )^{2}}}{8}-\frac {3 \,\operatorname {arcsinh}\left (d x +c \right )}{8}\right )}{d}+\frac {3 e a \,b^{2} \left (\frac {\operatorname {arcsinh}\left (d x +c \right )^{2} \left (1+\left (d x +c \right )^{2}\right )}{2}-\frac {\sqrt {1+\left (d x +c \right )^{2}}\, \left (d x +c \right ) \operatorname {arcsinh}\left (d x +c \right )}{2}-\frac {\operatorname {arcsinh}\left (d x +c \right )^{2}}{4}+\frac {\left (d x +c \right )^{2}}{4}+\frac {1}{4}\right )}{d}+\frac {3 e \,a^{2} b \left (\frac {\left (d x +c \right )^{2} \operatorname {arcsinh}\left (d x +c \right )}{2}-\frac {\left (d x +c \right ) \sqrt {1+\left (d x +c \right )^{2}}}{4}+\frac {\operatorname {arcsinh}\left (d x +c \right )}{4}\right )}{d}\)	\(250\)

[In]

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

[Out]

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

Fricas [B] (verification not implemented)

Leaf count of result is larger than twice the leaf count of optimal. 391 vs. \(2 (145) = 290\).

Time = 0.31 (sec) , antiderivative size = 391, normalized size of antiderivative = 2.43 \[ \int (c e+d e x) (a+b \text {arcsinh}(c+d x))^3 \, dx=\frac {2 \, {\left (2 \, a^{3} + 3 \, a b^{2}\right )} d^{2} e x^{2} + 4 \, {\left (2 \, a^{3} + 3 \, a b^{2}\right )} c d e x + 2 \, {\left (2 \, b^{3} d^{2} e x^{2} + 4 \, b^{3} c d e x + {\left (2 \, b^{3} c^{2} + b^{3}\right )} e\right )} \log \left (d x + c + \sqrt {d^{2} x^{2} + 2 \, c d x + c^{2} + 1}\right )^{3} + 6 \, {\left (2 \, a b^{2} d^{2} e x^{2} + 4 \, a b^{2} c d e x + {\left (2 \, a b^{2} c^{2} + a b^{2}\right )} e - {\left (b^{3} d e x + b^{3} c e\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 )^{2} + 3 \, {\left (2 \, {\left (2 \, a^{2} b + b^{3}\right )} d^{2} e x^{2} + 4 \, {\left (2 \, a^{2} b + b^{3}\right )} c d e x + {\left (2 \, a^{2} b + b^{3} + 2 \, {\left (2 \, a^{2} b + b^{3}\right )} c^{2}\right )} e - 4 \, {\left (a b^{2} d e x + a b^{2} c e\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 ) - 3 \, {\left ({\left (2 \, a^{2} b + b^{3}\right )} d e x + {\left (2 \, a^{2} b + b^{3}\right )} c e\right )} \sqrt {d^{2} x^{2} + 2 \, c d x + c^{2} + 1}}{8 \, d} \]

[In]

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

[Out]

1/8*(2*(2*a^3 + 3*a*b^2)*d^2*e*x^2 + 4*(2*a^3 + 3*a*b^2)*c*d*e*x + 2*(2*b^3*d^2*e*x^2 + 4*b^3*c*d*e*x + (2*b^3
*c^2 + b^3)*e)*log(d*x + c + sqrt(d^2*x^2 + 2*c*d*x + c^2 + 1))^3 + 6*(2*a*b^2*d^2*e*x^2 + 4*a*b^2*c*d*e*x + (
2*a*b^2*c^2 + a*b^2)*e - (b^3*d*e*x + b^3*c*e)*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 + 3*(2*(2*a^2*b + b^3)*d^2*e*x^2 + 4*(2*a^2*b + b^3)*c*d*e*x + (2*a^2*b + b^3 + 2*(2*a^
2*b + b^3)*c^2)*e - 4*(a*b^2*d*e*x + a*b^2*c*e)*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*((2*a^2*b + b^3)*d*e*x + (2*a^2*b + b^3)*c*e)*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. 685 vs. \(2 (148) = 296\).

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

[In]

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

[Out]

Piecewise((a**3*c*e*x + a**3*d*e*x**2/2 + 3*a**2*b*c**2*e*asinh(c + d*x)/(2*d) + 3*a**2*b*c*e*x*asinh(c + d*x)
- 3*a**2*b*c*e*sqrt(c**2 + 2*c*d*x + d**2*x**2 + 1)/(4*d) + 3*a**2*b*d*e*x**2*asinh(c + d*x)/2 - 3*a**2*b*e*x
*sqrt(c**2 + 2*c*d*x + d**2*x**2 + 1)/4 + 3*a**2*b*e*asinh(c + d*x)/(4*d) + 3*a*b**2*c**2*e*asinh(c + d*x)**2/
(2*d) + 3*a*b**2*c*e*x*asinh(c + d*x)**2 + 3*a*b**2*c*e*x/2 - 3*a*b**2*c*e*sqrt(c**2 + 2*c*d*x + d**2*x**2 + 1
)*asinh(c + d*x)/(2*d) + 3*a*b**2*d*e*x**2*asinh(c + d*x)**2/2 + 3*a*b**2*d*e*x**2/4 - 3*a*b**2*e*x*sqrt(c**2
+ 2*c*d*x + d**2*x**2 + 1)*asinh(c + d*x)/2 + 3*a*b**2*e*asinh(c + d*x)**2/(4*d) + b**3*c**2*e*asinh(c + d*x)*
*3/(2*d) + 3*b**3*c**2*e*asinh(c + d*x)/(4*d) + b**3*c*e*x*asinh(c + d*x)**3 + 3*b**3*c*e*x*asinh(c + d*x)/2 -
3*b**3*c*e*sqrt(c**2 + 2*c*d*x + d**2*x**2 + 1)*asinh(c + d*x)**2/(4*d) - 3*b**3*c*e*sqrt(c**2 + 2*c*d*x + d*
*2*x**2 + 1)/(8*d) + b**3*d*e*x**2*asinh(c + d*x)**3/2 + 3*b**3*d*e*x**2*asinh(c + d*x)/4 - 3*b**3*e*x*sqrt(c*
*2 + 2*c*d*x + d**2*x**2 + 1)*asinh(c + d*x)**2/4 - 3*b**3*e*x*sqrt(c**2 + 2*c*d*x + d**2*x**2 + 1)/8 + b**3*e
*asinh(c + d*x)**3/(4*d) + 3*b**3*e*asinh(c + d*x)/(8*d), Ne(d, 0)), (c*e*x*(a + b*asinh(c))**3, True))

Maxima [F]

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

[In]

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

[Out]

1/2*a^3*d*e*x^2 + 3/4*(2*x^2*arcsinh(d*x + c) - d*(3*c^2*arcsinh(2*(d^2*x + c*d)/sqrt(-4*c^2*d^2 + 4*(c^2 + 1)
*d^2))/d^3 + sqrt(d^2*x^2 + 2*c*d*x + c^2 + 1)*x/d^2 - (c^2 + 1)*arcsinh(2*(d^2*x + c*d)/sqrt(-4*c^2*d^2 + 4*(
c^2 + 1)*d^2))/d^3 - 3*sqrt(d^2*x^2 + 2*c*d*x + c^2 + 1)*c/d^3))*a^2*b*d*e + a^3*c*e*x + 3*((d*x + c)*arcsinh(
d*x + c) - sqrt((d*x + c)^2 + 1))*a^2*b*c*e/d + 1/2*(b^3*d*e*x^2 + 2*b^3*c*e*x)*log(d*x + c + sqrt(d^2*x^2 + 2
*c*d*x + c^2 + 1))^3 + integrate(3/2*((2*a*b^2*d^4*e - b^3*d^4*e)*x^4 + 2*(c^4*e + c^2*e)*a*b^2 + 4*(2*a*b^2*c
*d^3*e - b^3*c*d^3*e)*x^3 + (2*(6*c^2*d^2*e + d^2*e)*a*b^2 - (5*c^2*d^2*e + d^2*e)*b^3)*x^2 + 2*(2*(2*c^3*d*e
+ c*d*e)*a*b^2 - (c^3*d*e + c*d*e)*b^3)*x + (2*(c^3*e + c*e)*a*b^2 + (2*a*b^2*d^3*e - b^3*d^3*e)*x^3 + 3*(2*a*
b^2*c*d^2*e - b^3*c*d^2*e)*x^2 - 2*(b^3*c^2*d*e - (3*c^2*d*e + d*e)*a*b^2)*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))^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 (c e+d e x) (a+b \text {arcsinh}(c+d x))^3 \, dx=\int { {\left (d e x + c e\right )} {\left (b \operatorname {arsinh}\left (d x + c\right ) + a\right )}^{3} \,d x } \]

[In]

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

[Out]

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

Mupad [F(-1)]

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

[In]

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

[Out]

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