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

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 [B] (verification not implemented)
Mupad [F(-1)]
Reduce [B] (verification not implemented)

Optimal result

Integrand size = 12, antiderivative size = 174 \[ \int \frac {1}{(a+b \sinh (c+d x))^4} \, dx=-\frac {a \left (2 a^2-3 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 )^{7/2} d}-\frac {b \cosh (c+d x)}{3 \left (a^2+b^2\right ) d (a+b \sinh (c+d x))^3}-\frac {5 a b \cosh (c+d x)}{6 \left (a^2+b^2\right )^2 d (a+b \sinh (c+d x))^2}-\frac {b \left (11 a^2-4 b^2\right ) \cosh (c+d x)}{6 \left (a^2+b^2\right )^3 d (a+b \sinh (c+d x))} \] Output:

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

Mathematica [A] (verified)

Time = 0.50 (sec) , antiderivative size = 159, normalized size of antiderivative = 0.91 \[ \int \frac {1}{(a+b \sinh (c+d x))^4} \, dx=\frac {\frac {6 a \left (2 a^2-3 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 (-18 a^4-5 a^2 b^2-2 b^4+3 a b \left (-9 a^2+b^2\right ) \sinh (c+d x)+\left (-11 a^2 b^2+4 b^4\right ) \sinh ^2(c+d x)\right )}{(a+b \sinh (c+d x))^3}}{6 \left (a^2+b^2\right )^3 d} \] Input:

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

Output:

((6*a*(2*a^2 - 3*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]*(-18*a^4 - 5*a^2*b^2 - 2*b^4 + 3*a*b*( 
-9*a^2 + b^2)*Sinh[c + d*x] + (-11*a^2*b^2 + 4*b^4)*Sinh[c + d*x]^2))/(a + 
 b*Sinh[c + d*x])^3)/(6*(a^2 + b^2)^3*d)
 

Rubi [A] (warning: unable to verify)

Time = 0.87 (sec) , antiderivative size = 198, normalized size of antiderivative = 1.14, number of steps used = 14, number of rules used = 13, \(\frac {\text {number of rules}}{\text {integrand size}}\) = 1.083, Rules used = {3042, 3143, 25, 3042, 3233, 25, 3042, 3233, 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))^4} \, dx\)

\(\Big \downarrow \) 3042

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

\(\Big \downarrow \) 3143

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

\(\Big \downarrow \) 25

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

\(\Big \downarrow \) 3042

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

\(\Big \downarrow \) 3233

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

\(\Big \downarrow \) 25

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

\(\Big \downarrow \) 3042

\(\displaystyle -\frac {b \cosh (c+d x)}{3 d \left (a^2+b^2\right ) (a+b \sinh (c+d x))^3}+\frac {-\frac {5 a b \cosh (c+d x)}{2 d \left (a^2+b^2\right ) (a+b \sinh (c+d x))^2}+\frac {\int \frac {2 \left (3 a^2-2 b^2\right )+5 i a 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 )}}{3 \left (a^2+b^2\right )}\)

\(\Big \downarrow \) 3233

\(\displaystyle \frac {\frac {-\frac {\int -\frac {3 a \left (2 a^2-3 b^2\right )}{a+b \sinh (c+d x)}dx}{a^2+b^2}-\frac {b \left (11 a^2-4 b^2\right ) \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 {5 a b \cosh (c+d x)}{2 d \left (a^2+b^2\right ) (a+b \sinh (c+d x))^2}}{3 \left (a^2+b^2\right )}-\frac {b \cosh (c+d x)}{3 d \left (a^2+b^2\right ) (a+b \sinh (c+d x))^3}\)

\(\Big \downarrow \) 27

\(\displaystyle \frac {\frac {\frac {3 a \left (2 a^2-3 b^2\right ) \int \frac {1}{a+b \sinh (c+d x)}dx}{a^2+b^2}-\frac {b \left (11 a^2-4 b^2\right ) \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 {5 a b \cosh (c+d x)}{2 d \left (a^2+b^2\right ) (a+b \sinh (c+d x))^2}}{3 \left (a^2+b^2\right )}-\frac {b \cosh (c+d x)}{3 d \left (a^2+b^2\right ) (a+b \sinh (c+d x))^3}\)

\(\Big \downarrow \) 3042

\(\displaystyle -\frac {b \cosh (c+d x)}{3 d \left (a^2+b^2\right ) (a+b \sinh (c+d x))^3}+\frac {-\frac {5 a b \cosh (c+d x)}{2 d \left (a^2+b^2\right ) (a+b \sinh (c+d x))^2}+\frac {-\frac {b \left (11 a^2-4 b^2\right ) \cosh (c+d x)}{d \left (a^2+b^2\right ) (a+b \sinh (c+d x))}+\frac {3 a \left (2 a^2-3 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 )}}{3 \left (a^2+b^2\right )}\)

\(\Big \downarrow \) 3139

\(\displaystyle -\frac {b \cosh (c+d x)}{3 d \left (a^2+b^2\right ) (a+b \sinh (c+d x))^3}+\frac {-\frac {5 a b \cosh (c+d x)}{2 d \left (a^2+b^2\right ) (a+b \sinh (c+d x))^2}+\frac {-\frac {b \left (11 a^2-4 b^2\right ) \cosh (c+d x)}{d \left (a^2+b^2\right ) (a+b \sinh (c+d x))}-\frac {6 i a \left (2 a^2-3 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 )}}{3 \left (a^2+b^2\right )}\)

\(\Big \downarrow \) 1083

\(\displaystyle -\frac {b \cosh (c+d x)}{3 d \left (a^2+b^2\right ) (a+b \sinh (c+d x))^3}+\frac {-\frac {5 a b \cosh (c+d x)}{2 d \left (a^2+b^2\right ) (a+b \sinh (c+d x))^2}+\frac {-\frac {b \left (11 a^2-4 b^2\right ) \cosh (c+d x)}{d \left (a^2+b^2\right ) (a+b \sinh (c+d x))}+\frac {12 i a \left (2 a^2-3 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 )}}{3 \left (a^2+b^2\right )}\)

\(\Big \downarrow \) 217

\(\displaystyle \frac {\frac {\frac {6 a \left (2 a^2-3 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 {b \left (11 a^2-4 b^2\right ) \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 {5 a b \cosh (c+d x)}{2 d \left (a^2+b^2\right ) (a+b \sinh (c+d x))^2}}{3 \left (a^2+b^2\right )}-\frac {b \cosh (c+d x)}{3 d \left (a^2+b^2\right ) (a+b \sinh (c+d x))^3}\)

Input:

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

Output:

-1/3*(b*Cosh[c + d*x])/((a^2 + b^2)*d*(a + b*Sinh[c + d*x])^3) + ((-5*a*b* 
Cosh[c + d*x])/(2*(a^2 + b^2)*d*(a + b*Sinh[c + d*x])^2) + ((6*a*(2*a^2 - 
3*b^2)*ArcTanh[Tanh[(c + d*x)/2]/(2*Sqrt[a^2 + b^2])])/((a^2 + b^2)^(3/2)* 
d) - (b*(11*a^2 - 4*b^2)*Cosh[c + d*x])/((a^2 + b^2)*d*(a + b*Sinh[c + d*x 
])))/(2*(a^2 + b^2)))/(3*(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. \(493\) vs. \(2(163)=326\).

Time = 0.77 (sec) , antiderivative size = 494, normalized size of antiderivative = 2.84

method result size
derivativedivides \(\frac {-\frac {2 \left (-\frac {b^{2} \left (9 a^{4}+6 a^{2} b^{2}+2 b^{4}\right ) \tanh \left (\frac {d x}{2}+\frac {c}{2}\right )^{5}}{2 a \left (a^{6}+3 a^{4} b^{2}+3 a^{2} b^{4}+b^{6}\right )}-\frac {b \left (6 a^{6}-27 a^{4} b^{2}-12 a^{2} b^{4}-4 b^{6}\right ) \tanh \left (\frac {d x}{2}+\frac {c}{2}\right )^{4}}{2 a^{2} \left (a^{6}+3 a^{4} b^{2}+3 a^{2} b^{4}+b^{6}\right )}+\frac {b^{2} \left (54 a^{6}-21 a^{4} b^{2}-4 a^{2} b^{4}-4 b^{6}\right ) \tanh \left (\frac {d x}{2}+\frac {c}{2}\right )^{3}}{3 a^{3} \left (a^{6}+3 a^{4} b^{2}+3 a^{2} b^{4}+b^{6}\right )}+\frac {b \left (6 a^{6}-20 a^{4} b^{2}-3 a^{2} b^{4}-2 b^{6}\right ) \tanh \left (\frac {d x}{2}+\frac {c}{2}\right )^{2}}{a^{2} \left (a^{6}+3 a^{4} b^{2}+3 a^{2} b^{4}+b^{6}\right )}-\frac {b^{2} \left (27 a^{4}+4 a^{2} b^{2}+2 b^{4}\right ) \tanh \left (\frac {d x}{2}+\frac {c}{2}\right )}{2 a \left (a^{6}+3 a^{4} b^{2}+3 a^{2} b^{4}+b^{6}\right )}-\frac {b \left (18 a^{4}+5 a^{2} b^{2}+2 b^{4}\right )}{6 \left (a^{6}+3 a^{4} b^{2}+3 a^{2} b^{4}+b^{6}\right )}\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 )^{3}}+\frac {a \left (2 a^{2}-3 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^{6}+3 a^{4} b^{2}+3 a^{2} b^{4}+b^{6}\right ) \sqrt {a^{2}+b^{2}}}}{d}\) \(494\)
default \(\frac {-\frac {2 \left (-\frac {b^{2} \left (9 a^{4}+6 a^{2} b^{2}+2 b^{4}\right ) \tanh \left (\frac {d x}{2}+\frac {c}{2}\right )^{5}}{2 a \left (a^{6}+3 a^{4} b^{2}+3 a^{2} b^{4}+b^{6}\right )}-\frac {b \left (6 a^{6}-27 a^{4} b^{2}-12 a^{2} b^{4}-4 b^{6}\right ) \tanh \left (\frac {d x}{2}+\frac {c}{2}\right )^{4}}{2 a^{2} \left (a^{6}+3 a^{4} b^{2}+3 a^{2} b^{4}+b^{6}\right )}+\frac {b^{2} \left (54 a^{6}-21 a^{4} b^{2}-4 a^{2} b^{4}-4 b^{6}\right ) \tanh \left (\frac {d x}{2}+\frac {c}{2}\right )^{3}}{3 a^{3} \left (a^{6}+3 a^{4} b^{2}+3 a^{2} b^{4}+b^{6}\right )}+\frac {b \left (6 a^{6}-20 a^{4} b^{2}-3 a^{2} b^{4}-2 b^{6}\right ) \tanh \left (\frac {d x}{2}+\frac {c}{2}\right )^{2}}{a^{2} \left (a^{6}+3 a^{4} b^{2}+3 a^{2} b^{4}+b^{6}\right )}-\frac {b^{2} \left (27 a^{4}+4 a^{2} b^{2}+2 b^{4}\right ) \tanh \left (\frac {d x}{2}+\frac {c}{2}\right )}{2 a \left (a^{6}+3 a^{4} b^{2}+3 a^{2} b^{4}+b^{6}\right )}-\frac {b \left (18 a^{4}+5 a^{2} b^{2}+2 b^{4}\right )}{6 \left (a^{6}+3 a^{4} b^{2}+3 a^{2} b^{4}+b^{6}\right )}\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 )^{3}}+\frac {a \left (2 a^{2}-3 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^{6}+3 a^{4} b^{2}+3 a^{2} b^{4}+b^{6}\right ) \sqrt {a^{2}+b^{2}}}}{d}\) \(494\)
risch \(\frac {6 a^{3} b^{2} {\mathrm e}^{5 d x +5 c}-9 a \,b^{4} {\mathrm e}^{5 d x +5 c}+30 a^{4} b \,{\mathrm e}^{4 d x +4 c}-45 a^{2} b^{3} {\mathrm e}^{4 d x +4 c}+44 a^{5} {\mathrm e}^{3 d x +3 c}-82 a^{3} b^{2} {\mathrm e}^{3 d x +3 c}+24 a \,b^{4} {\mathrm e}^{3 d x +3 c}-102 a^{4} b \,{\mathrm e}^{2 d x +2 c}+36 a^{2} b^{3} {\mathrm e}^{2 d x +2 c}-12 b^{5} {\mathrm e}^{2 d x +2 c}+60 a^{3} b^{2} {\mathrm e}^{d x +c}-15 a \,b^{4} {\mathrm e}^{d x +c}-11 a^{2} b^{3}+4 b^{5}}{3 d \left (a^{2}+b^{2}\right )^{3} \left (b \,{\mathrm e}^{2 d x +2 c}+2 a \,{\mathrm e}^{d x +c}-b \right )^{3}}+\frac {a^{3} \ln \left ({\mathrm e}^{d x +c}+\frac {\left (a^{2}+b^{2}\right )^{\frac {7}{2}} a -a^{8}-4 a^{6} b^{2}-6 b^{4} a^{4}-4 a^{2} b^{6}-b^{8}}{b \left (a^{2}+b^{2}\right )^{\frac {7}{2}}}\right )}{\left (a^{2}+b^{2}\right )^{\frac {7}{2}} d}-\frac {3 a \ln \left ({\mathrm e}^{d x +c}+\frac {\left (a^{2}+b^{2}\right )^{\frac {7}{2}} a -a^{8}-4 a^{6} b^{2}-6 b^{4} a^{4}-4 a^{2} b^{6}-b^{8}}{b \left (a^{2}+b^{2}\right )^{\frac {7}{2}}}\right ) b^{2}}{2 \left (a^{2}+b^{2}\right )^{\frac {7}{2}} d}-\frac {a^{3} \ln \left ({\mathrm e}^{d x +c}+\frac {\left (a^{2}+b^{2}\right )^{\frac {7}{2}} a +a^{8}+4 a^{6} b^{2}+6 b^{4} a^{4}+4 a^{2} b^{6}+b^{8}}{b \left (a^{2}+b^{2}\right )^{\frac {7}{2}}}\right )}{\left (a^{2}+b^{2}\right )^{\frac {7}{2}} d}+\frac {3 a \ln \left ({\mathrm e}^{d x +c}+\frac {\left (a^{2}+b^{2}\right )^{\frac {7}{2}} a +a^{8}+4 a^{6} b^{2}+6 b^{4} a^{4}+4 a^{2} b^{6}+b^{8}}{b \left (a^{2}+b^{2}\right )^{\frac {7}{2}}}\right ) b^{2}}{2 \left (a^{2}+b^{2}\right )^{\frac {7}{2}} d}\) \(567\)

Input:

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

Output:

1/d*(-2*(-1/2*b^2*(9*a^4+6*a^2*b^2+2*b^4)/a/(a^6+3*a^4*b^2+3*a^2*b^4+b^6)* 
tanh(1/2*d*x+1/2*c)^5-1/2*b*(6*a^6-27*a^4*b^2-12*a^2*b^4-4*b^6)/a^2/(a^6+3 
*a^4*b^2+3*a^2*b^4+b^6)*tanh(1/2*d*x+1/2*c)^4+1/3/a^3*b^2*(54*a^6-21*a^4*b 
^2-4*a^2*b^4-4*b^6)/(a^6+3*a^4*b^2+3*a^2*b^4+b^6)*tanh(1/2*d*x+1/2*c)^3+1/ 
a^2*b*(6*a^6-20*a^4*b^2-3*a^2*b^4-2*b^6)/(a^6+3*a^4*b^2+3*a^2*b^4+b^6)*tan 
h(1/2*d*x+1/2*c)^2-1/2/a*b^2*(27*a^4+4*a^2*b^2+2*b^4)/(a^6+3*a^4*b^2+3*a^2 
*b^4+b^6)*tanh(1/2*d*x+1/2*c)-1/6*b*(18*a^4+5*a^2*b^2+2*b^4)/(a^6+3*a^4*b^ 
2+3*a^2*b^4+b^6))/(tanh(1/2*d*x+1/2*c)^2*a-2*b*tanh(1/2*d*x+1/2*c)-a)^3+a* 
(2*a^2-3*b^2)/(a^6+3*a^4*b^2+3*a^2*b^4+b^6)/(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. 2934 vs. \(2 (165) = 330\).

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

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

Output:

-1/6*(22*a^4*b^3 + 14*a^2*b^5 - 8*b^7 - 6*(2*a^5*b^2 - a^3*b^4 - 3*a*b^6)* 
cosh(d*x + c)^5 - 6*(2*a^5*b^2 - a^3*b^4 - 3*a*b^6)*sinh(d*x + c)^5 - 30*( 
2*a^6*b - a^4*b^3 - 3*a^2*b^5)*cosh(d*x + c)^4 - 30*(2*a^6*b - a^4*b^3 - 3 
*a^2*b^5 + (2*a^5*b^2 - a^3*b^4 - 3*a*b^6)*cosh(d*x + c))*sinh(d*x + c)^4 
- 4*(22*a^7 - 19*a^5*b^2 - 29*a^3*b^4 + 12*a*b^6)*cosh(d*x + c)^3 - 4*(22* 
a^7 - 19*a^5*b^2 - 29*a^3*b^4 + 12*a*b^6 + 15*(2*a^5*b^2 - a^3*b^4 - 3*a*b 
^6)*cosh(d*x + c)^2 + 30*(2*a^6*b - a^4*b^3 - 3*a^2*b^5)*cosh(d*x + c))*si 
nh(d*x + c)^3 + 12*(17*a^6*b + 11*a^4*b^3 - 4*a^2*b^5 + 2*b^7)*cosh(d*x + 
c)^2 + 12*(17*a^6*b + 11*a^4*b^3 - 4*a^2*b^5 + 2*b^7 - 5*(2*a^5*b^2 - a^3* 
b^4 - 3*a*b^6)*cosh(d*x + c)^3 - 15*(2*a^6*b - a^4*b^3 - 3*a^2*b^5)*cosh(d 
*x + c)^2 - (22*a^7 - 19*a^5*b^2 - 29*a^3*b^4 + 12*a*b^6)*cosh(d*x + c))*s 
inh(d*x + c)^2 + 3*((2*a^3*b^3 - 3*a*b^5)*cosh(d*x + c)^6 + (2*a^3*b^3 - 3 
*a*b^5)*sinh(d*x + c)^6 - 2*a^3*b^3 + 3*a*b^5 + 6*(2*a^4*b^2 - 3*a^2*b^4)* 
cosh(d*x + c)^5 + 6*(2*a^4*b^2 - 3*a^2*b^4 + (2*a^3*b^3 - 3*a*b^5)*cosh(d* 
x + c))*sinh(d*x + c)^5 + 3*(8*a^5*b - 14*a^3*b^3 + 3*a*b^5)*cosh(d*x + c) 
^4 + 3*(8*a^5*b - 14*a^3*b^3 + 3*a*b^5 + 5*(2*a^3*b^3 - 3*a*b^5)*cosh(d*x 
+ c)^2 + 10*(2*a^4*b^2 - 3*a^2*b^4)*cosh(d*x + c))*sinh(d*x + c)^4 + 4*(4* 
a^6 - 12*a^4*b^2 + 9*a^2*b^4)*cosh(d*x + c)^3 + 4*(4*a^6 - 12*a^4*b^2 + 9* 
a^2*b^4 + 5*(2*a^3*b^3 - 3*a*b^5)*cosh(d*x + c)^3 + 15*(2*a^4*b^2 - 3*a^2* 
b^4)*cosh(d*x + c)^2 + 3*(8*a^5*b - 14*a^3*b^3 + 3*a*b^5)*cosh(d*x + c)...
 

Sympy [F(-1)]

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

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

Output:

Timed out
 

Maxima [B] (verification not implemented)

Leaf count of result is larger than twice the leaf count of optimal. 551 vs. \(2 (165) = 330\).

Time = 0.15 (sec) , antiderivative size = 551, normalized size of antiderivative = 3.17 \[ \int \frac {1}{(a+b \sinh (c+d x))^4} \, dx=\frac {{\left (2 \, a^{2} - 3 \, b^{2}\right )} a \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^{6} + 3 \, a^{4} b^{2} + 3 \, a^{2} b^{4} + b^{6}\right )} \sqrt {a^{2} + b^{2}} d} - \frac {11 \, a^{2} b^{3} - 4 \, b^{5} + 15 \, {\left (4 \, a^{3} b^{2} - a b^{4}\right )} e^{\left (-d x - c\right )} + 6 \, {\left (17 \, a^{4} b - 6 \, a^{2} b^{3} + 2 \, b^{5}\right )} e^{\left (-2 \, d x - 2 \, c\right )} + 2 \, {\left (22 \, a^{5} - 41 \, a^{3} b^{2} + 12 \, a b^{4}\right )} e^{\left (-3 \, d x - 3 \, c\right )} - 15 \, {\left (2 \, a^{4} b - 3 \, a^{2} b^{3}\right )} e^{\left (-4 \, d x - 4 \, c\right )} + 3 \, {\left (2 \, a^{3} b^{2} - 3 \, a b^{4}\right )} e^{\left (-5 \, d x - 5 \, c\right )}}{3 \, {\left (a^{6} b^{3} + 3 \, a^{4} b^{5} + 3 \, a^{2} b^{7} + b^{9} + 6 \, {\left (a^{7} b^{2} + 3 \, a^{5} b^{4} + 3 \, a^{3} b^{6} + a b^{8}\right )} e^{\left (-d x - c\right )} + 3 \, {\left (4 \, a^{8} b + 11 \, a^{6} b^{3} + 9 \, a^{4} b^{5} + a^{2} b^{7} - b^{9}\right )} e^{\left (-2 \, d x - 2 \, c\right )} + 4 \, {\left (2 \, a^{9} + 3 \, a^{7} b^{2} - 3 \, a^{5} b^{4} - 7 \, a^{3} b^{6} - 3 \, a b^{8}\right )} e^{\left (-3 \, d x - 3 \, c\right )} - 3 \, {\left (4 \, a^{8} b + 11 \, a^{6} b^{3} + 9 \, a^{4} b^{5} + a^{2} b^{7} - b^{9}\right )} e^{\left (-4 \, d x - 4 \, c\right )} + 6 \, {\left (a^{7} b^{2} + 3 \, a^{5} b^{4} + 3 \, a^{3} b^{6} + a b^{8}\right )} e^{\left (-5 \, d x - 5 \, c\right )} - {\left (a^{6} b^{3} + 3 \, a^{4} b^{5} + 3 \, a^{2} b^{7} + b^{9}\right )} e^{\left (-6 \, d x - 6 \, c\right )}\right )} d} \] Input:

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

Output:

1/2*(2*a^2 - 3*b^2)*a*log((b*e^(-d*x - c) - a - sqrt(a^2 + b^2))/(b*e^(-d* 
x - c) - a + sqrt(a^2 + b^2)))/((a^6 + 3*a^4*b^2 + 3*a^2*b^4 + b^6)*sqrt(a 
^2 + b^2)*d) - 1/3*(11*a^2*b^3 - 4*b^5 + 15*(4*a^3*b^2 - a*b^4)*e^(-d*x - 
c) + 6*(17*a^4*b - 6*a^2*b^3 + 2*b^5)*e^(-2*d*x - 2*c) + 2*(22*a^5 - 41*a^ 
3*b^2 + 12*a*b^4)*e^(-3*d*x - 3*c) - 15*(2*a^4*b - 3*a^2*b^3)*e^(-4*d*x - 
4*c) + 3*(2*a^3*b^2 - 3*a*b^4)*e^(-5*d*x - 5*c))/((a^6*b^3 + 3*a^4*b^5 + 3 
*a^2*b^7 + b^9 + 6*(a^7*b^2 + 3*a^5*b^4 + 3*a^3*b^6 + a*b^8)*e^(-d*x - c) 
+ 3*(4*a^8*b + 11*a^6*b^3 + 9*a^4*b^5 + a^2*b^7 - b^9)*e^(-2*d*x - 2*c) + 
4*(2*a^9 + 3*a^7*b^2 - 3*a^5*b^4 - 7*a^3*b^6 - 3*a*b^8)*e^(-3*d*x - 3*c) - 
 3*(4*a^8*b + 11*a^6*b^3 + 9*a^4*b^5 + a^2*b^7 - b^9)*e^(-4*d*x - 4*c) + 6 
*(a^7*b^2 + 3*a^5*b^4 + 3*a^3*b^6 + a*b^8)*e^(-5*d*x - 5*c) - (a^6*b^3 + 3 
*a^4*b^5 + 3*a^2*b^7 + b^9)*e^(-6*d*x - 6*c))*d)
 

Giac [B] (verification not implemented)

Leaf count of result is larger than twice the leaf count of optimal. 357 vs. \(2 (165) = 330\).

Time = 0.14 (sec) , antiderivative size = 357, normalized size of antiderivative = 2.05 \[ \int \frac {1}{(a+b \sinh (c+d x))^4} \, dx=\frac {\frac {3 \, {\left (2 \, a^{3} - 3 \, a 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^{6} + 3 \, a^{4} b^{2} + 3 \, a^{2} b^{4} + b^{6}\right )} \sqrt {a^{2} + b^{2}}} + \frac {2 \, {\left (6 \, a^{3} b^{2} e^{\left (5 \, d x + 5 \, c\right )} - 9 \, a b^{4} e^{\left (5 \, d x + 5 \, c\right )} + 30 \, a^{4} b e^{\left (4 \, d x + 4 \, c\right )} - 45 \, a^{2} b^{3} e^{\left (4 \, d x + 4 \, c\right )} + 44 \, a^{5} e^{\left (3 \, d x + 3 \, c\right )} - 82 \, a^{3} b^{2} e^{\left (3 \, d x + 3 \, c\right )} + 24 \, a b^{4} e^{\left (3 \, d x + 3 \, c\right )} - 102 \, a^{4} b e^{\left (2 \, d x + 2 \, c\right )} + 36 \, a^{2} b^{3} e^{\left (2 \, d x + 2 \, c\right )} - 12 \, b^{5} e^{\left (2 \, d x + 2 \, c\right )} + 60 \, a^{3} b^{2} e^{\left (d x + c\right )} - 15 \, a b^{4} e^{\left (d x + c\right )} - 11 \, a^{2} b^{3} + 4 \, b^{5}\right )}}{{\left (a^{6} + 3 \, a^{4} b^{2} + 3 \, a^{2} b^{4} + b^{6}\right )} {\left (b e^{\left (2 \, d x + 2 \, c\right )} + 2 \, a e^{\left (d x + c\right )} - b\right )}^{3}}}{6 \, d} \] Input:

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

Output:

1/6*(3*(2*a^3 - 3*a*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^6 + 3*a^4*b^2 + 3*a^ 
2*b^4 + b^6)*sqrt(a^2 + b^2)) + 2*(6*a^3*b^2*e^(5*d*x + 5*c) - 9*a*b^4*e^( 
5*d*x + 5*c) + 30*a^4*b*e^(4*d*x + 4*c) - 45*a^2*b^3*e^(4*d*x + 4*c) + 44* 
a^5*e^(3*d*x + 3*c) - 82*a^3*b^2*e^(3*d*x + 3*c) + 24*a*b^4*e^(3*d*x + 3*c 
) - 102*a^4*b*e^(2*d*x + 2*c) + 36*a^2*b^3*e^(2*d*x + 2*c) - 12*b^5*e^(2*d 
*x + 2*c) + 60*a^3*b^2*e^(d*x + c) - 15*a*b^4*e^(d*x + c) - 11*a^2*b^3 + 4 
*b^5)/((a^6 + 3*a^4*b^2 + 3*a^2*b^4 + b^6)*(b*e^(2*d*x + 2*c) + 2*a*e^(d*x 
 + c) - b)^3))/d
 

Mupad [F(-1)]

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

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

Output:

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

Reduce [B] (verification not implemented)

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

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

Output:

(12*e**(6*c + 6*d*x)*sqrt(a**2 + b**2)*atan((e**(c + d*x)*b*i + a*i)/sqrt( 
a**2 + b**2))*a**3*b**3*i - 18*e**(6*c + 6*d*x)*sqrt(a**2 + b**2)*atan((e* 
*(c + d*x)*b*i + a*i)/sqrt(a**2 + b**2))*a*b**5*i + 72*e**(5*c + 5*d*x)*sq 
rt(a**2 + b**2)*atan((e**(c + d*x)*b*i + a*i)/sqrt(a**2 + b**2))*a**4*b**2 
*i - 108*e**(5*c + 5*d*x)*sqrt(a**2 + b**2)*atan((e**(c + d*x)*b*i + a*i)/ 
sqrt(a**2 + b**2))*a**2*b**4*i + 144*e**(4*c + 4*d*x)*sqrt(a**2 + b**2)*at 
an((e**(c + d*x)*b*i + a*i)/sqrt(a**2 + b**2))*a**5*b*i - 252*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**3*i + 54*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**5*i + 96*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**6*i - 288*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**2*i + 216*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**4*i - 144*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*b*i + 252*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**3*i - 54*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**5*i + 72*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**2*i - 108*e**( 
c + d*x)*sqrt(a**2 + b**2)*atan((e**(c + d*x)*b*i + a*i)/sqrt(a**2 + b*...