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

Optimal result
Mathematica [A] (verified)
Rubi [A] (verified)
Maple [A] (verified)
Fricas [A] (verification not implemented)
Sympy [B] (verification not implemented)
Maxima [A] (verification not implemented)
Giac [F(-2)]
Mupad [F(-1)]
Reduce [F]

Optimal result

Integrand size = 18, antiderivative size = 368 \[ \int (d+e x)^3 (a+b \text {arcsinh}(c x))^2 \, dx=2 b^2 d^3 x-\frac {4 b^2 d e^2 x}{3 c^2}+\frac {3}{4} b^2 d^2 e x^2-\frac {3 b^2 e^3 x^2}{32 c^2}+\frac {2}{9} b^2 d e^2 x^3+\frac {1}{32} b^2 e^3 x^4-\frac {2 b d^3 \sqrt {1+c^2 x^2} (a+b \text {arcsinh}(c x))}{c}+\frac {4 b d e^2 \sqrt {1+c^2 x^2} (a+b \text {arcsinh}(c x))}{3 c^3}-\frac {3 b d^2 e x \sqrt {1+c^2 x^2} (a+b \text {arcsinh}(c x))}{2 c}+\frac {3 b e^3 x \sqrt {1+c^2 x^2} (a+b \text {arcsinh}(c x))}{16 c^3}-\frac {2 b d e^2 x^2 \sqrt {1+c^2 x^2} (a+b \text {arcsinh}(c x))}{3 c}-\frac {b e^3 x^3 \sqrt {1+c^2 x^2} (a+b \text {arcsinh}(c x))}{8 c}-\frac {d^4 (a+b \text {arcsinh}(c x))^2}{4 e}+\frac {3 d^2 e (a+b \text {arcsinh}(c x))^2}{4 c^2}-\frac {3 e^3 (a+b \text {arcsinh}(c x))^2}{32 c^4}+\frac {(d+e x)^4 (a+b \text {arcsinh}(c x))^2}{4 e} \] Output:

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

Mathematica [A] (verified)

Time = 0.33 (sec) , antiderivative size = 354, normalized size of antiderivative = 0.96 \[ \int (d+e x)^3 (a+b \text {arcsinh}(c x))^2 \, dx=\frac {c \left (72 a^2 c^3 x \left (4 d^3+6 d^2 e x+4 d e^2 x^2+e^3 x^3\right )-6 a b \sqrt {1+c^2 x^2} \left (-e^2 (64 d+9 e x)+c^2 \left (96 d^3+72 d^2 e x+32 d e^2 x^2+6 e^3 x^3\right )\right )+b^2 c x \left (-3 e^2 (128 d+9 e x)+c^2 \left (576 d^3+216 d^2 e x+64 d e^2 x^2+9 e^3 x^3\right )\right )\right )-6 b \left (-3 a \left (24 c^2 d^2 e-3 e^3+8 c^4 x \left (4 d^3+6 d^2 e x+4 d e^2 x^2+e^3 x^3\right )\right )+b c \sqrt {1+c^2 x^2} \left (-e^2 (64 d+9 e x)+c^2 \left (96 d^3+72 d^2 e x+32 d e^2 x^2+6 e^3 x^3\right )\right )\right ) \text {arcsinh}(c x)+9 b^2 \left (24 c^2 d^2 e-3 e^3+8 c^4 x \left (4 d^3+6 d^2 e x+4 d e^2 x^2+e^3 x^3\right )\right ) \text {arcsinh}(c x)^2}{288 c^4} \] Input:

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

Output:

(c*(72*a^2*c^3*x*(4*d^3 + 6*d^2*e*x + 4*d*e^2*x^2 + e^3*x^3) - 6*a*b*Sqrt[ 
1 + c^2*x^2]*(-(e^2*(64*d + 9*e*x)) + c^2*(96*d^3 + 72*d^2*e*x + 32*d*e^2* 
x^2 + 6*e^3*x^3)) + b^2*c*x*(-3*e^2*(128*d + 9*e*x) + c^2*(576*d^3 + 216*d 
^2*e*x + 64*d*e^2*x^2 + 9*e^3*x^3))) - 6*b*(-3*a*(24*c^2*d^2*e - 3*e^3 + 8 
*c^4*x*(4*d^3 + 6*d^2*e*x + 4*d*e^2*x^2 + e^3*x^3)) + b*c*Sqrt[1 + c^2*x^2 
]*(-(e^2*(64*d + 9*e*x)) + c^2*(96*d^3 + 72*d^2*e*x + 32*d*e^2*x^2 + 6*e^3 
*x^3)))*ArcSinh[c*x] + 9*b^2*(24*c^2*d^2*e - 3*e^3 + 8*c^4*x*(4*d^3 + 6*d^ 
2*e*x + 4*d*e^2*x^2 + e^3*x^3))*ArcSinh[c*x]^2)/(288*c^4)
 

Rubi [A] (verified)

Time = 1.06 (sec) , antiderivative size = 387, normalized size of antiderivative = 1.05, number of steps used = 3, number of rules used = 3, \(\frac {\text {number of rules}}{\text {integrand size}}\) = 0.167, Rules used = {6243, 6253, 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 (d+e x)^3 (a+b \text {arcsinh}(c x))^2 \, dx\)

\(\Big \downarrow \) 6243

\(\displaystyle \frac {(d+e x)^4 (a+b \text {arcsinh}(c x))^2}{4 e}-\frac {b c \int \frac {(d+e x)^4 (a+b \text {arcsinh}(c x))}{\sqrt {c^2 x^2+1}}dx}{2 e}\)

\(\Big \downarrow \) 6253

\(\displaystyle \frac {(d+e x)^4 (a+b \text {arcsinh}(c x))^2}{4 e}-\frac {b c \int \left (\frac {(a+b \text {arcsinh}(c x)) d^4}{\sqrt {c^2 x^2+1}}+\frac {4 e x (a+b \text {arcsinh}(c x)) d^3}{\sqrt {c^2 x^2+1}}+\frac {6 e^2 x^2 (a+b \text {arcsinh}(c x)) d^2}{\sqrt {c^2 x^2+1}}+\frac {4 e^3 x^3 (a+b \text {arcsinh}(c x)) d}{\sqrt {c^2 x^2+1}}+\frac {e^4 x^4 (a+b \text {arcsinh}(c x))}{\sqrt {c^2 x^2+1}}\right )dx}{2 e}\)

\(\Big \downarrow \) 2009

\(\displaystyle \frac {(d+e x)^4 (a+b \text {arcsinh}(c x))^2}{4 e}-\frac {b c \left (\frac {3 e^4 (a+b \text {arcsinh}(c x))^2}{16 b c^5}-\frac {3 d^2 e^2 (a+b \text {arcsinh}(c x))^2}{2 b c^3}+\frac {4 d^3 e \sqrt {c^2 x^2+1} (a+b \text {arcsinh}(c x))}{c^2}+\frac {3 d^2 e^2 x \sqrt {c^2 x^2+1} (a+b \text {arcsinh}(c x))}{c^2}+\frac {4 d e^3 x^2 \sqrt {c^2 x^2+1} (a+b \text {arcsinh}(c x))}{3 c^2}+\frac {e^4 x^3 \sqrt {c^2 x^2+1} (a+b \text {arcsinh}(c x))}{4 c^2}-\frac {8 d e^3 \sqrt {c^2 x^2+1} (a+b \text {arcsinh}(c x))}{3 c^4}-\frac {3 e^4 x \sqrt {c^2 x^2+1} (a+b \text {arcsinh}(c x))}{8 c^4}+\frac {d^4 (a+b \text {arcsinh}(c x))^2}{2 b c}+\frac {8 b d e^3 x}{3 c^3}+\frac {3 b e^4 x^2}{16 c^3}-\frac {4 b d^3 e x}{c}-\frac {3 b d^2 e^2 x^2}{2 c}-\frac {4 b d e^3 x^3}{9 c}-\frac {b e^4 x^4}{16 c}\right )}{2 e}\)

Input:

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

Output:

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

Defintions of rubi rules used

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

rule 6243
Int[((a_.) + ArcSinh[(c_.)*(x_)]*(b_.))^(n_.)*((d_.) + (e_.)*(x_))^(m_.), x 
_Symbol] :> Simp[(d + e*x)^(m + 1)*((a + b*ArcSinh[c*x])^n/(e*(m + 1))), x] 
 - Simp[b*c*(n/(e*(m + 1)))   Int[(d + e*x)^(m + 1)*((a + b*ArcSinh[c*x])^( 
n - 1)/Sqrt[1 + c^2*x^2]), x], x] /; FreeQ[{a, b, c, d, e, m}, x] && IGtQ[n 
, 0] && NeQ[m, -1]
 

rule 6253
Int[((a_.) + ArcSinh[(c_.)*(x_)]*(b_.))^(n_.)*((f_) + (g_.)*(x_))^(m_.)*((d 
_) + (e_.)*(x_)^2)^(p_), x_Symbol] :> Int[ExpandIntegrand[(d + e*x^2)^p*(a 
+ b*ArcSinh[c*x])^n, (f + g*x)^m, x], x] /; FreeQ[{a, b, c, d, e, f, g}, x] 
 && EqQ[e, c^2*d] && IGtQ[m, 0] && IntegerQ[p + 1/2] && GtQ[d, 0] && IGtQ[n 
, 0] && ((EqQ[n, 1] && GtQ[p, -1]) || GtQ[p, 0] || EqQ[m, 1] || (EqQ[m, 2] 
&& LtQ[p, -2]))
 
Maple [A] (verified)

Time = 2.86 (sec) , antiderivative size = 526, normalized size of antiderivative = 1.43

method result size
derivativedivides \(\frac {\frac {a^{2} \left (c e x +c d \right )^{4}}{4 c^{3} e}+\frac {b^{2} \left (d^{3} c^{3} \left (\operatorname {arcsinh}\left (x c \right )^{2} x c -2 \,\operatorname {arcsinh}\left (x c \right ) \sqrt {c^{2} x^{2}+1}+2 x c \right )+\frac {3 d^{2} c^{2} e \left (2 \operatorname {arcsinh}\left (x c \right )^{2} x^{2} c^{2}-2 \,\operatorname {arcsinh}\left (x c \right ) \sqrt {c^{2} x^{2}+1}\, x c +\operatorname {arcsinh}\left (x c \right )^{2}+c^{2} x^{2}+1\right )}{4}+\frac {d c \,e^{2} \left (9 \operatorname {arcsinh}\left (x c \right )^{2} x^{3} c^{3}-6 \sqrt {c^{2} x^{2}+1}\, \operatorname {arcsinh}\left (x c \right ) x^{2} c^{2}+2 x^{3} c^{3}+12 \,\operatorname {arcsinh}\left (x c \right ) \sqrt {c^{2} x^{2}+1}-12 x c \right )}{9}+\frac {e^{3} \left (8 \operatorname {arcsinh}\left (x c \right )^{2} x^{4} c^{4}-4 \,\operatorname {arcsinh}\left (x c \right ) \sqrt {c^{2} x^{2}+1}\, x^{3} c^{3}+c^{4} x^{4}+6 \,\operatorname {arcsinh}\left (x c \right ) \sqrt {c^{2} x^{2}+1}\, x c -3 \operatorname {arcsinh}\left (x c \right )^{2}-3 c^{2} x^{2}-3\right )}{32}\right )}{c^{3}}+\frac {2 a b \left (\frac {\operatorname {arcsinh}\left (x c \right ) c^{4} d^{4}}{4 e}+\operatorname {arcsinh}\left (x c \right ) c^{4} d^{3} x +\frac {3 e \,\operatorname {arcsinh}\left (x c \right ) c^{4} d^{2} x^{2}}{2}+e^{2} \operatorname {arcsinh}\left (x c \right ) c^{4} d \,x^{3}+\frac {e^{3} \operatorname {arcsinh}\left (x c \right ) x^{4} c^{4}}{4}-\frac {c^{4} d^{4} \operatorname {arcsinh}\left (x c \right )+e^{4} \left (\frac {\sqrt {c^{2} x^{2}+1}\, c^{3} x^{3}}{4}-\frac {3 \sqrt {c^{2} x^{2}+1}\, x c}{8}+\frac {3 \,\operatorname {arcsinh}\left (x c \right )}{8}\right )+4 d c \,e^{3} \left (\frac {x^{2} c^{2} \sqrt {c^{2} x^{2}+1}}{3}-\frac {2 \sqrt {c^{2} x^{2}+1}}{3}\right )+6 d^{2} c^{2} e^{2} \left (\frac {\sqrt {c^{2} x^{2}+1}\, x c}{2}-\frac {\operatorname {arcsinh}\left (x c \right )}{2}\right )+4 d^{3} c^{3} e \sqrt {c^{2} x^{2}+1}}{4 e}\right )}{c^{3}}}{c}\) \(526\)
default \(\frac {\frac {a^{2} \left (c e x +c d \right )^{4}}{4 c^{3} e}+\frac {b^{2} \left (d^{3} c^{3} \left (\operatorname {arcsinh}\left (x c \right )^{2} x c -2 \,\operatorname {arcsinh}\left (x c \right ) \sqrt {c^{2} x^{2}+1}+2 x c \right )+\frac {3 d^{2} c^{2} e \left (2 \operatorname {arcsinh}\left (x c \right )^{2} x^{2} c^{2}-2 \,\operatorname {arcsinh}\left (x c \right ) \sqrt {c^{2} x^{2}+1}\, x c +\operatorname {arcsinh}\left (x c \right )^{2}+c^{2} x^{2}+1\right )}{4}+\frac {d c \,e^{2} \left (9 \operatorname {arcsinh}\left (x c \right )^{2} x^{3} c^{3}-6 \sqrt {c^{2} x^{2}+1}\, \operatorname {arcsinh}\left (x c \right ) x^{2} c^{2}+2 x^{3} c^{3}+12 \,\operatorname {arcsinh}\left (x c \right ) \sqrt {c^{2} x^{2}+1}-12 x c \right )}{9}+\frac {e^{3} \left (8 \operatorname {arcsinh}\left (x c \right )^{2} x^{4} c^{4}-4 \,\operatorname {arcsinh}\left (x c \right ) \sqrt {c^{2} x^{2}+1}\, x^{3} c^{3}+c^{4} x^{4}+6 \,\operatorname {arcsinh}\left (x c \right ) \sqrt {c^{2} x^{2}+1}\, x c -3 \operatorname {arcsinh}\left (x c \right )^{2}-3 c^{2} x^{2}-3\right )}{32}\right )}{c^{3}}+\frac {2 a b \left (\frac {\operatorname {arcsinh}\left (x c \right ) c^{4} d^{4}}{4 e}+\operatorname {arcsinh}\left (x c \right ) c^{4} d^{3} x +\frac {3 e \,\operatorname {arcsinh}\left (x c \right ) c^{4} d^{2} x^{2}}{2}+e^{2} \operatorname {arcsinh}\left (x c \right ) c^{4} d \,x^{3}+\frac {e^{3} \operatorname {arcsinh}\left (x c \right ) x^{4} c^{4}}{4}-\frac {c^{4} d^{4} \operatorname {arcsinh}\left (x c \right )+e^{4} \left (\frac {\sqrt {c^{2} x^{2}+1}\, c^{3} x^{3}}{4}-\frac {3 \sqrt {c^{2} x^{2}+1}\, x c}{8}+\frac {3 \,\operatorname {arcsinh}\left (x c \right )}{8}\right )+4 d c \,e^{3} \left (\frac {x^{2} c^{2} \sqrt {c^{2} x^{2}+1}}{3}-\frac {2 \sqrt {c^{2} x^{2}+1}}{3}\right )+6 d^{2} c^{2} e^{2} \left (\frac {\sqrt {c^{2} x^{2}+1}\, x c}{2}-\frac {\operatorname {arcsinh}\left (x c \right )}{2}\right )+4 d^{3} c^{3} e \sqrt {c^{2} x^{2}+1}}{4 e}\right )}{c^{3}}}{c}\) \(526\)
orering \(\frac {\left (111 e^{5} c^{4} x^{6}+699 e^{4} c^{4} x^{5} d +1928 e^{3} c^{4} x^{4} d^{2}+3480 e^{2} c^{4} x^{3} d^{3}+672 c^{4} d^{4} e \,x^{2}-63 e^{5} c^{2} x^{4}+192 c^{4} d^{5} x -1079 e^{4} c^{2} x^{3} d +1632 e^{3} c^{2} x^{2} d^{2}+3600 e^{2} c^{2} d^{3} x +720 c^{2} d^{4} e -180 e^{5} x^{2}-2010 e^{4} d x -402 e^{3} d^{2}\right ) \left (a +b \,\operatorname {arcsinh}\left (x c \right )\right )^{2}}{192 c^{4} \left (e x +d \right )^{2}}-\frac {\left (81 e^{4} c^{4} x^{6}+539 e^{3} c^{4} x^{5} d +1640 e^{2} c^{4} x^{4} d^{2}+3672 e \,c^{4} x^{3} d^{3}-99 e^{4} c^{2} x^{4}-1719 e^{3} c^{2} x^{3} d +1920 e^{2} c^{2} x^{2} d^{2}+4464 e \,c^{2} d^{3} x +576 c^{2} d^{4}-216 e^{4} x^{2}-2742 e^{3} d x -384 d^{2} e^{2}\right ) \left (3 \left (e x +d \right )^{2} \left (a +b \,\operatorname {arcsinh}\left (x c \right )\right )^{2} e +\frac {2 \left (e x +d \right )^{3} \left (a +b \,\operatorname {arcsinh}\left (x c \right )\right ) b c}{\sqrt {c^{2} x^{2}+1}}\right )}{576 c^{4} \left (e x +d \right )^{4}}+\frac {x \left (9 e^{3} c^{2} x^{3}+64 e^{2} c^{2} x^{2} d +216 e \,c^{2} d^{2} x +576 c^{2} d^{3}-27 e^{3} x -384 e^{2} d \right ) \left (c^{2} x^{2}+1\right ) \left (6 \left (e x +d \right ) \left (a +b \,\operatorname {arcsinh}\left (x c \right )\right )^{2} e^{2}+\frac {12 \left (e x +d \right )^{2} \left (a +b \,\operatorname {arcsinh}\left (x c \right )\right ) e b c}{\sqrt {c^{2} x^{2}+1}}+\frac {2 \left (e x +d \right )^{3} b^{2} c^{2}}{c^{2} x^{2}+1}-\frac {2 \left (e x +d \right )^{3} \left (a +b \,\operatorname {arcsinh}\left (x c \right )\right ) b \,c^{3} x}{\left (c^{2} x^{2}+1\right )^{\frac {3}{2}}}\right )}{576 c^{4} \left (e x +d \right )^{3}}\) \(556\)
parts \(\frac {a^{2} \left (e x +d \right )^{4}}{4 e}+\frac {b^{2} \left (72 \operatorname {arcsinh}\left (x c \right )^{2} x^{4} c^{4} e^{3}+288 \operatorname {arcsinh}\left (x c \right )^{2} x^{3} c^{4} d \,e^{2}+432 \operatorname {arcsinh}\left (x c \right )^{2} x^{2} c^{4} d^{2} e +288 \operatorname {arcsinh}\left (x c \right )^{2} x \,c^{4} d^{3}-36 \,\operatorname {arcsinh}\left (x c \right ) \sqrt {c^{2} x^{2}+1}\, x^{3} c^{3} e^{3}-192 \,\operatorname {arcsinh}\left (x c \right ) \sqrt {c^{2} x^{2}+1}\, x^{2} c^{3} d \,e^{2}-432 \,\operatorname {arcsinh}\left (x c \right ) \sqrt {c^{2} x^{2}+1}\, x \,c^{3} d^{2} e -576 \,\operatorname {arcsinh}\left (x c \right ) \sqrt {c^{2} x^{2}+1}\, c^{3} d^{3}+216 \operatorname {arcsinh}\left (x c \right )^{2} c^{2} d^{2} e +9 c^{4} e^{3} x^{4}+64 c^{4} d \,e^{2} x^{3}+216 c^{4} d^{2} e \,x^{2}+576 c^{4} d^{3} x +54 \,\operatorname {arcsinh}\left (x c \right ) \sqrt {c^{2} x^{2}+1}\, x c \,e^{3}+384 \,\operatorname {arcsinh}\left (x c \right ) \sqrt {c^{2} x^{2}+1}\, c d \,e^{2}-27 \operatorname {arcsinh}\left (x c \right )^{2} e^{3}-27 x^{2} c^{2} e^{3}-384 c^{2} d \,e^{2} x +216 c^{2} d^{2} e -27 e^{3}\right )}{288 c^{4}}+\frac {2 a b \left (\frac {c \,e^{3} \operatorname {arcsinh}\left (x c \right ) x^{4}}{4}+c \,e^{2} \operatorname {arcsinh}\left (x c \right ) x^{3} d +\frac {3 c e \,\operatorname {arcsinh}\left (x c \right ) x^{2} d^{2}}{2}+\operatorname {arcsinh}\left (x c \right ) x c \,d^{3}+\frac {c \,\operatorname {arcsinh}\left (x c \right ) d^{4}}{4 e}-\frac {c^{4} d^{4} \operatorname {arcsinh}\left (x c \right )+e^{4} \left (\frac {\sqrt {c^{2} x^{2}+1}\, c^{3} x^{3}}{4}-\frac {3 \sqrt {c^{2} x^{2}+1}\, x c}{8}+\frac {3 \,\operatorname {arcsinh}\left (x c \right )}{8}\right )+4 d c \,e^{3} \left (\frac {x^{2} c^{2} \sqrt {c^{2} x^{2}+1}}{3}-\frac {2 \sqrt {c^{2} x^{2}+1}}{3}\right )+6 d^{2} c^{2} e^{2} \left (\frac {\sqrt {c^{2} x^{2}+1}\, x c}{2}-\frac {\operatorname {arcsinh}\left (x c \right )}{2}\right )+4 d^{3} c^{3} e \sqrt {c^{2} x^{2}+1}}{4 c^{3} e}\right )}{c}\) \(573\)

Input:

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

Output:

1/c*(1/4*a^2/c^3*(c*e*x+c*d)^4/e+b^2/c^3*(d^3*c^3*(arcsinh(x*c)^2*x*c-2*ar 
csinh(x*c)*(c^2*x^2+1)^(1/2)+2*x*c)+3/4*d^2*c^2*e*(2*arcsinh(x*c)^2*x^2*c^ 
2-2*arcsinh(x*c)*(c^2*x^2+1)^(1/2)*x*c+arcsinh(x*c)^2+c^2*x^2+1)+1/9*d*c*e 
^2*(9*arcsinh(x*c)^2*x^3*c^3-6*(c^2*x^2+1)^(1/2)*arcsinh(x*c)*x^2*c^2+2*x^ 
3*c^3+12*arcsinh(x*c)*(c^2*x^2+1)^(1/2)-12*x*c)+1/32*e^3*(8*arcsinh(x*c)^2 
*x^4*c^4-4*arcsinh(x*c)*(c^2*x^2+1)^(1/2)*x^3*c^3+c^4*x^4+6*arcsinh(x*c)*( 
c^2*x^2+1)^(1/2)*x*c-3*arcsinh(x*c)^2-3*c^2*x^2-3))+2*a*b/c^3*(1/4/e*arcsi 
nh(x*c)*c^4*d^4+arcsinh(x*c)*c^4*d^3*x+3/2*e*arcsinh(x*c)*c^4*d^2*x^2+e^2* 
arcsinh(x*c)*c^4*d*x^3+1/4*e^3*arcsinh(x*c)*x^4*c^4-1/4/e*(c^4*d^4*arcsinh 
(x*c)+e^4*(1/4*(c^2*x^2+1)^(1/2)*c^3*x^3-3/8*(c^2*x^2+1)^(1/2)*x*c+3/8*arc 
sinh(x*c))+4*d*c*e^3*(1/3*x^2*c^2*(c^2*x^2+1)^(1/2)-2/3*(c^2*x^2+1)^(1/2)) 
+6*d^2*c^2*e^2*(1/2*(c^2*x^2+1)^(1/2)*x*c-1/2*arcsinh(x*c))+4*d^3*c^3*e*(c 
^2*x^2+1)^(1/2))))
 

Fricas [A] (verification not implemented)

Time = 0.10 (sec) , antiderivative size = 475, normalized size of antiderivative = 1.29 \[ \int (d+e x)^3 (a+b \text {arcsinh}(c x))^2 \, dx=\frac {9 \, {\left (8 \, a^{2} + b^{2}\right )} c^{4} e^{3} x^{4} + 32 \, {\left (9 \, a^{2} + 2 \, b^{2}\right )} c^{4} d e^{2} x^{3} + 27 \, {\left (8 \, {\left (2 \, a^{2} + b^{2}\right )} c^{4} d^{2} e - b^{2} c^{2} e^{3}\right )} x^{2} + 9 \, {\left (8 \, b^{2} c^{4} e^{3} x^{4} + 32 \, b^{2} c^{4} d e^{2} x^{3} + 48 \, b^{2} c^{4} d^{2} e x^{2} + 32 \, b^{2} c^{4} d^{3} x + 24 \, b^{2} c^{2} d^{2} e - 3 \, b^{2} e^{3}\right )} \log \left (c x + \sqrt {c^{2} x^{2} + 1}\right )^{2} + 96 \, {\left (3 \, {\left (a^{2} + 2 \, b^{2}\right )} c^{4} d^{3} - 4 \, b^{2} c^{2} d e^{2}\right )} x + 6 \, {\left (24 \, a b c^{4} e^{3} x^{4} + 96 \, a b c^{4} d e^{2} x^{3} + 144 \, a b c^{4} d^{2} e x^{2} + 96 \, a b c^{4} d^{3} x + 72 \, a b c^{2} d^{2} e - 9 \, a b e^{3} - {\left (6 \, b^{2} c^{3} e^{3} x^{3} + 32 \, b^{2} c^{3} d e^{2} x^{2} + 96 \, b^{2} c^{3} d^{3} - 64 \, b^{2} c d e^{2} + 9 \, {\left (8 \, b^{2} c^{3} d^{2} e - b^{2} c e^{3}\right )} x\right )} \sqrt {c^{2} x^{2} + 1}\right )} \log \left (c x + \sqrt {c^{2} x^{2} + 1}\right ) - 6 \, {\left (6 \, a b c^{3} e^{3} x^{3} + 32 \, a b c^{3} d e^{2} x^{2} + 96 \, a b c^{3} d^{3} - 64 \, a b c d e^{2} + 9 \, {\left (8 \, a b c^{3} d^{2} e - a b c e^{3}\right )} x\right )} \sqrt {c^{2} x^{2} + 1}}{288 \, c^{4}} \] Input:

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

Output:

1/288*(9*(8*a^2 + b^2)*c^4*e^3*x^4 + 32*(9*a^2 + 2*b^2)*c^4*d*e^2*x^3 + 27 
*(8*(2*a^2 + b^2)*c^4*d^2*e - b^2*c^2*e^3)*x^2 + 9*(8*b^2*c^4*e^3*x^4 + 32 
*b^2*c^4*d*e^2*x^3 + 48*b^2*c^4*d^2*e*x^2 + 32*b^2*c^4*d^3*x + 24*b^2*c^2* 
d^2*e - 3*b^2*e^3)*log(c*x + sqrt(c^2*x^2 + 1))^2 + 96*(3*(a^2 + 2*b^2)*c^ 
4*d^3 - 4*b^2*c^2*d*e^2)*x + 6*(24*a*b*c^4*e^3*x^4 + 96*a*b*c^4*d*e^2*x^3 
+ 144*a*b*c^4*d^2*e*x^2 + 96*a*b*c^4*d^3*x + 72*a*b*c^2*d^2*e - 9*a*b*e^3 
- (6*b^2*c^3*e^3*x^3 + 32*b^2*c^3*d*e^2*x^2 + 96*b^2*c^3*d^3 - 64*b^2*c*d* 
e^2 + 9*(8*b^2*c^3*d^2*e - b^2*c*e^3)*x)*sqrt(c^2*x^2 + 1))*log(c*x + sqrt 
(c^2*x^2 + 1)) - 6*(6*a*b*c^3*e^3*x^3 + 32*a*b*c^3*d*e^2*x^2 + 96*a*b*c^3* 
d^3 - 64*a*b*c*d*e^2 + 9*(8*a*b*c^3*d^2*e - a*b*c*e^3)*x)*sqrt(c^2*x^2 + 1 
))/c^4
 

Sympy [B] (verification not implemented)

Leaf count of result is larger than twice the leaf count of optimal. 743 vs. \(2 (364) = 728\).

Time = 0.50 (sec) , antiderivative size = 743, normalized size of antiderivative = 2.02 \[ \int (d+e x)^3 (a+b \text {arcsinh}(c x))^2 \, dx =\text {Too large to display} \] Input:

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

Output:

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

Maxima [A] (verification not implemented)

Time = 0.07 (sec) , antiderivative size = 590, normalized size of antiderivative = 1.60 \[ \int (d+e x)^3 (a+b \text {arcsinh}(c x))^2 \, dx =\text {Too large to display} \] Input:

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

Output:

1/4*b^2*e^3*x^4*arcsinh(c*x)^2 + b^2*d*e^2*x^3*arcsinh(c*x)^2 + 1/4*a^2*e^ 
3*x^4 + 3/2*b^2*d^2*e*x^2*arcsinh(c*x)^2 + a^2*d*e^2*x^3 + b^2*d^3*x*arcsi 
nh(c*x)^2 + 3/2*a^2*d^2*e*x^2 + 3/2*(2*x^2*arcsinh(c*x) - c*(sqrt(c^2*x^2 
+ 1)*x/c^2 - arcsinh(c*x)/c^3))*a*b*d^2*e + 3/4*(c^2*(x^2/c^2 - log(c*x + 
sqrt(c^2*x^2 + 1))^2/c^4) - 2*c*(sqrt(c^2*x^2 + 1)*x/c^2 - arcsinh(c*x)/c^ 
3)*arcsinh(c*x))*b^2*d^2*e + 2/3*(3*x^3*arcsinh(c*x) - c*(sqrt(c^2*x^2 + 1 
)*x^2/c^2 - 2*sqrt(c^2*x^2 + 1)/c^4))*a*b*d*e^2 - 2/9*(3*c*(sqrt(c^2*x^2 + 
 1)*x^2/c^2 - 2*sqrt(c^2*x^2 + 1)/c^4)*arcsinh(c*x) - (c^2*x^3 - 6*x)/c^2) 
*b^2*d*e^2 + 1/16*(8*x^4*arcsinh(c*x) - (2*sqrt(c^2*x^2 + 1)*x^3/c^2 - 3*s 
qrt(c^2*x^2 + 1)*x/c^4 + 3*arcsinh(c*x)/c^5)*c)*a*b*e^3 + 1/32*((x^4/c^2 - 
 3*x^2/c^4 + 3*log(c*x + sqrt(c^2*x^2 + 1))^2/c^6)*c^2 - 2*(2*sqrt(c^2*x^2 
 + 1)*x^3/c^2 - 3*sqrt(c^2*x^2 + 1)*x/c^4 + 3*arcsinh(c*x)/c^5)*c*arcsinh( 
c*x))*b^2*e^3 + 2*b^2*d^3*(x - sqrt(c^2*x^2 + 1)*arcsinh(c*x)/c) + a^2*d^3 
*x + 2*(c*x*arcsinh(c*x) - sqrt(c^2*x^2 + 1))*a*b*d^3/c
 

Giac [F(-2)]

Exception generated. \[ \int (d+e x)^3 (a+b \text {arcsinh}(c x))^2 \, dx=\text {Exception raised: RuntimeError} \] Input:

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

Output:

Exception raised: RuntimeError >> an error occurred running a Giac command 
:INPUT:sage2OUTPUT:sym2poly/r2sym(const gen & e,const index_m & i,const ve 
cteur & l) Error: Bad Argument Value
 

Mupad [F(-1)]

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

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

Output:

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

Reduce [F]

\[ \int (d+e x)^3 (a+b \text {arcsinh}(c x))^2 \, dx=\frac {144 \left (\int \mathit {asinh} \left (c x \right )^{2} x^{2}d x \right ) b^{2} c^{4} d \,e^{2}+144 \mathit {asinh} \left (c x \right ) a b \,c^{4} d^{2} e \,x^{2}+96 \mathit {asinh} \left (c x \right ) a b \,c^{4} d \,e^{2} x^{3}-72 \sqrt {c^{2} x^{2}+1}\, a b \,c^{3} d^{2} e x -32 \sqrt {c^{2} x^{2}+1}\, a b \,c^{3} d \,e^{2} x^{2}+48 \mathit {asinh} \left (c x \right )^{2} b^{2} c^{4} d^{3} x +36 \mathit {asinh} \left (c x \right )^{2} b^{2} c^{2} d^{2} e -96 \sqrt {c^{2} x^{2}+1}\, \mathit {asinh} \left (c x \right ) b^{2} c^{3} d^{3}-96 \sqrt {c^{2} x^{2}+1}\, a b \,c^{3} d^{3}+72 a^{2} c^{4} d^{2} e \,x^{2}+48 a^{2} c^{4} d \,e^{2} x^{3}+36 b^{2} c^{4} d^{2} e \,x^{2}+48 \left (\int \mathit {asinh} \left (c x \right )^{2} x^{3}d x \right ) b^{2} c^{4} e^{3}-9 \,\mathrm {log}\left (\sqrt {c^{2} x^{2}+1}+c x \right ) a b \,e^{3}+48 a^{2} c^{4} d^{3} x +12 a^{2} c^{4} e^{3} x^{4}+96 b^{2} c^{4} d^{3} x -72 \sqrt {c^{2} x^{2}+1}\, \mathit {asinh} \left (c x \right ) b^{2} c^{3} d^{2} e x +72 \mathit {asinh} \left (c x \right )^{2} b^{2} c^{4} d^{2} e \,x^{2}+96 \mathit {asinh} \left (c x \right ) a b \,c^{4} d^{3} x +24 \mathit {asinh} \left (c x \right ) a b \,c^{4} e^{3} x^{4}-6 \sqrt {c^{2} x^{2}+1}\, a b \,c^{3} e^{3} x^{3}+64 \sqrt {c^{2} x^{2}+1}\, a b c d \,e^{2}+9 \sqrt {c^{2} x^{2}+1}\, a b c \,e^{3} x +72 \,\mathrm {log}\left (\sqrt {c^{2} x^{2}+1}+c x \right ) a b \,c^{2} d^{2} e}{48 c^{4}} \] Input:

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

Output:

(48*asinh(c*x)**2*b**2*c**4*d**3*x + 72*asinh(c*x)**2*b**2*c**4*d**2*e*x** 
2 + 36*asinh(c*x)**2*b**2*c**2*d**2*e - 96*sqrt(c**2*x**2 + 1)*asinh(c*x)* 
b**2*c**3*d**3 - 72*sqrt(c**2*x**2 + 1)*asinh(c*x)*b**2*c**3*d**2*e*x + 96 
*asinh(c*x)*a*b*c**4*d**3*x + 144*asinh(c*x)*a*b*c**4*d**2*e*x**2 + 96*asi 
nh(c*x)*a*b*c**4*d*e**2*x**3 + 24*asinh(c*x)*a*b*c**4*e**3*x**4 - 96*sqrt( 
c**2*x**2 + 1)*a*b*c**3*d**3 - 72*sqrt(c**2*x**2 + 1)*a*b*c**3*d**2*e*x - 
32*sqrt(c**2*x**2 + 1)*a*b*c**3*d*e**2*x**2 - 6*sqrt(c**2*x**2 + 1)*a*b*c* 
*3*e**3*x**3 + 64*sqrt(c**2*x**2 + 1)*a*b*c*d*e**2 + 9*sqrt(c**2*x**2 + 1) 
*a*b*c*e**3*x + 48*int(asinh(c*x)**2*x**3,x)*b**2*c**4*e**3 + 144*int(asin 
h(c*x)**2*x**2,x)*b**2*c**4*d*e**2 + 72*log(sqrt(c**2*x**2 + 1) + c*x)*a*b 
*c**2*d**2*e - 9*log(sqrt(c**2*x**2 + 1) + c*x)*a*b*e**3 + 48*a**2*c**4*d* 
*3*x + 72*a**2*c**4*d**2*e*x**2 + 48*a**2*c**4*d*e**2*x**3 + 12*a**2*c**4* 
e**3*x**4 + 96*b**2*c**4*d**3*x + 36*b**2*c**4*d**2*e*x**2)/(48*c**4)