\(\int \frac {(e+f x) \cosh ^3(c+d x)}{a+b \sinh (c+d x)} \, dx\) [301]

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

Optimal result

Integrand size = 26, antiderivative size = 298 \[ \int \frac {(e+f x) \cosh ^3(c+d x)}{a+b \sinh (c+d x)} \, dx=\frac {f x}{4 b d}-\frac {\left (a^2+b^2\right ) (e+f x)^2}{2 b^3 f}+\frac {a f \cosh (c+d x)}{b^2 d^2}+\frac {\left (a^2+b^2\right ) (e+f x) \log \left (1+\frac {b e^{c+d x}}{a-\sqrt {a^2+b^2}}\right )}{b^3 d}+\frac {\left (a^2+b^2\right ) (e+f x) \log \left (1+\frac {b e^{c+d x}}{a+\sqrt {a^2+b^2}}\right )}{b^3 d}+\frac {\left (a^2+b^2\right ) f \operatorname {PolyLog}\left (2,-\frac {b e^{c+d x}}{a-\sqrt {a^2+b^2}}\right )}{b^3 d^2}+\frac {\left (a^2+b^2\right ) f \operatorname {PolyLog}\left (2,-\frac {b e^{c+d x}}{a+\sqrt {a^2+b^2}}\right )}{b^3 d^2}-\frac {a (e+f x) \sinh (c+d x)}{b^2 d}-\frac {f \cosh (c+d x) \sinh (c+d x)}{4 b d^2}+\frac {(e+f x) \sinh ^2(c+d x)}{2 b d} \] Output:

1/4*f*x/b/d-1/2*(a^2+b^2)*(f*x+e)^2/b^3/f+a*f*cosh(d*x+c)/b^2/d^2+(a^2+b^2 
)*(f*x+e)*ln(1+b*exp(d*x+c)/(a-(a^2+b^2)^(1/2)))/b^3/d+(a^2+b^2)*(f*x+e)*l 
n(1+b*exp(d*x+c)/(a+(a^2+b^2)^(1/2)))/b^3/d+(a^2+b^2)*f*polylog(2,-b*exp(d 
*x+c)/(a-(a^2+b^2)^(1/2)))/b^3/d^2+(a^2+b^2)*f*polylog(2,-b*exp(d*x+c)/(a+ 
(a^2+b^2)^(1/2)))/b^3/d^2-a*(f*x+e)*sinh(d*x+c)/b^2/d-1/4*f*cosh(d*x+c)*si 
nh(d*x+c)/b/d^2+1/2*(f*x+e)*sinh(d*x+c)^2/b/d
 

Mathematica [A] (warning: unable to verify)

Time = 3.29 (sec) , antiderivative size = 418, normalized size of antiderivative = 1.40 \[ \int \frac {(e+f x) \cosh ^3(c+d x)}{a+b \sinh (c+d x)} \, dx=\frac {8 a b f \cosh (c+d x)+2 b^2 d (e+f x) \cosh (2 (c+d x))+4 \left (a^2+b^2\right ) \left (-2 d e (c+d x)+2 c f (c+d x)-f (c+d x)^2+\frac {4 a \sqrt {a^2+b^2} d e \arctan \left (\frac {a+b e^{c+d x}}{\sqrt {-a^2-b^2}}\right )}{\sqrt {-\left (a^2+b^2\right )^2}}-\frac {4 a \sqrt {-\left (a^2+b^2\right )^2} d e \text {arctanh}\left (\frac {a+b e^{c+d x}}{\sqrt {a^2+b^2}}\right )}{\left (-a^2-b^2\right )^{3/2}}+2 f (c+d x) \log \left (1+\frac {b e^{c+d x}}{a-\sqrt {a^2+b^2}}\right )+2 f (c+d x) \log \left (1+\frac {b e^{c+d x}}{a+\sqrt {a^2+b^2}}\right )-2 c f \log \left (b-2 a e^{c+d x}-b e^{2 (c+d x)}\right )+2 d e \log \left (2 a e^{c+d x}+b \left (-1+e^{2 (c+d x)}\right )\right )+2 f \operatorname {PolyLog}\left (2,\frac {b e^{c+d x}}{-a+\sqrt {a^2+b^2}}\right )+2 f \operatorname {PolyLog}\left (2,-\frac {b e^{c+d x}}{a+\sqrt {a^2+b^2}}\right )\right )-8 a b d (e+f x) \sinh (c+d x)-b^2 f \sinh (2 (c+d x))}{8 b^3 d^2} \] Input:

Integrate[((e + f*x)*Cosh[c + d*x]^3)/(a + b*Sinh[c + d*x]),x]
 

Output:

(8*a*b*f*Cosh[c + d*x] + 2*b^2*d*(e + f*x)*Cosh[2*(c + d*x)] + 4*(a^2 + b^ 
2)*(-2*d*e*(c + d*x) + 2*c*f*(c + d*x) - f*(c + d*x)^2 + (4*a*Sqrt[a^2 + b 
^2]*d*e*ArcTan[(a + b*E^(c + d*x))/Sqrt[-a^2 - b^2]])/Sqrt[-(a^2 + b^2)^2] 
 - (4*a*Sqrt[-(a^2 + b^2)^2]*d*e*ArcTanh[(a + b*E^(c + d*x))/Sqrt[a^2 + b^ 
2]])/(-a^2 - b^2)^(3/2) + 2*f*(c + d*x)*Log[1 + (b*E^(c + d*x))/(a - Sqrt[ 
a^2 + b^2])] + 2*f*(c + d*x)*Log[1 + (b*E^(c + d*x))/(a + Sqrt[a^2 + b^2]) 
] - 2*c*f*Log[b - 2*a*E^(c + d*x) - b*E^(2*(c + d*x))] + 2*d*e*Log[2*a*E^( 
c + d*x) + b*(-1 + E^(2*(c + d*x)))] + 2*f*PolyLog[2, (b*E^(c + d*x))/(-a 
+ Sqrt[a^2 + b^2])] + 2*f*PolyLog[2, -((b*E^(c + d*x))/(a + Sqrt[a^2 + b^2 
]))]) - 8*a*b*d*(e + f*x)*Sinh[c + d*x] - b^2*f*Sinh[2*(c + d*x)])/(8*b^3* 
d^2)
 

Rubi [A] (verified)

Time = 1.37 (sec) , antiderivative size = 274, normalized size of antiderivative = 0.92, number of steps used = 17, number of rules used = 16, \(\frac {\text {number of rules}}{\text {integrand size}}\) = 0.615, Rules used = {6099, 3042, 3777, 26, 3042, 26, 3118, 5969, 3042, 25, 3115, 24, 6095, 2620, 2715, 2838}

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

\(\Big \downarrow \) 6099

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

\(\Big \downarrow \) 3042

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

\(\Big \downarrow \) 3777

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

\(\Big \downarrow \) 26

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

\(\Big \downarrow \) 3042

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

\(\Big \downarrow \) 26

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

\(\Big \downarrow \) 3118

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

\(\Big \downarrow \) 5969

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

\(\Big \downarrow \) 3042

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

\(\Big \downarrow \) 25

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

\(\Big \downarrow \) 3115

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

\(\Big \downarrow \) 24

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

\(\Big \downarrow \) 6095

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

\(\Big \downarrow \) 2620

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

\(\Big \downarrow \) 2715

\(\displaystyle \frac {\left (a^2+b^2\right ) \left (-\frac {f \int e^{-c-d x} \log \left (\frac {e^{c+d x} b}{a-\sqrt {a^2+b^2}}+1\right )de^{c+d x}}{b d^2}-\frac {f \int e^{-c-d x} \log \left (\frac {e^{c+d x} b}{a+\sqrt {a^2+b^2}}+1\right )de^{c+d x}}{b d^2}+\frac {(e+f x) \log \left (\frac {b e^{c+d x}}{a-\sqrt {a^2+b^2}}+1\right )}{b d}+\frac {(e+f x) \log \left (\frac {b e^{c+d x}}{\sqrt {a^2+b^2}+a}+1\right )}{b d}-\frac {(e+f x)^2}{2 b f}\right )}{b^2}-\frac {a \left (\frac {(e+f x) \sinh (c+d x)}{d}-\frac {f \cosh (c+d x)}{d^2}\right )}{b^2}+\frac {\frac {(e+f x) \sinh ^2(c+d x)}{2 d}+\frac {f \left (\frac {x}{2}-\frac {\sinh (c+d x) \cosh (c+d x)}{2 d}\right )}{2 d}}{b}\)

\(\Big \downarrow \) 2838

\(\displaystyle \frac {\left (a^2+b^2\right ) \left (\frac {f \operatorname {PolyLog}\left (2,-\frac {b e^{c+d x}}{a-\sqrt {a^2+b^2}}\right )}{b d^2}+\frac {f \operatorname {PolyLog}\left (2,-\frac {b e^{c+d x}}{a+\sqrt {a^2+b^2}}\right )}{b d^2}+\frac {(e+f x) \log \left (\frac {b e^{c+d x}}{a-\sqrt {a^2+b^2}}+1\right )}{b d}+\frac {(e+f x) \log \left (\frac {b e^{c+d x}}{\sqrt {a^2+b^2}+a}+1\right )}{b d}-\frac {(e+f x)^2}{2 b f}\right )}{b^2}-\frac {a \left (\frac {(e+f x) \sinh (c+d x)}{d}-\frac {f \cosh (c+d x)}{d^2}\right )}{b^2}+\frac {\frac {(e+f x) \sinh ^2(c+d x)}{2 d}+\frac {f \left (\frac {x}{2}-\frac {\sinh (c+d x) \cosh (c+d x)}{2 d}\right )}{2 d}}{b}\)

Input:

Int[((e + f*x)*Cosh[c + d*x]^3)/(a + b*Sinh[c + d*x]),x]
 

Output:

((a^2 + b^2)*(-1/2*(e + f*x)^2/(b*f) + ((e + f*x)*Log[1 + (b*E^(c + d*x))/ 
(a - Sqrt[a^2 + b^2])])/(b*d) + ((e + f*x)*Log[1 + (b*E^(c + d*x))/(a + Sq 
rt[a^2 + b^2])])/(b*d) + (f*PolyLog[2, -((b*E^(c + d*x))/(a - Sqrt[a^2 + b 
^2]))])/(b*d^2) + (f*PolyLog[2, -((b*E^(c + d*x))/(a + Sqrt[a^2 + b^2]))]) 
/(b*d^2)))/b^2 - (a*(-((f*Cosh[c + d*x])/d^2) + ((e + f*x)*Sinh[c + d*x])/ 
d))/b^2 + (((e + f*x)*Sinh[c + d*x]^2)/(2*d) + (f*(x/2 - (Cosh[c + d*x]*Si 
nh[c + d*x])/(2*d)))/(2*d))/b
 

Defintions of rubi rules used

rule 24
Int[a_, x_Symbol] :> Simp[a*x, x] /; FreeQ[a, x]
 

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

rule 26
Int[(Complex[0, a_])*(Fx_), x_Symbol] :> Simp[(Complex[Identity[0], a])   I 
nt[Fx, x], x] /; FreeQ[a, x] && EqQ[a^2, 1]
 

rule 2620
Int[(((F_)^((g_.)*((e_.) + (f_.)*(x_))))^(n_.)*((c_.) + (d_.)*(x_))^(m_.))/ 
((a_) + (b_.)*((F_)^((g_.)*((e_.) + (f_.)*(x_))))^(n_.)), x_Symbol] :> Simp 
[((c + d*x)^m/(b*f*g*n*Log[F]))*Log[1 + b*((F^(g*(e + f*x)))^n/a)], x] - Si 
mp[d*(m/(b*f*g*n*Log[F]))   Int[(c + d*x)^(m - 1)*Log[1 + b*((F^(g*(e + f*x 
)))^n/a)], x], x] /; FreeQ[{F, a, b, c, d, e, f, g, n}, x] && IGtQ[m, 0]
 

rule 2715
Int[Log[(a_) + (b_.)*((F_)^((e_.)*((c_.) + (d_.)*(x_))))^(n_.)], x_Symbol] 
:> Simp[1/(d*e*n*Log[F])   Subst[Int[Log[a + b*x]/x, x], x, (F^(e*(c + d*x) 
))^n], x] /; FreeQ[{F, a, b, c, d, e, n}, x] && GtQ[a, 0]
 

rule 2838
Int[Log[(c_.)*((d_) + (e_.)*(x_)^(n_.))]/(x_), x_Symbol] :> Simp[-PolyLog[2 
, (-c)*e*x^n]/n, x] /; FreeQ[{c, d, e, n}, x] && EqQ[c*d, 1]
 

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

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

rule 3118
Int[sin[(c_.) + (d_.)*(x_)], x_Symbol] :> Simp[-Cos[c + d*x]/d, x] /; FreeQ 
[{c, d}, x]
 

rule 3777
Int[((c_.) + (d_.)*(x_))^(m_.)*sin[(e_.) + (f_.)*(x_)], x_Symbol] :> Simp[( 
-(c + d*x)^m)*(Cos[e + f*x]/f), x] + Simp[d*(m/f)   Int[(c + d*x)^(m - 1)*C 
os[e + f*x], x], x] /; FreeQ[{c, d, e, f}, x] && GtQ[m, 0]
 

rule 5969
Int[Cosh[(a_.) + (b_.)*(x_)]*((c_.) + (d_.)*(x_))^(m_.)*Sinh[(a_.) + (b_.)* 
(x_)]^(n_.), x_Symbol] :> Simp[(c + d*x)^m*(Sinh[a + b*x]^(n + 1)/(b*(n + 1 
))), x] - Simp[d*(m/(b*(n + 1)))   Int[(c + d*x)^(m - 1)*Sinh[a + b*x]^(n + 
 1), x], x] /; FreeQ[{a, b, c, d, n}, x] && IGtQ[m, 0] && NeQ[n, -1]
 

rule 6095
Int[(Cosh[(c_.) + (d_.)*(x_)]*((e_.) + (f_.)*(x_))^(m_.))/((a_) + (b_.)*Sin 
h[(c_.) + (d_.)*(x_)]), x_Symbol] :> Simp[-(e + f*x)^(m + 1)/(b*f*(m + 1)), 
 x] + (Int[(e + f*x)^m*(E^(c + d*x)/(a - Rt[a^2 + b^2, 2] + b*E^(c + d*x))) 
, x] + Int[(e + f*x)^m*(E^(c + d*x)/(a + Rt[a^2 + b^2, 2] + b*E^(c + d*x))) 
, x]) /; FreeQ[{a, b, c, d, e, f}, x] && IGtQ[m, 0] && NeQ[a^2 + b^2, 0]
 

rule 6099
Int[(Cosh[(c_.) + (d_.)*(x_)]^(n_)*((e_.) + (f_.)*(x_))^(m_.))/((a_) + (b_. 
)*Sinh[(c_.) + (d_.)*(x_)]), x_Symbol] :> Simp[-a/b^2   Int[(e + f*x)^m*Cos 
h[c + d*x]^(n - 2), x], x] + (Simp[1/b   Int[(e + f*x)^m*Cosh[c + d*x]^(n - 
 2)*Sinh[c + d*x], x], x] + Simp[(a^2 + b^2)/b^2   Int[(e + f*x)^m*(Cosh[c 
+ d*x]^(n - 2)/(a + b*Sinh[c + d*x])), x], x]) /; FreeQ[{a, b, c, d, e, f}, 
 x] && IGtQ[n, 1] && NeQ[a^2 + b^2, 0] && IGtQ[m, 0]
 
Maple [B] (verified)

Leaf count of result is larger than twice the leaf count of optimal. \(974\) vs. \(2(278)=556\).

Time = 5.82 (sec) , antiderivative size = 975, normalized size of antiderivative = 3.27

method result size
risch \(\frac {a^{2} f \ln \left (\frac {b \,{\mathrm e}^{d x +c}+\sqrt {a^{2}+b^{2}}+a}{a +\sqrt {a^{2}+b^{2}}}\right ) c}{d^{2} b^{3}}+\frac {2 c \,a^{2} f \ln \left ({\mathrm e}^{d x +c}\right )}{d^{2} b^{3}}-\frac {a^{2} f \,c^{2}}{d^{2} b^{3}}-\frac {2 a^{2} e \ln \left ({\mathrm e}^{d x +c}\right )}{d \,b^{3}}+\frac {a^{2} e \ln \left (b \,{\mathrm e}^{2 d x +2 c}+2 a \,{\mathrm e}^{d x +c}-b \right )}{d \,b^{3}}+\frac {a^{2} f \operatorname {dilog}\left (\frac {-b \,{\mathrm e}^{d x +c}+\sqrt {a^{2}+b^{2}}-a}{-a +\sqrt {a^{2}+b^{2}}}\right )}{d^{2} b^{3}}+\frac {a^{2} f \operatorname {dilog}\left (\frac {b \,{\mathrm e}^{d x +c}+\sqrt {a^{2}+b^{2}}+a}{a +\sqrt {a^{2}+b^{2}}}\right )}{d^{2} b^{3}}-\frac {f \,c^{2}}{d^{2} b}+\frac {e \ln \left (b \,{\mathrm e}^{2 d x +2 c}+2 a \,{\mathrm e}^{d x +c}-b \right )}{d b}-\frac {2 e \ln \left ({\mathrm e}^{d x +c}\right )}{d b}+\frac {f \operatorname {dilog}\left (\frac {b \,{\mathrm e}^{d x +c}+\sqrt {a^{2}+b^{2}}+a}{a +\sqrt {a^{2}+b^{2}}}\right )}{d^{2} b}+\frac {f \operatorname {dilog}\left (\frac {-b \,{\mathrm e}^{d x +c}+\sqrt {a^{2}+b^{2}}-a}{-a +\sqrt {a^{2}+b^{2}}}\right )}{d^{2} b}-\frac {2 a^{2} f c x}{d \,b^{3}}+\frac {a^{2} f \ln \left (\frac {-b \,{\mathrm e}^{d x +c}+\sqrt {a^{2}+b^{2}}-a}{-a +\sqrt {a^{2}+b^{2}}}\right ) x}{d \,b^{3}}+\frac {a^{2} f \ln \left (\frac {b \,{\mathrm e}^{d x +c}+\sqrt {a^{2}+b^{2}}+a}{a +\sqrt {a^{2}+b^{2}}}\right ) x}{d \,b^{3}}+\frac {a^{2} f \ln \left (\frac {-b \,{\mathrm e}^{d x +c}+\sqrt {a^{2}+b^{2}}-a}{-a +\sqrt {a^{2}+b^{2}}}\right ) c}{d^{2} b^{3}}-\frac {a \left (d x f +d e -f \right ) {\mathrm e}^{d x +c}}{2 b^{2} d^{2}}-\frac {c \,a^{2} f \ln \left (b \,{\mathrm e}^{2 d x +2 c}+2 a \,{\mathrm e}^{d x +c}-b \right )}{d^{2} b^{3}}-\frac {2 f c x}{d b}+\frac {f \ln \left (\frac {-b \,{\mathrm e}^{d x +c}+\sqrt {a^{2}+b^{2}}-a}{-a +\sqrt {a^{2}+b^{2}}}\right ) c}{d^{2} b}+\frac {f \ln \left (\frac {b \,{\mathrm e}^{d x +c}+\sqrt {a^{2}+b^{2}}+a}{a +\sqrt {a^{2}+b^{2}}}\right ) c}{d^{2} b}+\frac {2 c f \ln \left ({\mathrm e}^{d x +c}\right )}{d^{2} b}-\frac {c f \ln \left (b \,{\mathrm e}^{2 d x +2 c}+2 a \,{\mathrm e}^{d x +c}-b \right )}{d^{2} b}+\frac {f \ln \left (\frac {-b \,{\mathrm e}^{d x +c}+\sqrt {a^{2}+b^{2}}-a}{-a +\sqrt {a^{2}+b^{2}}}\right ) x}{d b}+\frac {f \ln \left (\frac {b \,{\mathrm e}^{d x +c}+\sqrt {a^{2}+b^{2}}+a}{a +\sqrt {a^{2}+b^{2}}}\right ) x}{d b}-\frac {a^{2} f \,x^{2}}{2 b^{3}}-\frac {f \,x^{2}}{2 b}+\frac {a^{2} e x}{b^{3}}+\frac {\left (2 d x f +2 d e -f \right ) {\mathrm e}^{2 d x +2 c}}{16 b \,d^{2}}+\frac {\left (2 d x f +2 d e +f \right ) {\mathrm e}^{-2 d x -2 c}}{16 b \,d^{2}}+\frac {e x}{b}+\frac {a \left (d x f +d e +f \right ) {\mathrm e}^{-d x -c}}{2 b^{2} d^{2}}\) \(975\)

Input:

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

Output:

1/d^2/b^3*a^2*f*ln((b*exp(d*x+c)+(a^2+b^2)^(1/2)+a)/(a+(a^2+b^2)^(1/2)))*c 
+2/d^2/b^3*c*a^2*f*ln(exp(d*x+c))-1/d^2/b^3*a^2*f*c^2-2/d/b^3*a^2*e*ln(exp 
(d*x+c))+1/d/b^3*a^2*e*ln(b*exp(2*d*x+2*c)+2*a*exp(d*x+c)-b)+1/d^2/b^3*a^2 
*f*dilog((-b*exp(d*x+c)+(a^2+b^2)^(1/2)-a)/(-a+(a^2+b^2)^(1/2)))+1/d^2/b^3 
*a^2*f*dilog((b*exp(d*x+c)+(a^2+b^2)^(1/2)+a)/(a+(a^2+b^2)^(1/2)))-1/d^2/b 
*f*c^2+1/d/b*e*ln(b*exp(2*d*x+2*c)+2*a*exp(d*x+c)-b)-2/d/b*e*ln(exp(d*x+c) 
)+1/d^2/b*f*dilog((b*exp(d*x+c)+(a^2+b^2)^(1/2)+a)/(a+(a^2+b^2)^(1/2)))+1/ 
d^2/b*f*dilog((-b*exp(d*x+c)+(a^2+b^2)^(1/2)-a)/(-a+(a^2+b^2)^(1/2)))-2/d/ 
b^3*a^2*f*c*x+1/d/b^3*a^2*f*ln((-b*exp(d*x+c)+(a^2+b^2)^(1/2)-a)/(-a+(a^2+ 
b^2)^(1/2)))*x+1/d/b^3*a^2*f*ln((b*exp(d*x+c)+(a^2+b^2)^(1/2)+a)/(a+(a^2+b 
^2)^(1/2)))*x+1/d^2/b^3*a^2*f*ln((-b*exp(d*x+c)+(a^2+b^2)^(1/2)-a)/(-a+(a^ 
2+b^2)^(1/2)))*c-1/2*a*(d*f*x+d*e-f)/b^2/d^2*exp(d*x+c)-1/d^2/b^3*c*a^2*f* 
ln(b*exp(2*d*x+2*c)+2*a*exp(d*x+c)-b)-2/d/b*f*c*x+1/d^2/b*f*ln((-b*exp(d*x 
+c)+(a^2+b^2)^(1/2)-a)/(-a+(a^2+b^2)^(1/2)))*c+1/d^2/b*f*ln((b*exp(d*x+c)+ 
(a^2+b^2)^(1/2)+a)/(a+(a^2+b^2)^(1/2)))*c+2/d^2/b*c*f*ln(exp(d*x+c))-1/d^2 
/b*c*f*ln(b*exp(2*d*x+2*c)+2*a*exp(d*x+c)-b)+1/d/b*f*ln((-b*exp(d*x+c)+(a^ 
2+b^2)^(1/2)-a)/(-a+(a^2+b^2)^(1/2)))*x+1/d/b*f*ln((b*exp(d*x+c)+(a^2+b^2) 
^(1/2)+a)/(a+(a^2+b^2)^(1/2)))*x-1/2/b^3*a^2*f*x^2-1/2/b*f*x^2+1/b^3*a^2*e 
*x+1/16*(2*d*f*x+2*d*e-f)/b/d^2*exp(2*d*x+2*c)+1/16*(2*d*f*x+2*d*e+f)/b/d^ 
2*exp(-2*d*x-2*c)+1/b*e*x+1/2*a*(d*f*x+d*e+f)/b^2/d^2*exp(-d*x-c)
 

Fricas [B] (verification not implemented)

Leaf count of result is larger than twice the leaf count of optimal. 1416 vs. \(2 (276) = 552\).

Time = 0.13 (sec) , antiderivative size = 1416, normalized size of antiderivative = 4.75 \[ \int \frac {(e+f x) \cosh ^3(c+d x)}{a+b \sinh (c+d x)} \, dx=\text {Too large to display} \] Input:

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

Output:

1/16*(2*b^2*d*f*x + (2*b^2*d*f*x + 2*b^2*d*e - b^2*f)*cosh(d*x + c)^4 + (2 
*b^2*d*f*x + 2*b^2*d*e - b^2*f)*sinh(d*x + c)^4 + 2*b^2*d*e - 8*(a*b*d*f*x 
 + a*b*d*e - a*b*f)*cosh(d*x + c)^3 - 4*(2*a*b*d*f*x + 2*a*b*d*e - 2*a*b*f 
 - (2*b^2*d*f*x + 2*b^2*d*e - b^2*f)*cosh(d*x + c))*sinh(d*x + c)^3 + b^2* 
f - 8*((a^2 + b^2)*d^2*f*x^2 + 2*(a^2 + b^2)*d^2*e*x + 4*(a^2 + b^2)*c*d*e 
 - 2*(a^2 + b^2)*c^2*f)*cosh(d*x + c)^2 - 2*(4*(a^2 + b^2)*d^2*f*x^2 + 8*( 
a^2 + b^2)*d^2*e*x + 16*(a^2 + b^2)*c*d*e - 8*(a^2 + b^2)*c^2*f - 3*(2*b^2 
*d*f*x + 2*b^2*d*e - b^2*f)*cosh(d*x + c)^2 + 12*(a*b*d*f*x + a*b*d*e - a* 
b*f)*cosh(d*x + c))*sinh(d*x + c)^2 + 8*(a*b*d*f*x + a*b*d*e + a*b*f)*cosh 
(d*x + c) + 16*((a^2 + b^2)*f*cosh(d*x + c)^2 + 2*(a^2 + b^2)*f*cosh(d*x + 
 c)*sinh(d*x + c) + (a^2 + b^2)*f*sinh(d*x + c)^2)*dilog((a*cosh(d*x + c) 
+ a*sinh(d*x + c) + (b*cosh(d*x + c) + b*sinh(d*x + c))*sqrt((a^2 + b^2)/b 
^2) - b)/b + 1) + 16*((a^2 + b^2)*f*cosh(d*x + c)^2 + 2*(a^2 + b^2)*f*cosh 
(d*x + c)*sinh(d*x + c) + (a^2 + b^2)*f*sinh(d*x + c)^2)*dilog((a*cosh(d*x 
 + c) + a*sinh(d*x + c) - (b*cosh(d*x + c) + b*sinh(d*x + c))*sqrt((a^2 + 
b^2)/b^2) - b)/b + 1) + 16*(((a^2 + b^2)*d*e - (a^2 + b^2)*c*f)*cosh(d*x + 
 c)^2 + 2*((a^2 + b^2)*d*e - (a^2 + b^2)*c*f)*cosh(d*x + c)*sinh(d*x + c) 
+ ((a^2 + b^2)*d*e - (a^2 + b^2)*c*f)*sinh(d*x + c)^2)*log(2*b*cosh(d*x + 
c) + 2*b*sinh(d*x + c) + 2*b*sqrt((a^2 + b^2)/b^2) + 2*a) + 16*(((a^2 + b^ 
2)*d*e - (a^2 + b^2)*c*f)*cosh(d*x + c)^2 + 2*((a^2 + b^2)*d*e - (a^2 +...
                                                                                    
                                                                                    
 

Sympy [F(-1)]

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

integrate((f*x+e)*cosh(d*x+c)**3/(a+b*sinh(d*x+c)),x)
 

Output:

Timed out
 

Maxima [F]

\[ \int \frac {(e+f x) \cosh ^3(c+d x)}{a+b \sinh (c+d x)} \, dx=\int { \frac {{\left (f x + e\right )} \cosh \left (d x + c\right )^{3}}{b \sinh \left (d x + c\right ) + a} \,d x } \] Input:

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

Output:

-1/8*e*((4*a*e^(-d*x - c) - b)*e^(2*d*x + 2*c)/(b^2*d) - 8*(a^2 + b^2)*(d* 
x + c)/(b^3*d) - (4*a*e^(-d*x - c) + b*e^(-2*d*x - 2*c))/(b^2*d) - 8*(a^2 
+ b^2)*log(-2*a*e^(-d*x - c) + b*e^(-2*d*x - 2*c) - b)/(b^3*d)) + 1/16*f*( 
(8*(a^2*d^2*e^(2*c) + b^2*d^2*e^(2*c))*x^2 + (2*b^2*d*x*e^(4*c) - b^2*e^(4 
*c))*e^(2*d*x) - 8*(a*b*d*x*e^(3*c) - a*b*e^(3*c))*e^(d*x) + 8*(a*b*d*x*e^ 
c + a*b*e^c)*e^(-d*x) + (2*b^2*d*x + b^2)*e^(-2*d*x))*e^(-2*c)/(b^3*d^2) - 
 2*integrate(16*((a^3*e^c + a*b^2*e^c)*x*e^(d*x) - (a^2*b + b^3)*x)/(b^4*e 
^(2*d*x + 2*c) + 2*a*b^3*e^(d*x + c) - b^4), x))
 

Giac [F]

\[ \int \frac {(e+f x) \cosh ^3(c+d x)}{a+b \sinh (c+d x)} \, dx=\int { \frac {{\left (f x + e\right )} \cosh \left (d x + c\right )^{3}}{b \sinh \left (d x + c\right ) + a} \,d x } \] Input:

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

Output:

integrate((f*x + e)*cosh(d*x + c)^3/(b*sinh(d*x + c) + a), x)
 

Mupad [F(-1)]

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

int((cosh(c + d*x)^3*(e + f*x))/(a + b*sinh(c + d*x)),x)
 

Output:

int((cosh(c + d*x)^3*(e + f*x))/(a + b*sinh(c + d*x)), x)
 

Reduce [F]

\[ \int \frac {(e+f x) \cosh ^3(c+d x)}{a+b \sinh (c+d x)} \, dx=\frac {2 e^{4 d x +4 c} b^{4} d e -24 a^{2} b^{2} f +2 b^{4} d e -e^{4 d x +4 c} b^{4} f -48 a^{2} b^{2} d f x -8 e^{3 d x +3 c} a \,b^{3} d f x +16 e^{2 d x +2 c} \mathrm {log}\left (e^{2 d x +2 c} b +2 e^{d x +c} a -b \right ) a^{2} b^{2} d e -16 e^{2 d x +2 c} a^{2} b^{2} d^{2} e x +8 e^{2 d x +2 c} a^{2} b^{2} d^{2} f \,x^{2}+32 e^{d x +c} a^{3} b d f x +40 e^{d x +c} a \,b^{3} d f x +64 e^{2 d x +2 c} \left (\int \frac {x}{e^{4 d x +4 c} b +2 e^{3 d x +3 c} a -e^{2 d x +2 c} b}d x \right ) a^{4} b \,d^{2} f +96 e^{2 d x +2 c} \left (\int \frac {x}{e^{4 d x +4 c} b +2 e^{3 d x +3 c} a -e^{2 d x +2 c} b}d x \right ) a^{2} b^{3} d^{2} f -224 e^{2 d x +c} \left (\int \frac {x}{e^{3 d x +2 c} b +2 e^{2 d x +c} a -e^{d x} b}d x \right ) a^{3} b^{2} d^{2} f -96 e^{2 d x +c} \left (\int \frac {x}{e^{3 d x +2 c} b +2 e^{2 d x +c} a -e^{d x} b}d x \right ) a \,b^{4} d^{2} f +8 e^{3 d x +3 c} a \,b^{3} f +32 e^{d x +c} a^{3} b f +40 e^{d x +c} a \,b^{3} f -32 a^{4} d f x -14 b^{4} d f x -16 a^{4} f -7 b^{4} f +2 e^{4 d x +4 c} b^{4} d f x -8 e^{3 d x +3 c} a \,b^{3} d e +16 e^{2 d x +2 c} \mathrm {log}\left (e^{2 d x +2 c} b +2 e^{d x +c} a -b \right ) b^{4} d e -16 e^{2 d x +2 c} b^{4} d^{2} e x +8 e^{2 d x +2 c} b^{4} d^{2} f \,x^{2}+8 e^{d x +c} a \,b^{3} d e +32 e^{2 d x +2 c} \left (\int \frac {x}{e^{4 d x +4 c} b +2 e^{3 d x +3 c} a -e^{2 d x +2 c} b}d x \right ) b^{5} d^{2} f -128 e^{2 d x +c} \left (\int \frac {x}{e^{3 d x +2 c} b +2 e^{2 d x +c} a -e^{d x} b}d x \right ) a^{5} d^{2} f}{16 e^{2 d x +2 c} b^{5} d^{2}} \] Input:

int((f*x+e)*cosh(d*x+c)^3/(a+b*sinh(d*x+c)),x)
 

Output:

(2*e**(4*c + 4*d*x)*b**4*d*e + 2*e**(4*c + 4*d*x)*b**4*d*f*x - e**(4*c + 4 
*d*x)*b**4*f - 8*e**(3*c + 3*d*x)*a*b**3*d*e - 8*e**(3*c + 3*d*x)*a*b**3*d 
*f*x + 8*e**(3*c + 3*d*x)*a*b**3*f + 64*e**(2*c + 2*d*x)*int(x/(e**(4*c + 
4*d*x)*b + 2*e**(3*c + 3*d*x)*a - e**(2*c + 2*d*x)*b),x)*a**4*b*d**2*f + 9 
6*e**(2*c + 2*d*x)*int(x/(e**(4*c + 4*d*x)*b + 2*e**(3*c + 3*d*x)*a - e**( 
2*c + 2*d*x)*b),x)*a**2*b**3*d**2*f + 32*e**(2*c + 2*d*x)*int(x/(e**(4*c + 
 4*d*x)*b + 2*e**(3*c + 3*d*x)*a - e**(2*c + 2*d*x)*b),x)*b**5*d**2*f + 16 
*e**(2*c + 2*d*x)*log(e**(2*c + 2*d*x)*b + 2*e**(c + d*x)*a - b)*a**2*b**2 
*d*e + 16*e**(2*c + 2*d*x)*log(e**(2*c + 2*d*x)*b + 2*e**(c + d*x)*a - b)* 
b**4*d*e - 16*e**(2*c + 2*d*x)*a**2*b**2*d**2*e*x + 8*e**(2*c + 2*d*x)*a** 
2*b**2*d**2*f*x**2 - 16*e**(2*c + 2*d*x)*b**4*d**2*e*x + 8*e**(2*c + 2*d*x 
)*b**4*d**2*f*x**2 - 128*e**(c + 2*d*x)*int(x/(e**(2*c + 3*d*x)*b + 2*e**( 
c + 2*d*x)*a - e**(d*x)*b),x)*a**5*d**2*f - 224*e**(c + 2*d*x)*int(x/(e**( 
2*c + 3*d*x)*b + 2*e**(c + 2*d*x)*a - e**(d*x)*b),x)*a**3*b**2*d**2*f - 96 
*e**(c + 2*d*x)*int(x/(e**(2*c + 3*d*x)*b + 2*e**(c + 2*d*x)*a - e**(d*x)* 
b),x)*a*b**4*d**2*f + 32*e**(c + d*x)*a**3*b*d*f*x + 32*e**(c + d*x)*a**3* 
b*f + 8*e**(c + d*x)*a*b**3*d*e + 40*e**(c + d*x)*a*b**3*d*f*x + 40*e**(c 
+ d*x)*a*b**3*f - 32*a**4*d*f*x - 16*a**4*f - 48*a**2*b**2*d*f*x - 24*a**2 
*b**2*f + 2*b**4*d*e - 14*b**4*d*f*x - 7*b**4*f)/(16*e**(2*c + 2*d*x)*b**5 
*d**2)