\(\int \frac {\sinh ^2(x)}{a+b \sinh (x)+c \sinh ^2(x)} \, dx\) [828]

Optimal result

Integrand size = 19, antiderivative size = 309 \[ \int \frac {\sinh ^2(x)}{a+b \sinh (x)+c \sinh ^2(x)} \, dx=\frac {x}{c}-\frac {\sqrt {2} \left (i b+\frac {b^2-2 a c}{\sqrt {-b^2+4 a c}}\right ) \arctan \left (\frac {2 i c-\left (i b-\sqrt {-b^2+4 a c}\right ) \tanh \left (\frac {x}{2}\right )}{\sqrt {2} \sqrt {b^2-2 (a-c) c+i b \sqrt {-b^2+4 a c}}}\right )}{c \sqrt {b^2-2 (a-c) c+i b \sqrt {-b^2+4 a c}}}-\frac {\sqrt {2} \left (i b-\frac {b^2-2 a c}{\sqrt {-b^2+4 a c}}\right ) \arctan \left (\frac {2 i c-\left (i b+\sqrt {-b^2+4 a c}\right ) \tanh \left (\frac {x}{2}\right )}{\sqrt {2} \sqrt {b^2-2 (a-c) c-i b \sqrt {-b^2+4 a c}}}\right )}{c \sqrt {b^2-2 (a-c) c-i b \sqrt {-b^2+4 a c}}} \]

[Out]

x/c-arctan(1/2*(2*I*c-(I*b+(4*a*c-b^2)^(1/2))*tanh(1/2*x))*2^(1/2)/(b^2-2*(a-c)*c-I*b*(4*a*c-b^2)^(1/2))^(1/2)
)*2^(1/2)*(I*b+(2*a*c-b^2)/(4*a*c-b^2)^(1/2))/c/(b^2-2*(a-c)*c-I*b*(4*a*c-b^2)^(1/2))^(1/2)-arctan(1/2*(2*I*c-
(I*b-(4*a*c-b^2)^(1/2))*tanh(1/2*x))*2^(1/2)/(b^2-2*(a-c)*c+I*b*(4*a*c-b^2)^(1/2))^(1/2))*2^(1/2)*(I*b+(-2*a*c
+b^2)/(4*a*c-b^2)^(1/2))/c/(b^2-2*(a-c)*c+I*b*(4*a*c-b^2)^(1/2))^(1/2)

Rubi [A] (verified)

Time = 0.69 (sec) , antiderivative size = 309, normalized size of antiderivative = 1.00, number of steps used = 9, number of rules used = 5, \(\frac {\text {number of rules}}{\text {integrand size}}\) = 0.263, Rules used = {3337, 3373, 2739, 632, 210} \[ \int \frac {\sinh ^2(x)}{a+b \sinh (x)+c \sinh ^2(x)} \, dx=-\frac {\sqrt {2} \left (\frac {b^2-2 a c}{\sqrt {4 a c-b^2}}+i b\right ) \arctan \left (\frac {2 i c-\tanh \left (\frac {x}{2}\right ) \left (-\sqrt {4 a c-b^2}+i b\right )}{\sqrt {2} \sqrt {i b \sqrt {4 a c-b^2}-2 c (a-c)+b^2}}\right )}{c \sqrt {i b \sqrt {4 a c-b^2}-2 c (a-c)+b^2}}-\frac {\sqrt {2} \left (-\frac {b^2-2 a c}{\sqrt {4 a c-b^2}}+i b\right ) \arctan \left (\frac {2 i c-\tanh \left (\frac {x}{2}\right ) \left (\sqrt {4 a c-b^2}+i b\right )}{\sqrt {2} \sqrt {-i b \sqrt {4 a c-b^2}-2 c (a-c)+b^2}}\right )}{c \sqrt {-i b \sqrt {4 a c-b^2}-2 c (a-c)+b^2}}+\frac {x}{c} \]

[In]

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

[Out]

x/c - (Sqrt[2]*(I*b + (b^2 - 2*a*c)/Sqrt[-b^2 + 4*a*c])*ArcTan[((2*I)*c - (I*b - Sqrt[-b^2 + 4*a*c])*Tanh[x/2]
)/(Sqrt[2]*Sqrt[b^2 - 2*(a - c)*c + I*b*Sqrt[-b^2 + 4*a*c]])])/(c*Sqrt[b^2 - 2*(a - c)*c + I*b*Sqrt[-b^2 + 4*a
*c]]) - (Sqrt[2]*(I*b - (b^2 - 2*a*c)/Sqrt[-b^2 + 4*a*c])*ArcTan[((2*I)*c - (I*b + Sqrt[-b^2 + 4*a*c])*Tanh[x/
2])/(Sqrt[2]*Sqrt[b^2 - 2*(a - c)*c - I*b*Sqrt[-b^2 + 4*a*c]])])/(c*Sqrt[b^2 - 2*(a - c)*c - I*b*Sqrt[-b^2 + 4
*a*c]])

Rule 210

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 632

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

Rule 2739

Int[((a_) + (b_.)*sin[(c_.) + (d_.)*(x_)])^(-1), x_Symbol] :> With[{e = FreeFactors[Tan[(c + d*x)/2], x]}, Dis
t[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 3337

Int[sin[(d_.) + (e_.)*(x_)]^(m_.)*((a_.) + (b_.)*sin[(d_.) + (e_.)*(x_)]^(n_.) + (c_.)*sin[(d_.) + (e_.)*(x_)]
^(n2_.))^(p_), x_Symbol] :> Int[ExpandTrig[sin[d + e*x]^m*(a + b*sin[d + e*x]^n + c*sin[d + e*x]^(2*n))^p, x],
x] /; FreeQ[{a, b, c, d, e}, x] && EqQ[n2, 2*n] && NeQ[b^2 - 4*a*c, 0] && IntegersQ[m, n, p]

Rule 3373

Int[((A_) + (B_.)*sin[(d_.) + (e_.)*(x_)])/((a_.) + (b_.)*sin[(d_.) + (e_.)*(x_)] + (c_.)*sin[(d_.) + (e_.)*(x
_)]^2), x_Symbol] :> Module[{q = Rt[b^2 - 4*a*c, 2]}, Dist[B + (b*B - 2*A*c)/q, Int[1/(b + q + 2*c*Sin[d + e*x
]), x], x] + Dist[B - (b*B - 2*A*c)/q, Int[1/(b - q + 2*c*Sin[d + e*x]), x], x]] /; FreeQ[{a, b, c, d, e, A, B
}, x] && NeQ[b^2 - 4*a*c, 0]

Rubi steps \begin{align*} \text {integral}& = -\int \left (-\frac {1}{c}+\frac {-a-b \sinh (x)}{c \left (-a-b \sinh (x)-c \sinh ^2(x)\right )}\right ) \, dx \\ & = \frac {x}{c}-\frac {\int \frac {-a-b \sinh (x)}{-a-b \sinh (x)-c \sinh ^2(x)} \, dx}{c} \\ & = \frac {x}{c}-\frac {\left (i b-\frac {b^2-2 a c}{\sqrt {-b^2+4 a c}}\right ) \int \frac {1}{i b+\sqrt {-b^2+4 a c}+2 i c \sinh (x)} \, dx}{c}-\frac {\left (i b+\frac {b^2-2 a c}{\sqrt {-b^2+4 a c}}\right ) \int \frac {1}{i b-\sqrt {-b^2+4 a c}+2 i c \sinh (x)} \, dx}{c} \\ & = \frac {x}{c}-\frac {\left (2 \left (i b-\frac {b^2-2 a c}{\sqrt {-b^2+4 a c}}\right )\right ) \text {Subst}\left (\int \frac {1}{i b+\sqrt {-b^2+4 a c}+4 i c x-\left (i b+\sqrt {-b^2+4 a c}\right ) x^2} \, dx,x,\tanh \left (\frac {x}{2}\right )\right )}{c}-\frac {\left (2 \left (i b+\frac {b^2-2 a c}{\sqrt {-b^2+4 a c}}\right )\right ) \text {Subst}\left (\int \frac {1}{i b-\sqrt {-b^2+4 a c}+4 i c x-\left (i b-\sqrt {-b^2+4 a c}\right ) x^2} \, dx,x,\tanh \left (\frac {x}{2}\right )\right )}{c} \\ & = \frac {x}{c}+\frac {\left (4 \left (i b-\frac {b^2-2 a c}{\sqrt {-b^2+4 a c}}\right )\right ) \text {Subst}\left (\int \frac {1}{-8 \left (b^2-2 (a-c) c-i b \sqrt {-b^2+4 a c}\right )-x^2} \, dx,x,4 i c+2 \left (-i b-\sqrt {-b^2+4 a c}\right ) \tanh \left (\frac {x}{2}\right )\right )}{c}+\frac {\left (4 \left (i b+\frac {b^2-2 a c}{\sqrt {-b^2+4 a c}}\right )\right ) \text {Subst}\left (\int \frac {1}{-8 \left (b^2-2 (a-c) c+i b \sqrt {-b^2+4 a c}\right )-x^2} \, dx,x,4 i c+2 \left (-i b+\sqrt {-b^2+4 a c}\right ) \tanh \left (\frac {x}{2}\right )\right )}{c} \\ & = \frac {x}{c}-\frac {\sqrt {2} \left (i b+\frac {b^2-2 a c}{\sqrt {-b^2+4 a c}}\right ) \arctan \left (\frac {2 i c-\left (i b-\sqrt {-b^2+4 a c}\right ) \tanh \left (\frac {x}{2}\right )}{\sqrt {2} \sqrt {b^2-2 (a-c) c+i b \sqrt {-b^2+4 a c}}}\right )}{c \sqrt {b^2-2 (a-c) c+i b \sqrt {-b^2+4 a c}}}-\frac {\sqrt {2} \left (i b-\frac {b^2-2 a c}{\sqrt {-b^2+4 a c}}\right ) \arctan \left (\frac {2 i c-\left (i b+\sqrt {-b^2+4 a c}\right ) \tanh \left (\frac {x}{2}\right )}{\sqrt {2} \sqrt {b^2-2 (a-c) c-i b \sqrt {-b^2+4 a c}}}\right )}{c \sqrt {b^2-2 (a-c) c-i b \sqrt {-b^2+4 a c}}} \\ \end{align*}

Mathematica [A] (verified)

Time = 2.60 (sec) , antiderivative size = 283, normalized size of antiderivative = 0.92 \[ \int \frac {\sinh ^2(x)}{a+b \sinh (x)+c \sinh ^2(x)} \, dx=\frac {x-\frac {\sqrt {2} \left (-b^2+2 a c+b \sqrt {b^2-4 a c}\right ) \arctan \left (\frac {2 c+\left (-b+\sqrt {b^2-4 a c}\right ) \tanh \left (\frac {x}{2}\right )}{\sqrt {-2 b^2+4 (a-c) c+2 b \sqrt {b^2-4 a c}}}\right )}{\sqrt {b^2-4 a c} \sqrt {-b^2+2 (a-c) c+b \sqrt {b^2-4 a c}}}-\frac {\sqrt {2} \left (b^2-2 a c+b \sqrt {b^2-4 a c}\right ) \arctan \left (\frac {2 c-\left (b+\sqrt {b^2-4 a c}\right ) \tanh \left (\frac {x}{2}\right )}{\sqrt {2} \sqrt {-b^2+2 (a-c) c-b \sqrt {b^2-4 a c}}}\right )}{\sqrt {b^2-4 a c} \sqrt {-b^2+2 (a-c) c-b \sqrt {b^2-4 a c}}}}{c} \]

[In]

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

[Out]

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

Maple [C] (verified)

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

Time = 0.90 (sec) , antiderivative size = 108, normalized size of antiderivative = 0.35

method	result	size
default	\(\frac {\munderset {\textit {\_R} =\operatorname {RootOf}\left (a \,\textit {\_Z}^{4}-2 b \,\textit {\_Z}^{3}+\left (-2 a +4 c \right ) \textit {\_Z}^{2}+2 b \textit {\_Z} +a \right )}{\sum }\frac {\left (\textit {\_R}^{2} a -2 \textit {\_R} b -a \right ) \ln \left (\tanh \left (\frac {x}{2}\right )-\textit {\_R} \right )}{2 \textit {\_R}^{3} a -3 \textit {\_R}^{2} b -2 \textit {\_R} a +4 \textit {\_R} c +b}}{c}-\frac {\ln \left (\tanh \left (\frac {x}{2}\right )-1\right )}{c}+\frac {\ln \left (1+\tanh \left (\frac {x}{2}\right )\right )}{c}\)	\(108\)
risch	\(\text {Expression too large to display}\)	\(1126\)

[In]

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

[Out]

1/c*sum((_R^2*a-2*_R*b-a)/(2*_R^3*a-3*_R^2*b-2*_R*a+4*_R*c+b)*ln(tanh(1/2*x)-_R),_R=RootOf(a*_Z^4-2*b*_Z^3+(-2
*a+4*c)*_Z^2+2*b*_Z+a))-1/c*ln(tanh(1/2*x)-1)+1/c*ln(1+tanh(1/2*x))

Fricas [B] (verification not implemented)

Leaf count of result is larger than twice the leaf count of optimal. 4943 vs. \(2 (253) = 506\).

Time = 0.60 (sec) , antiderivative size = 4943, normalized size of antiderivative = 16.00 \[ \int \frac {\sinh ^2(x)}{a+b \sinh (x)+c \sinh ^2(x)} \, dx=\text {Too large to display} \]

[In]

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

[Out]

1/2*(sqrt(2)*c*sqrt(-(a^2*b^2 + b^4 + 2*a^2*c^2 - 2*(a^3 + 2*a*b^2)*c + (4*a*c^5 - (8*a^2 + b^2)*c^4 + 2*(2*a^
3 + 3*a*b^2)*c^3 - (a^2*b^2 + b^4)*c^2)*sqrt(-(a^4*b^2 + 2*a^2*b^4 + b^6 + 4*a^2*b^2*c^2 - 4*(a^3*b^2 + a*b^4)
*c)/(4*a*c^9 - (16*a^2 + b^2)*c^8 + 12*(2*a^3 + a*b^2)*c^7 - 2*(8*a^4 + 11*a^2*b^2 + b^4)*c^6 + 4*(a^5 + 3*a^3
*b^2 + 2*a*b^4)*c^5 - (a^4*b^2 + 2*a^2*b^4 + b^6)*c^4)))/(4*a*c^5 - (8*a^2 + b^2)*c^4 + 2*(2*a^3 + 3*a*b^2)*c^
3 - (a^2*b^2 + b^4)*c^2))*log(-2*a^4*b^2 - 2*a^2*b^4 + 4*a^3*b^2*c + sqrt(2)*(8*a^2*b^2*c^3 - 2*(2*a^3*b^2 + 3
*a*b^4)*c^2 + (a^2*b^4 + b^6)*c + (8*a^2*c^7 - 6*(4*a^3 + a*b^2)*c^6 + (24*a^4 + 22*a^2*b^2 + b^4)*c^5 - 2*(4*
a^5 + 9*a^3*b^2 + 4*a*b^4)*c^4 + (2*a^4*b^2 + 3*a^2*b^4 + b^6)*c^3)*sqrt(-(a^4*b^2 + 2*a^2*b^4 + b^6 + 4*a^2*b
^2*c^2 - 4*(a^3*b^2 + a*b^4)*c)/(4*a*c^9 - (16*a^2 + b^2)*c^8 + 12*(2*a^3 + a*b^2)*c^7 - 2*(8*a^4 + 11*a^2*b^2
+ b^4)*c^6 + 4*(a^5 + 3*a^3*b^2 + 2*a*b^4)*c^5 - (a^4*b^2 + 2*a^2*b^4 + b^6)*c^4)))*sqrt(-(a^2*b^2 + b^4 + 2*
a^2*c^2 - 2*(a^3 + 2*a*b^2)*c + (4*a*c^5 - (8*a^2 + b^2)*c^4 + 2*(2*a^3 + 3*a*b^2)*c^3 - (a^2*b^2 + b^4)*c^2)*
sqrt(-(a^4*b^2 + 2*a^2*b^4 + b^6 + 4*a^2*b^2*c^2 - 4*(a^3*b^2 + a*b^4)*c)/(4*a*c^9 - (16*a^2 + b^2)*c^8 + 12*(
2*a^3 + a*b^2)*c^7 - 2*(8*a^4 + 11*a^2*b^2 + b^4)*c^6 + 4*(a^5 + 3*a^3*b^2 + 2*a*b^4)*c^5 - (a^4*b^2 + 2*a^2*b
^4 + b^6)*c^4)))/(4*a*c^5 - (8*a^2 + b^2)*c^4 + 2*(2*a^3 + 3*a*b^2)*c^3 - (a^2*b^2 + b^4)*c^2)) + 4*(2*a^3*b*c
^2 - (a^4*b + a^2*b^3)*c)*cosh(x) + 4*(2*a^3*b*c^2 - (a^4*b + a^2*b^3)*c)*sinh(x) - 2*(4*a^3*c^5 - (8*a^4 + a^
2*b^2)*c^4 + 2*(2*a^5 + 3*a^3*b^2)*c^3 - (a^4*b^2 + a^2*b^4)*c^2)*sqrt(-(a^4*b^2 + 2*a^2*b^4 + b^6 + 4*a^2*b^2
*c^2 - 4*(a^3*b^2 + a*b^4)*c)/(4*a*c^9 - (16*a^2 + b^2)*c^8 + 12*(2*a^3 + a*b^2)*c^7 - 2*(8*a^4 + 11*a^2*b^2 +
b^4)*c^6 + 4*(a^5 + 3*a^3*b^2 + 2*a*b^4)*c^5 - (a^4*b^2 + 2*a^2*b^4 + b^6)*c^4))) - sqrt(2)*c*sqrt(-(a^2*b^2
+ b^4 + 2*a^2*c^2 - 2*(a^3 + 2*a*b^2)*c + (4*a*c^5 - (8*a^2 + b^2)*c^4 + 2*(2*a^3 + 3*a*b^2)*c^3 - (a^2*b^2 +
b^4)*c^2)*sqrt(-(a^4*b^2 + 2*a^2*b^4 + b^6 + 4*a^2*b^2*c^2 - 4*(a^3*b^2 + a*b^4)*c)/(4*a*c^9 - (16*a^2 + b^2)*
c^8 + 12*(2*a^3 + a*b^2)*c^7 - 2*(8*a^4 + 11*a^2*b^2 + b^4)*c^6 + 4*(a^5 + 3*a^3*b^2 + 2*a*b^4)*c^5 - (a^4*b^2
+ 2*a^2*b^4 + b^6)*c^4)))/(4*a*c^5 - (8*a^2 + b^2)*c^4 + 2*(2*a^3 + 3*a*b^2)*c^3 - (a^2*b^2 + b^4)*c^2))*log(
-2*a^4*b^2 - 2*a^2*b^4 + 4*a^3*b^2*c - sqrt(2)*(8*a^2*b^2*c^3 - 2*(2*a^3*b^2 + 3*a*b^4)*c^2 + (a^2*b^4 + b^6)*
c + (8*a^2*c^7 - 6*(4*a^3 + a*b^2)*c^6 + (24*a^4 + 22*a^2*b^2 + b^4)*c^5 - 2*(4*a^5 + 9*a^3*b^2 + 4*a*b^4)*c^4
+ (2*a^4*b^2 + 3*a^2*b^4 + b^6)*c^3)*sqrt(-(a^4*b^2 + 2*a^2*b^4 + b^6 + 4*a^2*b^2*c^2 - 4*(a^3*b^2 + a*b^4)*c
)/(4*a*c^9 - (16*a^2 + b^2)*c^8 + 12*(2*a^3 + a*b^2)*c^7 - 2*(8*a^4 + 11*a^2*b^2 + b^4)*c^6 + 4*(a^5 + 3*a^3*b
^2 + 2*a*b^4)*c^5 - (a^4*b^2 + 2*a^2*b^4 + b^6)*c^4)))*sqrt(-(a^2*b^2 + b^4 + 2*a^2*c^2 - 2*(a^3 + 2*a*b^2)*c
+ (4*a*c^5 - (8*a^2 + b^2)*c^4 + 2*(2*a^3 + 3*a*b^2)*c^3 - (a^2*b^2 + b^4)*c^2)*sqrt(-(a^4*b^2 + 2*a^2*b^4 + b
^6 + 4*a^2*b^2*c^2 - 4*(a^3*b^2 + a*b^4)*c)/(4*a*c^9 - (16*a^2 + b^2)*c^8 + 12*(2*a^3 + a*b^2)*c^7 - 2*(8*a^4
+ 11*a^2*b^2 + b^4)*c^6 + 4*(a^5 + 3*a^3*b^2 + 2*a*b^4)*c^5 - (a^4*b^2 + 2*a^2*b^4 + b^6)*c^4)))/(4*a*c^5 - (8
*a^2 + b^2)*c^4 + 2*(2*a^3 + 3*a*b^2)*c^3 - (a^2*b^2 + b^4)*c^2)) + 4*(2*a^3*b*c^2 - (a^4*b + a^2*b^3)*c)*cosh
(x) + 4*(2*a^3*b*c^2 - (a^4*b + a^2*b^3)*c)*sinh(x) - 2*(4*a^3*c^5 - (8*a^4 + a^2*b^2)*c^4 + 2*(2*a^5 + 3*a^3*
b^2)*c^3 - (a^4*b^2 + a^2*b^4)*c^2)*sqrt(-(a^4*b^2 + 2*a^2*b^4 + b^6 + 4*a^2*b^2*c^2 - 4*(a^3*b^2 + a*b^4)*c)/
(4*a*c^9 - (16*a^2 + b^2)*c^8 + 12*(2*a^3 + a*b^2)*c^7 - 2*(8*a^4 + 11*a^2*b^2 + b^4)*c^6 + 4*(a^5 + 3*a^3*b^2
+ 2*a*b^4)*c^5 - (a^4*b^2 + 2*a^2*b^4 + b^6)*c^4))) + sqrt(2)*c*sqrt(-(a^2*b^2 + b^4 + 2*a^2*c^2 - 2*(a^3 + 2
*a*b^2)*c - (4*a*c^5 - (8*a^2 + b^2)*c^4 + 2*(2*a^3 + 3*a*b^2)*c^3 - (a^2*b^2 + b^4)*c^2)*sqrt(-(a^4*b^2 + 2*a
^2*b^4 + b^6 + 4*a^2*b^2*c^2 - 4*(a^3*b^2 + a*b^4)*c)/(4*a*c^9 - (16*a^2 + b^2)*c^8 + 12*(2*a^3 + a*b^2)*c^7 -
2*(8*a^4 + 11*a^2*b^2 + b^4)*c^6 + 4*(a^5 + 3*a^3*b^2 + 2*a*b^4)*c^5 - (a^4*b^2 + 2*a^2*b^4 + b^6)*c^4)))/(4*
a*c^5 - (8*a^2 + b^2)*c^4 + 2*(2*a^3 + 3*a*b^2)*c^3 - (a^2*b^2 + b^4)*c^2))*log(-2*a^4*b^2 - 2*a^2*b^4 + 4*a^3
*b^2*c + sqrt(2)*(8*a^2*b^2*c^3 - 2*(2*a^3*b^2 + 3*a*b^4)*c^2 + (a^2*b^4 + b^6)*c - (8*a^2*c^7 - 6*(4*a^3 + a*
b^2)*c^6 + (24*a^4 + 22*a^2*b^2 + b^4)*c^5 - 2*(4*a^5 + 9*a^3*b^2 + 4*a*b^4)*c^4 + (2*a^4*b^2 + 3*a^2*b^4 + b^
6)*c^3)*sqrt(-(a^4*b^2 + 2*a^2*b^4 + b^6 + 4*a^2*b^2*c^2 - 4*(a^3*b^2 + a*b^4)*c)/(4*a*c^9 - (16*a^2 + b^2)*c^
8 + 12*(2*a^3 + a*b^2)*c^7 - 2*(8*a^4 + 11*a^2*b^2 + b^4)*c^6 + 4*(a^5 + 3*a^3*b^2 + 2*a*b^4)*c^5 - (a^4*b^2 +
2*a^2*b^4 + b^6)*c^4)))*sqrt(-(a^2*b^2 + b^4 + 2*a^2*c^2 - 2*(a^3 + 2*a*b^2)*c - (4*a*c^5 - (8*a^2 + b^2)*c^4
+ 2*(2*a^3 + 3*a*b^2)*c^3 - (a^2*b^2 + b^4)*c^2)*sqrt(-(a^4*b^2 + 2*a^2*b^4 + b^6 + 4*a^2*b^2*c^2 - 4*(a^3*b^
2 + a*b^4)*c)/(4*a*c^9 - (16*a^2 + b^2)*c^8 + 12*(2*a^3 + a*b^2)*c^7 - 2*(8*a^4 + 11*a^2*b^2 + b^4)*c^6 + 4*(a
^5 + 3*a^3*b^2 + 2*a*b^4)*c^5 - (a^4*b^2 + 2*a^2*b^4 + b^6)*c^4)))/(4*a*c^5 - (8*a^2 + b^2)*c^4 + 2*(2*a^3 + 3
*a*b^2)*c^3 - (a^2*b^2 + b^4)*c^2)) + 4*(2*a^3*b*c^2 - (a^4*b + a^2*b^3)*c)*cosh(x) + 4*(2*a^3*b*c^2 - (a^4*b
+ a^2*b^3)*c)*sinh(x) + 2*(4*a^3*c^5 - (8*a^4 + a^2*b^2)*c^4 + 2*(2*a^5 + 3*a^3*b^2)*c^3 - (a^4*b^2 + a^2*b^4)
*c^2)*sqrt(-(a^4*b^2 + 2*a^2*b^4 + b^6 + 4*a^2*b^2*c^2 - 4*(a^3*b^2 + a*b^4)*c)/(4*a*c^9 - (16*a^2 + b^2)*c^8
+ 12*(2*a^3 + a*b^2)*c^7 - 2*(8*a^4 + 11*a^2*b^2 + b^4)*c^6 + 4*(a^5 + 3*a^3*b^2 + 2*a*b^4)*c^5 - (a^4*b^2 + 2
*a^2*b^4 + b^6)*c^4))) - sqrt(2)*c*sqrt(-(a^2*b^2 + b^4 + 2*a^2*c^2 - 2*(a^3 + 2*a*b^2)*c - (4*a*c^5 - (8*a^2
+ b^2)*c^4 + 2*(2*a^3 + 3*a*b^2)*c^3 - (a^2*b^2 + b^4)*c^2)*sqrt(-(a^4*b^2 + 2*a^2*b^4 + b^6 + 4*a^2*b^2*c^2 -
4*(a^3*b^2 + a*b^4)*c)/(4*a*c^9 - (16*a^2 + b^2)*c^8 + 12*(2*a^3 + a*b^2)*c^7 - 2*(8*a^4 + 11*a^2*b^2 + b^4)*
c^6 + 4*(a^5 + 3*a^3*b^2 + 2*a*b^4)*c^5 - (a^4*b^2 + 2*a^2*b^4 + b^6)*c^4)))/(4*a*c^5 - (8*a^2 + b^2)*c^4 + 2*
(2*a^3 + 3*a*b^2)*c^3 - (a^2*b^2 + b^4)*c^2))*log(-2*a^4*b^2 - 2*a^2*b^4 + 4*a^3*b^2*c - sqrt(2)*(8*a^2*b^2*c^
3 - 2*(2*a^3*b^2 + 3*a*b^4)*c^2 + (a^2*b^4 + b^6)*c - (8*a^2*c^7 - 6*(4*a^3 + a*b^2)*c^6 + (24*a^4 + 22*a^2*b^
2 + b^4)*c^5 - 2*(4*a^5 + 9*a^3*b^2 + 4*a*b^4)*c^4 + (2*a^4*b^2 + 3*a^2*b^4 + b^6)*c^3)*sqrt(-(a^4*b^2 + 2*a^2
*b^4 + b^6 + 4*a^2*b^2*c^2 - 4*(a^3*b^2 + a*b^4)*c)/(4*a*c^9 - (16*a^2 + b^2)*c^8 + 12*(2*a^3 + a*b^2)*c^7 - 2
*(8*a^4 + 11*a^2*b^2 + b^4)*c^6 + 4*(a^5 + 3*a^3*b^2 + 2*a*b^4)*c^5 - (a^4*b^2 + 2*a^2*b^4 + b^6)*c^4)))*sqrt(
-(a^2*b^2 + b^4 + 2*a^2*c^2 - 2*(a^3 + 2*a*b^2)*c - (4*a*c^5 - (8*a^2 + b^2)*c^4 + 2*(2*a^3 + 3*a*b^2)*c^3 - (
a^2*b^2 + b^4)*c^2)*sqrt(-(a^4*b^2 + 2*a^2*b^4 + b^6 + 4*a^2*b^2*c^2 - 4*(a^3*b^2 + a*b^4)*c)/(4*a*c^9 - (16*a
^2 + b^2)*c^8 + 12*(2*a^3 + a*b^2)*c^7 - 2*(8*a^4 + 11*a^2*b^2 + b^4)*c^6 + 4*(a^5 + 3*a^3*b^2 + 2*a*b^4)*c^5
- (a^4*b^2 + 2*a^2*b^4 + b^6)*c^4)))/(4*a*c^5 - (8*a^2 + b^2)*c^4 + 2*(2*a^3 + 3*a*b^2)*c^3 - (a^2*b^2 + b^4)*
c^2)) + 4*(2*a^3*b*c^2 - (a^4*b + a^2*b^3)*c)*cosh(x) + 4*(2*a^3*b*c^2 - (a^4*b + a^2*b^3)*c)*sinh(x) + 2*(4*a
^3*c^5 - (8*a^4 + a^2*b^2)*c^4 + 2*(2*a^5 + 3*a^3*b^2)*c^3 - (a^4*b^2 + a^2*b^4)*c^2)*sqrt(-(a^4*b^2 + 2*a^2*b
^4 + b^6 + 4*a^2*b^2*c^2 - 4*(a^3*b^2 + a*b^4)*c)/(4*a*c^9 - (16*a^2 + b^2)*c^8 + 12*(2*a^3 + a*b^2)*c^7 - 2*(
8*a^4 + 11*a^2*b^2 + b^4)*c^6 + 4*(a^5 + 3*a^3*b^2 + 2*a*b^4)*c^5 - (a^4*b^2 + 2*a^2*b^4 + b^6)*c^4))) + 2*x)/
c

Sympy [F(-1)]

Timed out. \[ \int \frac {\sinh ^2(x)}{a+b \sinh (x)+c \sinh ^2(x)} \, dx=\text {Timed out} \]

[In]

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

[Out]

Timed out

Maxima [F]

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

[In]

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

[Out]

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

Giac [A] (verification not implemented)

none

Time = 0.76 (sec) , antiderivative size = 5, normalized size of antiderivative = 0.02 \[ \int \frac {\sinh ^2(x)}{a+b \sinh (x)+c \sinh ^2(x)} \, dx=\frac {x}{c} \]

[In]

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

[Out]

x/c

Mupad [F(-1)]

Timed out. \[ \int \frac {\sinh ^2(x)}{a+b \sinh (x)+c \sinh ^2(x)} \, dx=\text {Hanged} \]

[In]

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

[Out]

\text{Hanged}