\(\int (a+b \sinh (x))^{3/2} \, dx\) [106]

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

Optimal result

Integrand size = 10, antiderivative size = 150 \[ \int (a+b \sinh (x))^{3/2} \, dx=\frac {2}{3} b \cosh (x) \sqrt {a+b \sinh (x)}+\frac {8 i a E\left (\frac {\pi }{4}-\frac {i x}{2}|\frac {2 b}{i a+b}\right ) \sqrt {a+b \sinh (x)}}{3 \sqrt {\frac {a+b \sinh (x)}{a-i b}}}-\frac {2 i \left (a^2+b^2\right ) \operatorname {EllipticF}\left (\frac {\pi }{4}-\frac {i x}{2},\frac {2 b}{i a+b}\right ) \sqrt {\frac {a+b \sinh (x)}{a-i b}}}{3 \sqrt {a+b \sinh (x)}} \]

[Out]

2/3*b*cosh(x)*(a+b*sinh(x))^(1/2)+8/3*I*a*(sin(1/4*Pi+1/2*I*x)^2)^(1/2)/sin(1/4*Pi+1/2*I*x)*EllipticE(cos(1/4*
Pi+1/2*I*x),2^(1/2)*(b/(I*a+b))^(1/2))*(a+b*sinh(x))^(1/2)/((a+b*sinh(x))/(a-I*b))^(1/2)-2/3*I*(a^2+b^2)*(sin(
1/4*Pi+1/2*I*x)^2)^(1/2)/sin(1/4*Pi+1/2*I*x)*EllipticF(cos(1/4*Pi+1/2*I*x),2^(1/2)*(b/(I*a+b))^(1/2))*((a+b*si
nh(x))/(a-I*b))^(1/2)/(a+b*sinh(x))^(1/2)

Rubi [A] (verified)

Time = 0.12 (sec) , antiderivative size = 150, normalized size of antiderivative = 1.00, number of steps used = 6, number of rules used = 6, \(\frac {\text {number of rules}}{\text {integrand size}}\) = 0.600, Rules used = {2735, 2831, 2742, 2740, 2734, 2732} \[ \int (a+b \sinh (x))^{3/2} \, dx=-\frac {2 i \left (a^2+b^2\right ) \sqrt {\frac {a+b \sinh (x)}{a-i b}} \operatorname {EllipticF}\left (\frac {\pi }{4}-\frac {i x}{2},\frac {2 b}{i a+b}\right )}{3 \sqrt {a+b \sinh (x)}}+\frac {2}{3} b \cosh (x) \sqrt {a+b \sinh (x)}+\frac {8 i a \sqrt {a+b \sinh (x)} E\left (\frac {\pi }{4}-\frac {i x}{2}|\frac {2 b}{i a+b}\right )}{3 \sqrt {\frac {a+b \sinh (x)}{a-i b}}} \]

[In]

Int[(a + b*Sinh[x])^(3/2),x]

[Out]

(2*b*Cosh[x]*Sqrt[a + b*Sinh[x]])/3 + (((8*I)/3)*a*EllipticE[Pi/4 - (I/2)*x, (2*b)/(I*a + b)]*Sqrt[a + b*Sinh[
x]])/Sqrt[(a + b*Sinh[x])/(a - I*b)] - (((2*I)/3)*(a^2 + b^2)*EllipticF[Pi/4 - (I/2)*x, (2*b)/(I*a + b)]*Sqrt[
(a + b*Sinh[x])/(a - I*b)])/Sqrt[a + b*Sinh[x]]

Rule 2732

Int[Sqrt[(a_) + (b_.)*sin[(c_.) + (d_.)*(x_)]], x_Symbol] :> Simp[2*(Sqrt[a + b]/d)*EllipticE[(1/2)*(c - Pi/2
+ d*x), 2*(b/(a + b))], x] /; FreeQ[{a, b, c, d}, x] && NeQ[a^2 - b^2, 0] && GtQ[a + b, 0]

Rule 2734

Int[Sqrt[(a_) + (b_.)*sin[(c_.) + (d_.)*(x_)]], x_Symbol] :> Dist[Sqrt[a + b*Sin[c + d*x]]/Sqrt[(a + b*Sin[c +
 d*x])/(a + b)], Int[Sqrt[a/(a + b) + (b/(a + b))*Sin[c + d*x]], x], x] /; FreeQ[{a, b, c, d}, x] && NeQ[a^2 -
 b^2, 0] &&  !GtQ[a + b, 0]

Rule 2735

Int[((a_) + (b_.)*sin[(c_.) + (d_.)*(x_)])^(n_), x_Symbol] :> Simp[(-b)*Cos[c + d*x]*((a + b*Sin[c + d*x])^(n
- 1)/(d*n)), x] + Dist[1/n, Int[(a + b*Sin[c + d*x])^(n - 2)*Simp[a^2*n + b^2*(n - 1) + a*b*(2*n - 1)*Sin[c +
d*x], x], x], x] /; FreeQ[{a, b, c, d}, x] && NeQ[a^2 - b^2, 0] && GtQ[n, 1] && IntegerQ[2*n]

Rule 2740

Int[1/Sqrt[(a_) + (b_.)*sin[(c_.) + (d_.)*(x_)]], x_Symbol] :> Simp[(2/(d*Sqrt[a + b]))*EllipticF[(1/2)*(c - P
i/2 + d*x), 2*(b/(a + b))], x] /; FreeQ[{a, b, c, d}, x] && NeQ[a^2 - b^2, 0] && GtQ[a + b, 0]

Rule 2742

Int[1/Sqrt[(a_) + (b_.)*sin[(c_.) + (d_.)*(x_)]], x_Symbol] :> Dist[Sqrt[(a + b*Sin[c + d*x])/(a + b)]/Sqrt[a
+ b*Sin[c + d*x]], Int[1/Sqrt[a/(a + b) + (b/(a + b))*Sin[c + d*x]], x], x] /; FreeQ[{a, b, c, d}, x] && NeQ[a
^2 - b^2, 0] &&  !GtQ[a + b, 0]

Rule 2831

Int[((c_.) + (d_.)*sin[(e_.) + (f_.)*(x_)])/Sqrt[(a_) + (b_.)*sin[(e_.) + (f_.)*(x_)]], x_Symbol] :> Dist[(b*c
 - a*d)/b, Int[1/Sqrt[a + b*Sin[e + f*x]], x], x] + Dist[d/b, Int[Sqrt[a + b*Sin[e + f*x]], x], x] /; FreeQ[{a
, b, c, d, e, f}, x] && NeQ[b*c - a*d, 0] && NeQ[a^2 - b^2, 0]

Rubi steps \begin{align*} \text {integral}& = \frac {2}{3} b \cosh (x) \sqrt {a+b \sinh (x)}+\frac {2}{3} \int \frac {\frac {1}{2} \left (3 a^2-b^2\right )+2 a b \sinh (x)}{\sqrt {a+b \sinh (x)}} \, dx \\ & = \frac {2}{3} b \cosh (x) \sqrt {a+b \sinh (x)}+\frac {1}{3} (4 a) \int \sqrt {a+b \sinh (x)} \, dx+\frac {1}{3} \left (-a^2-b^2\right ) \int \frac {1}{\sqrt {a+b \sinh (x)}} \, dx \\ & = \frac {2}{3} b \cosh (x) \sqrt {a+b \sinh (x)}+\frac {\left (4 a \sqrt {a+b \sinh (x)}\right ) \int \sqrt {\frac {a}{a-i b}+\frac {b \sinh (x)}{a-i b}} \, dx}{3 \sqrt {\frac {a+b \sinh (x)}{a-i b}}}+\frac {\left (\left (-a^2-b^2\right ) \sqrt {\frac {a+b \sinh (x)}{a-i b}}\right ) \int \frac {1}{\sqrt {\frac {a}{a-i b}+\frac {b \sinh (x)}{a-i b}}} \, dx}{3 \sqrt {a+b \sinh (x)}} \\ & = \frac {2}{3} b \cosh (x) \sqrt {a+b \sinh (x)}+\frac {8 i a E\left (\frac {\pi }{4}-\frac {i x}{2}|\frac {2 b}{i a+b}\right ) \sqrt {a+b \sinh (x)}}{3 \sqrt {\frac {a+b \sinh (x)}{a-i b}}}-\frac {2 i \left (a^2+b^2\right ) \operatorname {EllipticF}\left (\frac {\pi }{4}-\frac {i x}{2},\frac {2 b}{i a+b}\right ) \sqrt {\frac {a+b \sinh (x)}{a-i b}}}{3 \sqrt {a+b \sinh (x)}} \\ \end{align*}

Mathematica [A] (verified)

Time = 0.27 (sec) , antiderivative size = 139, normalized size of antiderivative = 0.93 \[ \int (a+b \sinh (x))^{3/2} \, dx=\frac {2 b \cosh (x) (a+b \sinh (x))+8 a (i a+b) E\left (\frac {1}{4} (\pi -2 i x)|-\frac {2 i b}{a-i b}\right ) \sqrt {\frac {a+b \sinh (x)}{a-i b}}-2 i \left (a^2+b^2\right ) \operatorname {EllipticF}\left (\frac {1}{4} (\pi -2 i x),-\frac {2 i b}{a-i b}\right ) \sqrt {\frac {a+b \sinh (x)}{a-i b}}}{3 \sqrt {a+b \sinh (x)}} \]

[In]

Integrate[(a + b*Sinh[x])^(3/2),x]

[Out]

(2*b*Cosh[x]*(a + b*Sinh[x]) + 8*a*(I*a + b)*EllipticE[(Pi - (2*I)*x)/4, ((-2*I)*b)/(a - I*b)]*Sqrt[(a + b*Sin
h[x])/(a - I*b)] - (2*I)*(a^2 + b^2)*EllipticF[(Pi - (2*I)*x)/4, ((-2*I)*b)/(a - I*b)]*Sqrt[(a + b*Sinh[x])/(a
 - I*b)])/(3*Sqrt[a + b*Sinh[x]])

Maple [B] (verified)

Both result and optimal contain complex but leaf count of result is larger than twice the leaf count of optimal. 675 vs. \(2 (178 ) = 356\).

Time = 1.68 (sec) , antiderivative size = 676, normalized size of antiderivative = 4.51

method result size
default \(\frac {\frac {2 i \sqrt {-\frac {a +b \sinh \left (x \right )}{i b -a}}\, \sqrt {\frac {\left (i-\sinh \left (x \right )\right ) b}{i b +a}}\, \sqrt {\frac {\left (i+\sinh \left (x \right )\right ) b}{i b -a}}\, \operatorname {EllipticF}\left (\sqrt {-\frac {a +b \sinh \left (x \right )}{i b -a}}, \sqrt {-\frac {i b -a}{i b +a}}\right ) a^{2} b}{3}+\frac {2 i \sqrt {-\frac {a +b \sinh \left (x \right )}{i b -a}}\, \sqrt {\frac {\left (i-\sinh \left (x \right )\right ) b}{i b +a}}\, \sqrt {\frac {\left (i+\sinh \left (x \right )\right ) b}{i b -a}}\, \operatorname {EllipticF}\left (\sqrt {-\frac {a +b \sinh \left (x \right )}{i b -a}}, \sqrt {-\frac {i b -a}{i b +a}}\right ) b^{3}}{3}+2 \sqrt {-\frac {a +b \sinh \left (x \right )}{i b -a}}\, \sqrt {\frac {\left (i-\sinh \left (x \right )\right ) b}{i b +a}}\, \sqrt {\frac {\left (i+\sinh \left (x \right )\right ) b}{i b -a}}\, \operatorname {EllipticF}\left (\sqrt {-\frac {a +b \sinh \left (x \right )}{i b -a}}, \sqrt {-\frac {i b -a}{i b +a}}\right ) a^{3}+2 \sqrt {-\frac {a +b \sinh \left (x \right )}{i b -a}}\, \sqrt {\frac {\left (i-\sinh \left (x \right )\right ) b}{i b +a}}\, \sqrt {\frac {\left (i+\sinh \left (x \right )\right ) b}{i b -a}}\, \operatorname {EllipticF}\left (\sqrt {-\frac {a +b \sinh \left (x \right )}{i b -a}}, \sqrt {-\frac {i b -a}{i b +a}}\right ) a \,b^{2}-\frac {8 \sqrt {-\frac {a +b \sinh \left (x \right )}{i b -a}}\, \sqrt {\frac {\left (i-\sinh \left (x \right )\right ) b}{i b +a}}\, \sqrt {\frac {\left (i+\sinh \left (x \right )\right ) b}{i b -a}}\, \operatorname {EllipticE}\left (\sqrt {-\frac {a +b \sinh \left (x \right )}{i b -a}}, \sqrt {-\frac {i b -a}{i b +a}}\right ) a^{3}}{3}-\frac {8 \sqrt {-\frac {a +b \sinh \left (x \right )}{i b -a}}\, \sqrt {\frac {\left (i-\sinh \left (x \right )\right ) b}{i b +a}}\, \sqrt {\frac {\left (i+\sinh \left (x \right )\right ) b}{i b -a}}\, \operatorname {EllipticE}\left (\sqrt {-\frac {a +b \sinh \left (x \right )}{i b -a}}, \sqrt {-\frac {i b -a}{i b +a}}\right ) a \,b^{2}}{3}+\frac {2 b^{3} \sinh \left (x \right )^{3}}{3}+\frac {2 a \,b^{2} \sinh \left (x \right )^{2}}{3}+\frac {2 b^{3} \sinh \left (x \right )}{3}+\frac {2 a \,b^{2}}{3}}{b \cosh \left (x \right ) \sqrt {a +b \sinh \left (x \right )}}\) \(676\)

[In]

int((a+b*sinh(x))^(3/2),x,method=_RETURNVERBOSE)

[Out]

2/3*(I*(-(a+b*sinh(x))/(I*b-a))^(1/2)*((I-sinh(x))*b/(I*b+a))^(1/2)*((I+sinh(x))*b/(I*b-a))^(1/2)*EllipticF((-
(a+b*sinh(x))/(I*b-a))^(1/2),(-(I*b-a)/(I*b+a))^(1/2))*a^2*b+I*(-(a+b*sinh(x))/(I*b-a))^(1/2)*((I-sinh(x))*b/(
I*b+a))^(1/2)*((I+sinh(x))*b/(I*b-a))^(1/2)*EllipticF((-(a+b*sinh(x))/(I*b-a))^(1/2),(-(I*b-a)/(I*b+a))^(1/2))
*b^3+3*(-(a+b*sinh(x))/(I*b-a))^(1/2)*((I-sinh(x))*b/(I*b+a))^(1/2)*((I+sinh(x))*b/(I*b-a))^(1/2)*EllipticF((-
(a+b*sinh(x))/(I*b-a))^(1/2),(-(I*b-a)/(I*b+a))^(1/2))*a^3+3*(-(a+b*sinh(x))/(I*b-a))^(1/2)*((I-sinh(x))*b/(I*
b+a))^(1/2)*((I+sinh(x))*b/(I*b-a))^(1/2)*EllipticF((-(a+b*sinh(x))/(I*b-a))^(1/2),(-(I*b-a)/(I*b+a))^(1/2))*a
*b^2-4*(-(a+b*sinh(x))/(I*b-a))^(1/2)*((I-sinh(x))*b/(I*b+a))^(1/2)*((I+sinh(x))*b/(I*b-a))^(1/2)*EllipticE((-
(a+b*sinh(x))/(I*b-a))^(1/2),(-(I*b-a)/(I*b+a))^(1/2))*a^3-4*(-(a+b*sinh(x))/(I*b-a))^(1/2)*((I-sinh(x))*b/(I*
b+a))^(1/2)*((I+sinh(x))*b/(I*b-a))^(1/2)*EllipticE((-(a+b*sinh(x))/(I*b-a))^(1/2),(-(I*b-a)/(I*b+a))^(1/2))*a
*b^2+b^3*sinh(x)^3+a*b^2*sinh(x)^2+b^3*sinh(x)+a*b^2)/b/cosh(x)/(a+b*sinh(x))^(1/2)

Fricas [C] (verification not implemented)

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

Time = 0.11 (sec) , antiderivative size = 263, normalized size of antiderivative = 1.75 \[ \int (a+b \sinh (x))^{3/2} \, dx=\frac {2 \, {\left (\sqrt {2} {\left (a^{2} - 3 \, b^{2}\right )} \cosh \left (x\right ) + \sqrt {2} {\left (a^{2} - 3 \, b^{2}\right )} \sinh \left (x\right )\right )} \sqrt {b} {\rm weierstrassPInverse}\left (\frac {4 \, {\left (4 \, a^{2} + 3 \, b^{2}\right )}}{3 \, b^{2}}, -\frac {8 \, {\left (8 \, a^{3} + 9 \, a b^{2}\right )}}{27 \, b^{3}}, \frac {3 \, b \cosh \left (x\right ) + 3 \, b \sinh \left (x\right ) + 2 \, a}{3 \, b}\right ) - 24 \, {\left (\sqrt {2} a b \cosh \left (x\right ) + \sqrt {2} a b \sinh \left (x\right )\right )} \sqrt {b} {\rm weierstrassZeta}\left (\frac {4 \, {\left (4 \, a^{2} + 3 \, b^{2}\right )}}{3 \, b^{2}}, -\frac {8 \, {\left (8 \, a^{3} + 9 \, a b^{2}\right )}}{27 \, b^{3}}, {\rm weierstrassPInverse}\left (\frac {4 \, {\left (4 \, a^{2} + 3 \, b^{2}\right )}}{3 \, b^{2}}, -\frac {8 \, {\left (8 \, a^{3} + 9 \, a b^{2}\right )}}{27 \, b^{3}}, \frac {3 \, b \cosh \left (x\right ) + 3 \, b \sinh \left (x\right ) + 2 \, a}{3 \, b}\right )\right ) + 3 \, {\left (b^{2} \cosh \left (x\right )^{2} + b^{2} \sinh \left (x\right )^{2} - 8 \, a b \cosh \left (x\right ) + b^{2} + 2 \, {\left (b^{2} \cosh \left (x\right ) - 4 \, a b\right )} \sinh \left (x\right )\right )} \sqrt {b \sinh \left (x\right ) + a}}{9 \, {\left (b \cosh \left (x\right ) + b \sinh \left (x\right )\right )}} \]

[In]

integrate((a+b*sinh(x))^(3/2),x, algorithm="fricas")

[Out]

1/9*(2*(sqrt(2)*(a^2 - 3*b^2)*cosh(x) + sqrt(2)*(a^2 - 3*b^2)*sinh(x))*sqrt(b)*weierstrassPInverse(4/3*(4*a^2
+ 3*b^2)/b^2, -8/27*(8*a^3 + 9*a*b^2)/b^3, 1/3*(3*b*cosh(x) + 3*b*sinh(x) + 2*a)/b) - 24*(sqrt(2)*a*b*cosh(x)
+ sqrt(2)*a*b*sinh(x))*sqrt(b)*weierstrassZeta(4/3*(4*a^2 + 3*b^2)/b^2, -8/27*(8*a^3 + 9*a*b^2)/b^3, weierstra
ssPInverse(4/3*(4*a^2 + 3*b^2)/b^2, -8/27*(8*a^3 + 9*a*b^2)/b^3, 1/3*(3*b*cosh(x) + 3*b*sinh(x) + 2*a)/b)) + 3
*(b^2*cosh(x)^2 + b^2*sinh(x)^2 - 8*a*b*cosh(x) + b^2 + 2*(b^2*cosh(x) - 4*a*b)*sinh(x))*sqrt(b*sinh(x) + a))/
(b*cosh(x) + b*sinh(x))

Sympy [F]

\[ \int (a+b \sinh (x))^{3/2} \, dx=\int \left (a + b \sinh {\left (x \right )}\right )^{\frac {3}{2}}\, dx \]

[In]

integrate((a+b*sinh(x))**(3/2),x)

[Out]

Integral((a + b*sinh(x))**(3/2), x)

Maxima [F]

\[ \int (a+b \sinh (x))^{3/2} \, dx=\int { {\left (b \sinh \left (x\right ) + a\right )}^{\frac {3}{2}} \,d x } \]

[In]

integrate((a+b*sinh(x))^(3/2),x, algorithm="maxima")

[Out]

integrate((b*sinh(x) + a)^(3/2), x)

Giac [F]

\[ \int (a+b \sinh (x))^{3/2} \, dx=\int { {\left (b \sinh \left (x\right ) + a\right )}^{\frac {3}{2}} \,d x } \]

[In]

integrate((a+b*sinh(x))^(3/2),x, algorithm="giac")

[Out]

integrate((b*sinh(x) + a)^(3/2), x)

Mupad [F(-1)]

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

[In]

int((a + b*sinh(x))^(3/2),x)

[Out]

int((a + b*sinh(x))^(3/2), x)