\(\int \frac {1}{(a+b \sinh (c+d x))^3} \, dx\) [103]

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

Optimal result

Integrand size = 12, antiderivative size = 127 \[ \int \frac {1}{(a+b \sinh (c+d x))^3} \, dx=-\frac {\left (2 a^2-b^2\right ) \text {arctanh}\left (\frac {b-a \tanh \left (\frac {1}{2} (c+d x)\right )}{\sqrt {a^2+b^2}}\right )}{\left (a^2+b^2\right )^{5/2} d}-\frac {b \cosh (c+d x)}{2 \left (a^2+b^2\right ) d (a+b \sinh (c+d x))^2}-\frac {3 a b \cosh (c+d x)}{2 \left (a^2+b^2\right )^2 d (a+b \sinh (c+d x))} \] Output:

-(2*a^2-b^2)*arctanh((b-a*tanh(1/2*d*x+1/2*c))/(a^2+b^2)^(1/2))/(a^2+b^2)^ 
(5/2)/d-1/2*b*cosh(d*x+c)/(a^2+b^2)/d/(a+b*sinh(d*x+c))^2-3/2*a*b*cosh(d*x 
+c)/(a^2+b^2)^2/d/(a+b*sinh(d*x+c))
 

Mathematica [A] (verified)

Time = 0.23 (sec) , antiderivative size = 117, normalized size of antiderivative = 0.92 \[ \int \frac {1}{(a+b \sinh (c+d x))^3} \, dx=\frac {\frac {2 \left (2 a^2-b^2\right ) \arctan \left (\frac {b-a \tanh \left (\frac {1}{2} (c+d x)\right )}{\sqrt {-a^2-b^2}}\right )}{\sqrt {-a^2-b^2}}-\frac {b \cosh (c+d x) \left (4 a^2+b^2+3 a b \sinh (c+d x)\right )}{(a+b \sinh (c+d x))^2}}{2 \left (a^2+b^2\right )^2 d} \] Input:

Integrate[(a + b*Sinh[c + d*x])^(-3),x]
 

Output:

((2*(2*a^2 - b^2)*ArcTan[(b - a*Tanh[(c + d*x)/2])/Sqrt[-a^2 - b^2]])/Sqrt 
[-a^2 - b^2] - (b*Cosh[c + d*x]*(4*a^2 + b^2 + 3*a*b*Sinh[c + d*x]))/(a + 
b*Sinh[c + d*x])^2)/(2*(a^2 + b^2)^2*d)
 

Rubi [A] (warning: unable to verify)

Time = 0.58 (sec) , antiderivative size = 137, normalized size of antiderivative = 1.08, number of steps used = 12, number of rules used = 11, \(\frac {\text {number of rules}}{\text {integrand size}}\) = 0.917, Rules used = {3042, 3143, 25, 3042, 3233, 25, 27, 3042, 3139, 1083, 217}

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 \frac {1}{(a+b \sinh (c+d x))^3} \, dx\)

\(\Big \downarrow \) 3042

\(\displaystyle \int \frac {1}{(a-i b \sin (i c+i d x))^3}dx\)

\(\Big \downarrow \) 3143

\(\displaystyle -\frac {\int -\frac {2 a-b \sinh (c+d x)}{(a+b \sinh (c+d x))^2}dx}{2 \left (a^2+b^2\right )}-\frac {b \cosh (c+d x)}{2 d \left (a^2+b^2\right ) (a+b \sinh (c+d x))^2}\)

\(\Big \downarrow \) 25

\(\displaystyle \frac {\int \frac {2 a-b \sinh (c+d x)}{(a+b \sinh (c+d x))^2}dx}{2 \left (a^2+b^2\right )}-\frac {b \cosh (c+d x)}{2 d \left (a^2+b^2\right ) (a+b \sinh (c+d x))^2}\)

\(\Big \downarrow \) 3042

\(\displaystyle -\frac {b \cosh (c+d x)}{2 d \left (a^2+b^2\right ) (a+b \sinh (c+d x))^2}+\frac {\int \frac {2 a+i b \sin (i c+i d x)}{(a-i b \sin (i c+i d x))^2}dx}{2 \left (a^2+b^2\right )}\)

\(\Big \downarrow \) 3233

\(\displaystyle \frac {-\frac {\int -\frac {2 a^2-b^2}{a+b \sinh (c+d x)}dx}{a^2+b^2}-\frac {3 a b \cosh (c+d x)}{d \left (a^2+b^2\right ) (a+b \sinh (c+d x))}}{2 \left (a^2+b^2\right )}-\frac {b \cosh (c+d x)}{2 d \left (a^2+b^2\right ) (a+b \sinh (c+d x))^2}\)

\(\Big \downarrow \) 25

\(\displaystyle \frac {\frac {\int \frac {2 a^2-b^2}{a+b \sinh (c+d x)}dx}{a^2+b^2}-\frac {3 a b \cosh (c+d x)}{d \left (a^2+b^2\right ) (a+b \sinh (c+d x))}}{2 \left (a^2+b^2\right )}-\frac {b \cosh (c+d x)}{2 d \left (a^2+b^2\right ) (a+b \sinh (c+d x))^2}\)

\(\Big \downarrow \) 27

\(\displaystyle \frac {\frac {\left (2 a^2-b^2\right ) \int \frac {1}{a+b \sinh (c+d x)}dx}{a^2+b^2}-\frac {3 a b \cosh (c+d x)}{d \left (a^2+b^2\right ) (a+b \sinh (c+d x))}}{2 \left (a^2+b^2\right )}-\frac {b \cosh (c+d x)}{2 d \left (a^2+b^2\right ) (a+b \sinh (c+d x))^2}\)

\(\Big \downarrow \) 3042

\(\displaystyle -\frac {b \cosh (c+d x)}{2 d \left (a^2+b^2\right ) (a+b \sinh (c+d x))^2}+\frac {-\frac {3 a b \cosh (c+d x)}{d \left (a^2+b^2\right ) (a+b \sinh (c+d x))}+\frac {\left (2 a^2-b^2\right ) \int \frac {1}{a-i b \sin (i c+i d x)}dx}{a^2+b^2}}{2 \left (a^2+b^2\right )}\)

\(\Big \downarrow \) 3139

\(\displaystyle -\frac {b \cosh (c+d x)}{2 d \left (a^2+b^2\right ) (a+b \sinh (c+d x))^2}+\frac {-\frac {3 a b \cosh (c+d x)}{d \left (a^2+b^2\right ) (a+b \sinh (c+d x))}-\frac {2 i \left (2 a^2-b^2\right ) \int \frac {1}{-a \tanh ^2\left (\frac {1}{2} (c+d x)\right )+2 b \tanh \left (\frac {1}{2} (c+d x)\right )+a}d\left (i \tanh \left (\frac {1}{2} (c+d x)\right )\right )}{d \left (a^2+b^2\right )}}{2 \left (a^2+b^2\right )}\)

\(\Big \downarrow \) 1083

\(\displaystyle -\frac {b \cosh (c+d x)}{2 d \left (a^2+b^2\right ) (a+b \sinh (c+d x))^2}+\frac {-\frac {3 a b \cosh (c+d x)}{d \left (a^2+b^2\right ) (a+b \sinh (c+d x))}+\frac {4 i \left (2 a^2-b^2\right ) \int \frac {1}{\tanh ^2\left (\frac {1}{2} (c+d x)\right )-4 \left (a^2+b^2\right )}d\left (2 i a \tanh \left (\frac {1}{2} (c+d x)\right )-2 i b\right )}{d \left (a^2+b^2\right )}}{2 \left (a^2+b^2\right )}\)

\(\Big \downarrow \) 217

\(\displaystyle \frac {\frac {2 \left (2 a^2-b^2\right ) \text {arctanh}\left (\frac {\tanh \left (\frac {1}{2} (c+d x)\right )}{2 \sqrt {a^2+b^2}}\right )}{d \left (a^2+b^2\right )^{3/2}}-\frac {3 a b \cosh (c+d x)}{d \left (a^2+b^2\right ) (a+b \sinh (c+d x))}}{2 \left (a^2+b^2\right )}-\frac {b \cosh (c+d x)}{2 d \left (a^2+b^2\right ) (a+b \sinh (c+d x))^2}\)

Input:

Int[(a + b*Sinh[c + d*x])^(-3),x]
 

Output:

-1/2*(b*Cosh[c + d*x])/((a^2 + b^2)*d*(a + b*Sinh[c + d*x])^2) + ((2*(2*a^ 
2 - b^2)*ArcTanh[Tanh[(c + d*x)/2]/(2*Sqrt[a^2 + b^2])])/((a^2 + b^2)^(3/2 
)*d) - (3*a*b*Cosh[c + d*x])/((a^2 + b^2)*d*(a + b*Sinh[c + d*x])))/(2*(a^ 
2 + b^2))
 

Defintions of rubi rules used

rule 25
Int[-(Fx_), x_Symbol] :> Simp[Identity[-1]   Int[Fx, x], x]
 

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 217
Int[((a_) + (b_.)*(x_)^2)^(-1), x_Symbol] :> Simp[(-(Rt[-a, 2]*Rt[-b, 2])^( 
-1))*ArcTan[Rt[-b, 2]*(x/Rt[-a, 2])], x] /; FreeQ[{a, b}, x] && PosQ[a/b] & 
& (LtQ[a, 0] || LtQ[b, 0])
 

rule 1083
Int[((a_) + (b_.)*(x_) + (c_.)*(x_)^2)^(-1), x_Symbol] :> Simp[-2   Subst[I 
nt[1/Simp[b^2 - 4*a*c - x^2, x], x], x, b + 2*c*x], x] /; FreeQ[{a, b, c}, 
x]
 

rule 3042
Int[u_, x_Symbol] :> Int[DeactivateTrig[u, x], x] /; FunctionOfTrigOfLinear 
Q[u, x]
 

rule 3139
Int[((a_) + (b_.)*sin[(c_.) + (d_.)*(x_)])^(-1), x_Symbol] :> With[{e = Fre 
eFactors[Tan[(c + d*x)/2], x]}, Simp[2*(e/d)   Subst[Int[1/(a + 2*b*e*x + a 
*e^2*x^2), x], x, Tan[(c + d*x)/2]/e], x]] /; FreeQ[{a, b, c, d}, x] && NeQ 
[a^2 - b^2, 0]
 

rule 3143
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 + 1)*(a^2 - b^2))), x] + Simp 
[1/((n + 1)*(a^2 - b^2))   Int[(a + b*Sin[c + d*x])^(n + 1)*Simp[a*(n + 1) 
- b*(n + 2)*Sin[c + d*x], x], x], x] /; FreeQ[{a, b, c, d}, x] && NeQ[a^2 - 
 b^2, 0] && LtQ[n, -1] && IntegerQ[2*n]
 

rule 3233
Int[((a_) + (b_.)*sin[(e_.) + (f_.)*(x_)])^(m_)*((c_.) + (d_.)*sin[(e_.) + 
(f_.)*(x_)]), x_Symbol] :> Simp[(-(b*c - a*d))*Cos[e + f*x]*((a + b*Sin[e + 
 f*x])^(m + 1)/(f*(m + 1)*(a^2 - b^2))), x] + Simp[1/((m + 1)*(a^2 - b^2)) 
  Int[(a + b*Sin[e + f*x])^(m + 1)*Simp[(a*c - b*d)*(m + 1) - (b*c - a*d)*( 
m + 2)*Sin[e + f*x], x], x], x] /; FreeQ[{a, b, c, d, e, f}, x] && NeQ[b*c 
- a*d, 0] && NeQ[a^2 - b^2, 0] && LtQ[m, -1] && IntegerQ[2*m]
 
Maple [B] (verified)

Leaf count of result is larger than twice the leaf count of optimal. \(279\) vs. \(2(118)=236\).

Time = 0.36 (sec) , antiderivative size = 280, normalized size of antiderivative = 2.20

method result size
derivativedivides \(\frac {-\frac {2 \left (-\frac {b^{2} \left (5 a^{2}+2 b^{2}\right ) \tanh \left (\frac {d x}{2}+\frac {c}{2}\right )^{3}}{2 a \left (a^{4}+2 a^{2} b^{2}+b^{4}\right )}-\frac {b \left (4 a^{4}-7 a^{2} b^{2}-2 b^{4}\right ) \tanh \left (\frac {d x}{2}+\frac {c}{2}\right )^{2}}{2 \left (a^{4}+2 a^{2} b^{2}+b^{4}\right ) a^{2}}+\frac {b^{2} \left (11 a^{2}+2 b^{2}\right ) \tanh \left (\frac {d x}{2}+\frac {c}{2}\right )}{2 \left (a^{4}+2 a^{2} b^{2}+b^{4}\right ) a}+\frac {b \left (4 a^{2}+b^{2}\right )}{2 a^{4}+4 a^{2} b^{2}+2 b^{4}}\right )}{\left (\tanh \left (\frac {d x}{2}+\frac {c}{2}\right )^{2} a -2 b \tanh \left (\frac {d x}{2}+\frac {c}{2}\right )-a \right )^{2}}+\frac {\left (2 a^{2}-b^{2}\right ) \operatorname {arctanh}\left (\frac {2 a \tanh \left (\frac {d x}{2}+\frac {c}{2}\right )-2 b}{2 \sqrt {a^{2}+b^{2}}}\right )}{\left (a^{4}+2 a^{2} b^{2}+b^{4}\right ) \sqrt {a^{2}+b^{2}}}}{d}\) \(280\)
default \(\frac {-\frac {2 \left (-\frac {b^{2} \left (5 a^{2}+2 b^{2}\right ) \tanh \left (\frac {d x}{2}+\frac {c}{2}\right )^{3}}{2 a \left (a^{4}+2 a^{2} b^{2}+b^{4}\right )}-\frac {b \left (4 a^{4}-7 a^{2} b^{2}-2 b^{4}\right ) \tanh \left (\frac {d x}{2}+\frac {c}{2}\right )^{2}}{2 \left (a^{4}+2 a^{2} b^{2}+b^{4}\right ) a^{2}}+\frac {b^{2} \left (11 a^{2}+2 b^{2}\right ) \tanh \left (\frac {d x}{2}+\frac {c}{2}\right )}{2 \left (a^{4}+2 a^{2} b^{2}+b^{4}\right ) a}+\frac {b \left (4 a^{2}+b^{2}\right )}{2 a^{4}+4 a^{2} b^{2}+2 b^{4}}\right )}{\left (\tanh \left (\frac {d x}{2}+\frac {c}{2}\right )^{2} a -2 b \tanh \left (\frac {d x}{2}+\frac {c}{2}\right )-a \right )^{2}}+\frac {\left (2 a^{2}-b^{2}\right ) \operatorname {arctanh}\left (\frac {2 a \tanh \left (\frac {d x}{2}+\frac {c}{2}\right )-2 b}{2 \sqrt {a^{2}+b^{2}}}\right )}{\left (a^{4}+2 a^{2} b^{2}+b^{4}\right ) \sqrt {a^{2}+b^{2}}}}{d}\) \(280\)
risch \(\frac {2 a^{2} b \,{\mathrm e}^{3 d x +3 c}-b^{3} {\mathrm e}^{3 d x +3 c}+6 a^{3} {\mathrm e}^{2 d x +2 c}-3 a \,b^{2} {\mathrm e}^{2 d x +2 c}-10 a^{2} b \,{\mathrm e}^{d x +c}-{\mathrm e}^{d x +c} b^{3}+3 a \,b^{2}}{d \left (a^{2}+b^{2}\right )^{2} \left (b \,{\mathrm e}^{2 d x +2 c}+2 a \,{\mathrm e}^{d x +c}-b \right )^{2}}+\frac {\ln \left ({\mathrm e}^{d x +c}+\frac {\left (a^{2}+b^{2}\right )^{\frac {5}{2}} a -a^{6}-3 a^{4} b^{2}-3 a^{2} b^{4}-b^{6}}{b \left (a^{2}+b^{2}\right )^{\frac {5}{2}}}\right ) a^{2}}{\left (a^{2}+b^{2}\right )^{\frac {5}{2}} d}-\frac {\ln \left ({\mathrm e}^{d x +c}+\frac {\left (a^{2}+b^{2}\right )^{\frac {5}{2}} a -a^{6}-3 a^{4} b^{2}-3 a^{2} b^{4}-b^{6}}{b \left (a^{2}+b^{2}\right )^{\frac {5}{2}}}\right ) b^{2}}{2 \left (a^{2}+b^{2}\right )^{\frac {5}{2}} d}-\frac {\ln \left ({\mathrm e}^{d x +c}+\frac {\left (a^{2}+b^{2}\right )^{\frac {5}{2}} a +a^{6}+3 a^{4} b^{2}+3 a^{2} b^{4}+b^{6}}{b \left (a^{2}+b^{2}\right )^{\frac {5}{2}}}\right ) a^{2}}{\left (a^{2}+b^{2}\right )^{\frac {5}{2}} d}+\frac {\ln \left ({\mathrm e}^{d x +c}+\frac {\left (a^{2}+b^{2}\right )^{\frac {5}{2}} a +a^{6}+3 a^{4} b^{2}+3 a^{2} b^{4}+b^{6}}{b \left (a^{2}+b^{2}\right )^{\frac {5}{2}}}\right ) b^{2}}{2 \left (a^{2}+b^{2}\right )^{\frac {5}{2}} d}\) \(424\)

Input:

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

Output:

1/d*(-2*(-1/2*b^2*(5*a^2+2*b^2)/a/(a^4+2*a^2*b^2+b^4)*tanh(1/2*d*x+1/2*c)^ 
3-1/2*b*(4*a^4-7*a^2*b^2-2*b^4)/(a^4+2*a^2*b^2+b^4)/a^2*tanh(1/2*d*x+1/2*c 
)^2+1/2*b^2*(11*a^2+2*b^2)/(a^4+2*a^2*b^2+b^4)/a*tanh(1/2*d*x+1/2*c)+1/2*b 
*(4*a^2+b^2)/(a^4+2*a^2*b^2+b^4))/(tanh(1/2*d*x+1/2*c)^2*a-2*b*tanh(1/2*d* 
x+1/2*c)-a)^2+(2*a^2-b^2)/(a^4+2*a^2*b^2+b^4)/(a^2+b^2)^(1/2)*arctanh(1/2* 
(2*a*tanh(1/2*d*x+1/2*c)-2*b)/(a^2+b^2)^(1/2)))
 

Fricas [B] (verification not implemented)

Leaf count of result is larger than twice the leaf count of optimal. 1347 vs. \(2 (120) = 240\).

Time = 0.11 (sec) , antiderivative size = 1347, normalized size of antiderivative = 10.61 \[ \int \frac {1}{(a+b \sinh (c+d x))^3} \, dx=\text {Too large to display} \] Input:

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

Output:

1/2*(6*a^3*b^2 + 6*a*b^4 + 2*(2*a^4*b + a^2*b^3 - b^5)*cosh(d*x + c)^3 + 2 
*(2*a^4*b + a^2*b^3 - b^5)*sinh(d*x + c)^3 + 6*(2*a^5 + a^3*b^2 - a*b^4)*c 
osh(d*x + c)^2 + 6*(2*a^5 + a^3*b^2 - a*b^4 + (2*a^4*b + a^2*b^3 - b^5)*co 
sh(d*x + c))*sinh(d*x + c)^2 - ((2*a^2*b^2 - b^4)*cosh(d*x + c)^4 + (2*a^2 
*b^2 - b^4)*sinh(d*x + c)^4 + 2*a^2*b^2 - b^4 + 4*(2*a^3*b - a*b^3)*cosh(d 
*x + c)^3 + 4*(2*a^3*b - a*b^3 + (2*a^2*b^2 - b^4)*cosh(d*x + c))*sinh(d*x 
 + c)^3 + 2*(4*a^4 - 4*a^2*b^2 + b^4)*cosh(d*x + c)^2 + 2*(4*a^4 - 4*a^2*b 
^2 + b^4 + 3*(2*a^2*b^2 - b^4)*cosh(d*x + c)^2 + 6*(2*a^3*b - a*b^3)*cosh( 
d*x + c))*sinh(d*x + c)^2 - 4*(2*a^3*b - a*b^3)*cosh(d*x + c) - 4*(2*a^3*b 
 - a*b^3 - (2*a^2*b^2 - b^4)*cosh(d*x + c)^3 - 3*(2*a^3*b - a*b^3)*cosh(d* 
x + c)^2 - (4*a^4 - 4*a^2*b^2 + b^4)*cosh(d*x + c))*sinh(d*x + c))*sqrt(a^ 
2 + b^2)*log((b^2*cosh(d*x + c)^2 + b^2*sinh(d*x + c)^2 + 2*a*b*cosh(d*x + 
 c) + 2*a^2 + b^2 + 2*(b^2*cosh(d*x + c) + a*b)*sinh(d*x + c) + 2*sqrt(a^2 
 + b^2)*(b*cosh(d*x + c) + b*sinh(d*x + c) + a))/(b*cosh(d*x + c)^2 + b*si 
nh(d*x + c)^2 + 2*a*cosh(d*x + c) + 2*(b*cosh(d*x + c) + a)*sinh(d*x + c) 
- b)) - 2*(10*a^4*b + 11*a^2*b^3 + b^5)*cosh(d*x + c) - 2*(10*a^4*b + 11*a 
^2*b^3 + b^5 - 3*(2*a^4*b + a^2*b^3 - b^5)*cosh(d*x + c)^2 - 6*(2*a^5 + a^ 
3*b^2 - a*b^4)*cosh(d*x + c))*sinh(d*x + c))/((a^6*b^2 + 3*a^4*b^4 + 3*a^2 
*b^6 + b^8)*d*cosh(d*x + c)^4 + (a^6*b^2 + 3*a^4*b^4 + 3*a^2*b^6 + b^8)*d* 
sinh(d*x + c)^4 + 4*(a^7*b + 3*a^5*b^3 + 3*a^3*b^5 + a*b^7)*d*cosh(d*x ...
 

Sympy [F(-1)]

Timed out. \[ \int \frac {1}{(a+b \sinh (c+d x))^3} \, dx=\text {Timed out} \] Input:

integrate(1/(a+b*sinh(d*x+c))**3,x)
 

Output:

Timed out
 

Maxima [B] (verification not implemented)

Leaf count of result is larger than twice the leaf count of optimal. 315 vs. \(2 (120) = 240\).

Time = 0.15 (sec) , antiderivative size = 315, normalized size of antiderivative = 2.48 \[ \int \frac {1}{(a+b \sinh (c+d x))^3} \, dx=\frac {{\left (2 \, a^{2} - b^{2}\right )} \log \left (\frac {b e^{\left (-d x - c\right )} - a - \sqrt {a^{2} + b^{2}}}{b e^{\left (-d x - c\right )} - a + \sqrt {a^{2} + b^{2}}}\right )}{2 \, {\left (a^{4} + 2 \, a^{2} b^{2} + b^{4}\right )} \sqrt {a^{2} + b^{2}} d} - \frac {3 \, a b^{2} + {\left (10 \, a^{2} b + b^{3}\right )} e^{\left (-d x - c\right )} + 3 \, {\left (2 \, a^{3} - a b^{2}\right )} e^{\left (-2 \, d x - 2 \, c\right )} - {\left (2 \, a^{2} b - b^{3}\right )} e^{\left (-3 \, d x - 3 \, c\right )}}{{\left (a^{4} b^{2} + 2 \, a^{2} b^{4} + b^{6} + 4 \, {\left (a^{5} b + 2 \, a^{3} b^{3} + a b^{5}\right )} e^{\left (-d x - c\right )} + 2 \, {\left (2 \, a^{6} + 3 \, a^{4} b^{2} - b^{6}\right )} e^{\left (-2 \, d x - 2 \, c\right )} - 4 \, {\left (a^{5} b + 2 \, a^{3} b^{3} + a b^{5}\right )} e^{\left (-3 \, d x - 3 \, c\right )} + {\left (a^{4} b^{2} + 2 \, a^{2} b^{4} + b^{6}\right )} e^{\left (-4 \, d x - 4 \, c\right )}\right )} d} \] Input:

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

Output:

1/2*(2*a^2 - b^2)*log((b*e^(-d*x - c) - a - sqrt(a^2 + b^2))/(b*e^(-d*x - 
c) - a + sqrt(a^2 + b^2)))/((a^4 + 2*a^2*b^2 + b^4)*sqrt(a^2 + b^2)*d) - ( 
3*a*b^2 + (10*a^2*b + b^3)*e^(-d*x - c) + 3*(2*a^3 - a*b^2)*e^(-2*d*x - 2* 
c) - (2*a^2*b - b^3)*e^(-3*d*x - 3*c))/((a^4*b^2 + 2*a^2*b^4 + b^6 + 4*(a^ 
5*b + 2*a^3*b^3 + a*b^5)*e^(-d*x - c) + 2*(2*a^6 + 3*a^4*b^2 - b^6)*e^(-2* 
d*x - 2*c) - 4*(a^5*b + 2*a^3*b^3 + a*b^5)*e^(-3*d*x - 3*c) + (a^4*b^2 + 2 
*a^2*b^4 + b^6)*e^(-4*d*x - 4*c))*d)
 

Giac [A] (verification not implemented)

Time = 0.13 (sec) , antiderivative size = 231, normalized size of antiderivative = 1.82 \[ \int \frac {1}{(a+b \sinh (c+d x))^3} \, dx=\frac {\frac {{\left (2 \, a^{2} - b^{2}\right )} \log \left (\frac {{\left | 2 \, b e^{\left (d x + c\right )} + 2 \, a - 2 \, \sqrt {a^{2} + b^{2}} \right |}}{{\left | 2 \, b e^{\left (d x + c\right )} + 2 \, a + 2 \, \sqrt {a^{2} + b^{2}} \right |}}\right )}{{\left (a^{4} + 2 \, a^{2} b^{2} + b^{4}\right )} \sqrt {a^{2} + b^{2}}} + \frac {2 \, {\left (2 \, a^{2} b e^{\left (3 \, d x + 3 \, c\right )} - b^{3} e^{\left (3 \, d x + 3 \, c\right )} + 6 \, a^{3} e^{\left (2 \, d x + 2 \, c\right )} - 3 \, a b^{2} e^{\left (2 \, d x + 2 \, c\right )} - 10 \, a^{2} b e^{\left (d x + c\right )} - b^{3} e^{\left (d x + c\right )} + 3 \, a b^{2}\right )}}{{\left (a^{4} + 2 \, a^{2} b^{2} + b^{4}\right )} {\left (b e^{\left (2 \, d x + 2 \, c\right )} + 2 \, a e^{\left (d x + c\right )} - b\right )}^{2}}}{2 \, d} \] Input:

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

Output:

1/2*((2*a^2 - b^2)*log(abs(2*b*e^(d*x + c) + 2*a - 2*sqrt(a^2 + b^2))/abs( 
2*b*e^(d*x + c) + 2*a + 2*sqrt(a^2 + b^2)))/((a^4 + 2*a^2*b^2 + b^4)*sqrt( 
a^2 + b^2)) + 2*(2*a^2*b*e^(3*d*x + 3*c) - b^3*e^(3*d*x + 3*c) + 6*a^3*e^( 
2*d*x + 2*c) - 3*a*b^2*e^(2*d*x + 2*c) - 10*a^2*b*e^(d*x + c) - b^3*e^(d*x 
 + c) + 3*a*b^2)/((a^4 + 2*a^2*b^2 + b^4)*(b*e^(2*d*x + 2*c) + 2*a*e^(d*x 
+ c) - b)^2))/d
 

Mupad [F(-1)]

Timed out. \[ \int \frac {1}{(a+b \sinh (c+d x))^3} \, dx=\int \frac {1}{{\left (a+b\,\mathrm {sinh}\left (c+d\,x\right )\right )}^3} \,d x \] Input:

int(1/(a + b*sinh(c + d*x))^3,x)
 

Output:

int(1/(a + b*sinh(c + d*x))^3, x)
 

Reduce [B] (verification not implemented)

Time = 0.16 (sec) , antiderivative size = 1033, normalized size of antiderivative = 8.13 \[ \int \frac {1}{(a+b \sinh (c+d x))^3} \, dx =\text {Too large to display} \] Input:

int(1/(a+b*sinh(d*x+c))^3,x)
 

Output:

(8*e**(4*c + 4*d*x)*sqrt(a**2 + b**2)*atan((e**(c + d*x)*b*i + a*i)/sqrt(a 
**2 + b**2))*a**3*b**2*i - 4*e**(4*c + 4*d*x)*sqrt(a**2 + b**2)*atan((e**( 
c + d*x)*b*i + a*i)/sqrt(a**2 + b**2))*a*b**4*i + 32*e**(3*c + 3*d*x)*sqrt 
(a**2 + b**2)*atan((e**(c + d*x)*b*i + a*i)/sqrt(a**2 + b**2))*a**4*b*i - 
16*e**(3*c + 3*d*x)*sqrt(a**2 + b**2)*atan((e**(c + d*x)*b*i + a*i)/sqrt(a 
**2 + b**2))*a**2*b**3*i + 32*e**(2*c + 2*d*x)*sqrt(a**2 + b**2)*atan((e** 
(c + d*x)*b*i + a*i)/sqrt(a**2 + b**2))*a**5*i - 32*e**(2*c + 2*d*x)*sqrt( 
a**2 + b**2)*atan((e**(c + d*x)*b*i + a*i)/sqrt(a**2 + b**2))*a**3*b**2*i 
+ 8*e**(2*c + 2*d*x)*sqrt(a**2 + b**2)*atan((e**(c + d*x)*b*i + a*i)/sqrt( 
a**2 + b**2))*a*b**4*i - 32*e**(c + d*x)*sqrt(a**2 + b**2)*atan((e**(c + d 
*x)*b*i + a*i)/sqrt(a**2 + b**2))*a**4*b*i + 16*e**(c + d*x)*sqrt(a**2 + b 
**2)*atan((e**(c + d*x)*b*i + a*i)/sqrt(a**2 + b**2))*a**2*b**3*i + 8*sqrt 
(a**2 + b**2)*atan((e**(c + d*x)*b*i + a*i)/sqrt(a**2 + b**2))*a**3*b**2*i 
 - 4*sqrt(a**2 + b**2)*atan((e**(c + d*x)*b*i + a*i)/sqrt(a**2 + b**2))*a* 
b**4*i - 2*e**(4*c + 4*d*x)*a**4*b**2 - e**(4*c + 4*d*x)*a**2*b**4 + e**(4 
*c + 4*d*x)*b**6 + 16*e**(2*c + 2*d*x)*a**6 + 12*e**(2*c + 2*d*x)*a**4*b** 
2 - 6*e**(2*c + 2*d*x)*a**2*b**4 - 2*e**(2*c + 2*d*x)*b**6 - 32*e**(c + d* 
x)*a**5*b - 40*e**(c + d*x)*a**3*b**3 - 8*e**(c + d*x)*a*b**5 + 10*a**4*b* 
*2 + 11*a**2*b**4 + b**6)/(4*a*d*(e**(4*c + 4*d*x)*a**6*b**2 + 3*e**(4*c + 
 4*d*x)*a**4*b**4 + 3*e**(4*c + 4*d*x)*a**2*b**6 + e**(4*c + 4*d*x)*b**...