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

3.3.30.1 Optimal result
3.3.30.2 Mathematica [C] (verified)
3.3.30.3 Rubi [A] (verified)
3.3.30.4 Maple [C] (verified)
3.3.30.5 Fricas [C] (verification not implemented)
3.3.30.6 Sympy [F]
3.3.30.7 Maxima [F(-2)]
3.3.30.8 Giac [F]
3.3.30.9 Mupad [F(-1)]

3.3.30.1 Optimal result

Integrand size = 23, antiderivative size = 261 \[ \int (c e+d e x)^{3/2} (a+b \text {arcsinh}(c+d x)) \, dx=-\frac {4 b (e (c+d x))^{3/2} \sqrt {1+(c+d x)^2}}{25 d}+\frac {12 b e \sqrt {e (c+d x)} \sqrt {1+(c+d x)^2}}{25 d (1+c+d x)}+\frac {2 (e (c+d x))^{5/2} (a+b \text {arcsinh}(c+d x))}{5 d e}-\frac {12 b e^{3/2} (1+c+d x) \sqrt {\frac {1+(c+d x)^2}{(1+c+d x)^2}} E\left (2 \arctan \left (\frac {\sqrt {e (c+d x)}}{\sqrt {e}}\right )|\frac {1}{2}\right )}{25 d \sqrt {1+(c+d x)^2}}+\frac {6 b e^{3/2} (1+c+d x) \sqrt {\frac {1+(c+d x)^2}{(1+c+d x)^2}} \operatorname {EllipticF}\left (2 \arctan \left (\frac {\sqrt {e (c+d x)}}{\sqrt {e}}\right ),\frac {1}{2}\right )}{25 d \sqrt {1+(c+d x)^2}} \]

output
2/5*(e*(d*x+c))^(5/2)*(a+b*arcsinh(d*x+c))/d/e-4/25*b*(e*(d*x+c))^(3/2)*(1 
+(d*x+c)^2)^(1/2)/d+12/25*b*e*(e*(d*x+c))^(1/2)*(1+(d*x+c)^2)^(1/2)/d/(d*x 
+c+1)-12/25*b*e^(3/2)*(d*x+c+1)*(cos(2*arctan((e*(d*x+c))^(1/2)/e^(1/2)))^ 
2)^(1/2)/cos(2*arctan((e*(d*x+c))^(1/2)/e^(1/2)))*EllipticE(sin(2*arctan(( 
e*(d*x+c))^(1/2)/e^(1/2))),1/2*2^(1/2))*((1+(d*x+c)^2)/(d*x+c+1)^2)^(1/2)/ 
d/(1+(d*x+c)^2)^(1/2)+6/25*b*e^(3/2)*(d*x+c+1)*(cos(2*arctan((e*(d*x+c))^( 
1/2)/e^(1/2)))^2)^(1/2)/cos(2*arctan((e*(d*x+c))^(1/2)/e^(1/2)))*EllipticF 
(sin(2*arctan((e*(d*x+c))^(1/2)/e^(1/2))),1/2*2^(1/2))*((1+(d*x+c)^2)/(d*x 
+c+1)^2)^(1/2)/d/(1+(d*x+c)^2)^(1/2)
 
3.3.30.2 Mathematica [C] (verified)

Result contains higher order function than in optimal. Order 5 vs. order 4 in optimal.

Time = 0.04 (sec) , antiderivative size = 87, normalized size of antiderivative = 0.33 \[ \int (c e+d e x)^{3/2} (a+b \text {arcsinh}(c+d x)) \, dx=\frac {2 (e (c+d x))^{3/2} \left (5 a c+5 a d x-2 b \sqrt {1+(c+d x)^2}+5 b c \text {arcsinh}(c+d x)+5 b d x \text {arcsinh}(c+d x)+2 b \operatorname {Hypergeometric2F1}\left (\frac {1}{2},\frac {3}{4},\frac {7}{4},-(c+d x)^2\right )\right )}{25 d} \]

input
Integrate[(c*e + d*e*x)^(3/2)*(a + b*ArcSinh[c + d*x]),x]
 
output
(2*(e*(c + d*x))^(3/2)*(5*a*c + 5*a*d*x - 2*b*Sqrt[1 + (c + d*x)^2] + 5*b* 
c*ArcSinh[c + d*x] + 5*b*d*x*ArcSinh[c + d*x] + 2*b*Hypergeometric2F1[1/2, 
 3/4, 7/4, -(c + d*x)^2]))/(25*d)
 
3.3.30.3 Rubi [A] (verified)

Time = 0.43 (sec) , antiderivative size = 286, normalized size of antiderivative = 1.10, number of steps used = 9, number of rules used = 8, \(\frac {\text {number of rules}}{\text {integrand size}}\) = 0.348, Rules used = {6274, 6191, 262, 266, 834, 27, 761, 1510}

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 (c e+d e x)^{3/2} (a+b \text {arcsinh}(c+d x)) \, dx\)

\(\Big \downarrow \) 6274

\(\displaystyle \frac {\int (e (c+d x))^{3/2} (a+b \text {arcsinh}(c+d x))d(c+d x)}{d}\)

\(\Big \downarrow \) 6191

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

\(\Big \downarrow \) 262

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

\(\Big \downarrow \) 266

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

\(\Big \downarrow \) 834

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

\(\Big \downarrow \) 27

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

\(\Big \downarrow \) 761

\(\displaystyle \frac {\frac {2 (e (c+d x))^{5/2} (a+b \text {arcsinh}(c+d x))}{5 e}-\frac {2 b \left (\frac {2}{5} e \sqrt {(c+d x)^2+1} (e (c+d x))^{3/2}-\frac {6}{5} e \left (\frac {\sqrt {e} (e (c+d x)+e) \sqrt {\frac {e^2 (c+d x)^2+e^2}{(e (c+d x)+e)^2}} \operatorname {EllipticF}\left (2 \arctan \left (\frac {\sqrt {e (c+d x)}}{\sqrt {e}}\right ),\frac {1}{2}\right )}{2 \sqrt {(c+d x)^2+1}}-\int \frac {e-e (c+d x)}{\sqrt {(c+d x)^2+1}}d\sqrt {e (c+d x)}\right )\right )}{5 e}}{d}\)

\(\Big \downarrow \) 1510

\(\displaystyle \frac {\frac {2 (e (c+d x))^{5/2} (a+b \text {arcsinh}(c+d x))}{5 e}-\frac {2 b \left (\frac {2}{5} e \sqrt {(c+d x)^2+1} (e (c+d x))^{3/2}-\frac {6}{5} e \left (\frac {\sqrt {e} (e (c+d x)+e) \sqrt {\frac {e^2 (c+d x)^2+e^2}{(e (c+d x)+e)^2}} \operatorname {EllipticF}\left (2 \arctan \left (\frac {\sqrt {e (c+d x)}}{\sqrt {e}}\right ),\frac {1}{2}\right )}{2 \sqrt {(c+d x)^2+1}}-\frac {\sqrt {e} (e (c+d x)+e) \sqrt {\frac {e^2 (c+d x)^2+e^2}{(e (c+d x)+e)^2}} E\left (2 \arctan \left (\frac {\sqrt {e (c+d x)}}{\sqrt {e}}\right )|\frac {1}{2}\right )}{\sqrt {(c+d x)^2+1}}+\frac {e^2 \sqrt {(c+d x)^2+1} \sqrt {e (c+d x)}}{e (c+d x)+e}\right )\right )}{5 e}}{d}\)

input
Int[(c*e + d*e*x)^(3/2)*(a + b*ArcSinh[c + d*x]),x]
 
output
((2*(e*(c + d*x))^(5/2)*(a + b*ArcSinh[c + d*x]))/(5*e) - (2*b*((2*e*(e*(c 
 + d*x))^(3/2)*Sqrt[1 + (c + d*x)^2])/5 - (6*e*((e^2*Sqrt[e*(c + d*x)]*Sqr 
t[1 + (c + d*x)^2])/(e + e*(c + d*x)) - (Sqrt[e]*(e + e*(c + d*x))*Sqrt[(e 
^2 + e^2*(c + d*x)^2)/(e + e*(c + d*x))^2]*EllipticE[2*ArcTan[Sqrt[e*(c + 
d*x)]/Sqrt[e]], 1/2])/Sqrt[1 + (c + d*x)^2] + (Sqrt[e]*(e + e*(c + d*x))*S 
qrt[(e^2 + e^2*(c + d*x)^2)/(e + e*(c + d*x))^2]*EllipticF[2*ArcTan[Sqrt[e 
*(c + d*x)]/Sqrt[e]], 1/2])/(2*Sqrt[1 + (c + d*x)^2])))/5))/(5*e))/d
 

3.3.30.3.1 Defintions of rubi rules used

rule 27
Int[(a_)*(Fx_), x_Symbol] :> Simp[a   Int[Fx, x], x] /; FreeQ[a, x] &&  !Ma 
tchQ[Fx, (b_)*(Gx_) /; FreeQ[b, x]]
 

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

rule 266
Int[((c_.)*(x_))^(m_)*((a_) + (b_.)*(x_)^2)^(p_), x_Symbol] :> With[{k = De 
nominator[m]}, Simp[k/c   Subst[Int[x^(k*(m + 1) - 1)*(a + b*(x^(2*k)/c^2)) 
^p, x], x, (c*x)^(1/k)], x]] /; FreeQ[{a, b, c, p}, x] && FractionQ[m] && I 
ntBinomialQ[a, b, c, 2, m, p, x]
 

rule 761
Int[1/Sqrt[(a_) + (b_.)*(x_)^4], x_Symbol] :> With[{q = Rt[b/a, 4]}, Simp[( 
1 + q^2*x^2)*(Sqrt[(a + b*x^4)/(a*(1 + q^2*x^2)^2)]/(2*q*Sqrt[a + b*x^4]))* 
EllipticF[2*ArcTan[q*x], 1/2], x]] /; FreeQ[{a, b}, x] && PosQ[b/a]
 

rule 834
Int[(x_)^2/Sqrt[(a_) + (b_.)*(x_)^4], x_Symbol] :> With[{q = Rt[b/a, 2]}, S 
imp[1/q   Int[1/Sqrt[a + b*x^4], x], x] - Simp[1/q   Int[(1 - q*x^2)/Sqrt[a 
 + b*x^4], x], x]] /; FreeQ[{a, b}, x] && PosQ[b/a]
 

rule 1510
Int[((d_) + (e_.)*(x_)^2)/Sqrt[(a_) + (c_.)*(x_)^4], x_Symbol] :> With[{q = 
 Rt[c/a, 4]}, Simp[(-d)*x*(Sqrt[a + c*x^4]/(a*(1 + q^2*x^2))), x] + Simp[d* 
(1 + q^2*x^2)*(Sqrt[(a + c*x^4)/(a*(1 + q^2*x^2)^2)]/(q*Sqrt[a + c*x^4]))*E 
llipticE[2*ArcTan[q*x], 1/2], x] /; EqQ[e + d*q^2, 0]] /; FreeQ[{a, c, d, e 
}, x] && PosQ[c/a]
 

rule 6191
Int[((a_.) + ArcSinh[(c_.)*(x_)]*(b_.))^(n_.)*((d_.)*(x_))^(m_.), x_Symbol] 
 :> Simp[(d*x)^(m + 1)*((a + b*ArcSinh[c*x])^n/(d*(m + 1))), x] - Simp[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 6274
Int[((a_.) + ArcSinh[(c_) + (d_.)*(x_)]*(b_.))^(n_.)*((e_.) + (f_.)*(x_))^( 
m_.), x_Symbol] :> Simp[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]
 
3.3.30.4 Maple [C] (verified)

Result contains complex when optimal does not.

Time = 0.91 (sec) , antiderivative size = 205, normalized size of antiderivative = 0.79

method result size
derivativedivides \(\frac {\frac {2 \left (d e x +c e \right )^{\frac {5}{2}} a}{5}+2 b \left (\frac {\left (d e x +c e \right )^{\frac {5}{2}} \operatorname {arcsinh}\left (\frac {d e x +c e}{e}\right )}{5}-\frac {2 \left (\frac {e^{2} \left (d e x +c e \right )^{\frac {3}{2}} \sqrt {\frac {\left (d e x +c e \right )^{2}}{e^{2}}+1}}{5}-\frac {3 i e^{3} \sqrt {1-\frac {i \left (d e x +c e \right )}{e}}\, \sqrt {1+\frac {i \left (d e x +c e \right )}{e}}\, \left (\operatorname {EllipticF}\left (\sqrt {d e x +c e}\, \sqrt {\frac {i}{e}}, i\right )-\operatorname {EllipticE}\left (\sqrt {d e x +c e}\, \sqrt {\frac {i}{e}}, i\right )\right )}{5 \sqrt {\frac {i}{e}}\, \sqrt {\frac {\left (d e x +c e \right )^{2}}{e^{2}}+1}}\right )}{5 e}\right )}{d e}\) \(205\)
default \(\frac {\frac {2 \left (d e x +c e \right )^{\frac {5}{2}} a}{5}+2 b \left (\frac {\left (d e x +c e \right )^{\frac {5}{2}} \operatorname {arcsinh}\left (\frac {d e x +c e}{e}\right )}{5}-\frac {2 \left (\frac {e^{2} \left (d e x +c e \right )^{\frac {3}{2}} \sqrt {\frac {\left (d e x +c e \right )^{2}}{e^{2}}+1}}{5}-\frac {3 i e^{3} \sqrt {1-\frac {i \left (d e x +c e \right )}{e}}\, \sqrt {1+\frac {i \left (d e x +c e \right )}{e}}\, \left (\operatorname {EllipticF}\left (\sqrt {d e x +c e}\, \sqrt {\frac {i}{e}}, i\right )-\operatorname {EllipticE}\left (\sqrt {d e x +c e}\, \sqrt {\frac {i}{e}}, i\right )\right )}{5 \sqrt {\frac {i}{e}}\, \sqrt {\frac {\left (d e x +c e \right )^{2}}{e^{2}}+1}}\right )}{5 e}\right )}{d e}\) \(205\)
parts \(\frac {2 a \left (d e x +c e \right )^{\frac {5}{2}}}{5 d e}+\frac {2 b \left (\frac {\left (d e x +c e \right )^{\frac {5}{2}} \operatorname {arcsinh}\left (\frac {d e x +c e}{e}\right )}{5}-\frac {2 \left (\frac {e^{2} \left (d e x +c e \right )^{\frac {3}{2}} \sqrt {\frac {\left (d e x +c e \right )^{2}}{e^{2}}+1}}{5}-\frac {3 i e^{3} \sqrt {1-\frac {i \left (d e x +c e \right )}{e}}\, \sqrt {1+\frac {i \left (d e x +c e \right )}{e}}\, \left (\operatorname {EllipticF}\left (\sqrt {d e x +c e}\, \sqrt {\frac {i}{e}}, i\right )-\operatorname {EllipticE}\left (\sqrt {d e x +c e}\, \sqrt {\frac {i}{e}}, i\right )\right )}{5 \sqrt {\frac {i}{e}}\, \sqrt {\frac {\left (d e x +c e \right )^{2}}{e^{2}}+1}}\right )}{5 e}\right )}{d e}\) \(210\)

input
int((d*e*x+c*e)^(3/2)*(a+b*arcsinh(d*x+c)),x,method=_RETURNVERBOSE)
 
output
2/d/e*(1/5*(d*e*x+c*e)^(5/2)*a+b*(1/5*(d*e*x+c*e)^(5/2)*arcsinh(1/e*(d*e*x 
+c*e))-2/5/e*(1/5*e^2*(d*e*x+c*e)^(3/2)*(1/e^2*(d*e*x+c*e)^2+1)^(1/2)-3/5* 
I*e^3/(I/e)^(1/2)*(1-I/e*(d*e*x+c*e))^(1/2)*(1+I/e*(d*e*x+c*e))^(1/2)/(1/e 
^2*(d*e*x+c*e)^2+1)^(1/2)*(EllipticF((d*e*x+c*e)^(1/2)*(I/e)^(1/2),I)-Elli 
pticE((d*e*x+c*e)^(1/2)*(I/e)^(1/2),I)))))
 
3.3.30.5 Fricas [C] (verification not implemented)

Result contains higher order function than in optimal. Order 9 vs. order 4.

Time = 0.10 (sec) , antiderivative size = 185, normalized size of antiderivative = 0.71 \[ \int (c e+d e x)^{3/2} (a+b \text {arcsinh}(c+d x)) \, dx=-\frac {2 \, {\left (6 \, \sqrt {d^{3} e} b e {\rm weierstrassZeta}\left (-\frac {4}{d^{2}}, 0, {\rm weierstrassPInverse}\left (-\frac {4}{d^{2}}, 0, \frac {d x + c}{d}\right )\right ) - 5 \, {\left (b d^{3} e x^{2} + 2 \, b c d^{2} e x + b c^{2} d e\right )} \sqrt {d e x + c e} \log \left (d x + c + \sqrt {d^{2} x^{2} + 2 \, c d x + c^{2} + 1}\right ) + 2 \, {\left (b d^{2} e x + b c d e\right )} \sqrt {d^{2} x^{2} + 2 \, c d x + c^{2} + 1} \sqrt {d e x + c e} - 5 \, {\left (a d^{3} e x^{2} + 2 \, a c d^{2} e x + a c^{2} d e\right )} \sqrt {d e x + c e}\right )}}{25 \, d^{2}} \]

input
integrate((d*e*x+c*e)^(3/2)*(a+b*arcsinh(d*x+c)),x, algorithm="fricas")
 
output
-2/25*(6*sqrt(d^3*e)*b*e*weierstrassZeta(-4/d^2, 0, weierstrassPInverse(-4 
/d^2, 0, (d*x + c)/d)) - 5*(b*d^3*e*x^2 + 2*b*c*d^2*e*x + b*c^2*d*e)*sqrt( 
d*e*x + c*e)*log(d*x + c + sqrt(d^2*x^2 + 2*c*d*x + c^2 + 1)) + 2*(b*d^2*e 
*x + b*c*d*e)*sqrt(d^2*x^2 + 2*c*d*x + c^2 + 1)*sqrt(d*e*x + c*e) - 5*(a*d 
^3*e*x^2 + 2*a*c*d^2*e*x + a*c^2*d*e)*sqrt(d*e*x + c*e))/d^2
 
3.3.30.6 Sympy [F]

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

input
integrate((d*e*x+c*e)**(3/2)*(a+b*asinh(d*x+c)),x)
 
output
Integral((e*(c + d*x))**(3/2)*(a + b*asinh(c + d*x)), x)
 
3.3.30.7 Maxima [F(-2)]

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

input
integrate((d*e*x+c*e)^(3/2)*(a+b*arcsinh(d*x+c)),x, algorithm="maxima")
 
output
Exception raised: ValueError >> Computation failed since Maxima requested 
additional constraints; using the 'assume' command before evaluation *may* 
 help (example of legal syntax is 'assume(e>0)', see `assume?` for more de 
tails)Is e
 
3.3.30.8 Giac [F]

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

input
integrate((d*e*x+c*e)^(3/2)*(a+b*arcsinh(d*x+c)),x, algorithm="giac")
 
output
integrate((d*e*x + c*e)^(3/2)*(b*arcsinh(d*x + c) + a), x)
 
3.3.30.9 Mupad [F(-1)]

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

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