\(\int x^m \cosh ^4(a+b \log (c x^n)) \, dx\) [246]

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

Optimal result

Integrand size = 17, antiderivative size = 266 \[ \int x^m \cosh ^4\left (a+b \log \left (c x^n\right )\right ) \, dx=\frac {24 b^4 n^4 x^{1+m}}{(1+m) \left ((1+m)^2-16 b^2 n^2\right ) \left ((1+m)^2-4 b^2 n^2\right )}-\frac {12 b^2 (1+m) n^2 x^{1+m} \cosh ^2\left (a+b \log \left (c x^n\right )\right )}{\left ((1+m)^2-16 b^2 n^2\right ) \left ((1+m)^2-4 b^2 n^2\right )}+\frac {(1+m) x^{1+m} \cosh ^4\left (a+b \log \left (c x^n\right )\right )}{(1+m)^2-16 b^2 n^2}+\frac {24 b^3 n^3 x^{1+m} \cosh \left (a+b \log \left (c x^n\right )\right ) \sinh \left (a+b \log \left (c x^n\right )\right )}{\left ((1+m)^2-16 b^2 n^2\right ) \left ((1+m)^2-4 b^2 n^2\right )}-\frac {4 b n x^{1+m} \cosh ^3\left (a+b \log \left (c x^n\right )\right ) \sinh \left (a+b \log \left (c x^n\right )\right )}{(1+m)^2-16 b^2 n^2} \]

[Out]

24*b^4*n^4*x^(1+m)/(1+m)/((1+m)^2-16*b^2*n^2)/((1+m)^2-4*b^2*n^2)-12*b^2*(1+m)*n^2*x^(1+m)*cosh(a+b*ln(c*x^n))
^2/((1+m)^2-16*b^2*n^2)/((1+m)^2-4*b^2*n^2)+(1+m)*x^(1+m)*cosh(a+b*ln(c*x^n))^4/((1+m)^2-16*b^2*n^2)+24*b^3*n^
3*x^(1+m)*cosh(a+b*ln(c*x^n))*sinh(a+b*ln(c*x^n))/((1+m)^2-16*b^2*n^2)/((1+m)^2-4*b^2*n^2)-4*b*n*x^(1+m)*cosh(
a+b*ln(c*x^n))^3*sinh(a+b*ln(c*x^n))/((1+m)^2-16*b^2*n^2)

Rubi [A] (verified)

Time = 0.10 (sec) , antiderivative size = 267, normalized size of antiderivative = 1.00, number of steps used = 3, number of rules used = 2, \(\frac {\text {number of rules}}{\text {integrand size}}\) = 0.118, Rules used = {5641, 30} \[ \int x^m \cosh ^4\left (a+b \log \left (c x^n\right )\right ) \, dx=\frac {(m+1) x^{m+1} \cosh ^4\left (a+b \log \left (c x^n\right )\right )}{-16 b^2 n^2+m^2+2 m+1}-\frac {12 b^2 (m+1) n^2 x^{m+1} \cosh ^2\left (a+b \log \left (c x^n\right )\right )}{\left ((m+1)^2-16 b^2 n^2\right ) \left (-4 b^2 n^2+m^2+2 m+1\right )}-\frac {4 b n x^{m+1} \sinh \left (a+b \log \left (c x^n\right )\right ) \cosh ^3\left (a+b \log \left (c x^n\right )\right )}{(m+1)^2-16 b^2 n^2}+\frac {24 b^3 n^3 x^{m+1} \sinh \left (a+b \log \left (c x^n\right )\right ) \cosh \left (a+b \log \left (c x^n\right )\right )}{64 b^4 n^4-20 b^2 (m+1)^2 n^2+(m+1)^4}+\frac {24 b^4 n^4 x^{m+1}}{(m+1) \left ((m+1)^2-16 b^2 n^2\right ) \left ((m+1)^2-4 b^2 n^2\right )} \]

[In]

Int[x^m*Cosh[a + b*Log[c*x^n]]^4,x]

[Out]

(24*b^4*n^4*x^(1 + m))/((1 + m)*((1 + m)^2 - 16*b^2*n^2)*((1 + m)^2 - 4*b^2*n^2)) - (12*b^2*(1 + m)*n^2*x^(1 +
 m)*Cosh[a + b*Log[c*x^n]]^2)/(((1 + m)^2 - 16*b^2*n^2)*(1 + 2*m + m^2 - 4*b^2*n^2)) + ((1 + m)*x^(1 + m)*Cosh
[a + b*Log[c*x^n]]^4)/(1 + 2*m + m^2 - 16*b^2*n^2) + (24*b^3*n^3*x^(1 + m)*Cosh[a + b*Log[c*x^n]]*Sinh[a + b*L
og[c*x^n]])/((1 + m)^4 - 20*b^2*(1 + m)^2*n^2 + 64*b^4*n^4) - (4*b*n*x^(1 + m)*Cosh[a + b*Log[c*x^n]]^3*Sinh[a
 + b*Log[c*x^n]])/((1 + m)^2 - 16*b^2*n^2)

Rule 30

Int[(x_)^(m_.), x_Symbol] :> Simp[x^(m + 1)/(m + 1), x] /; FreeQ[m, x] && NeQ[m, -1]

Rule 5641

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

Rubi steps \begin{align*} \text {integral}& = \frac {(1+m) x^{1+m} \cosh ^4\left (a+b \log \left (c x^n\right )\right )}{1+2 m+m^2-16 b^2 n^2}-\frac {4 b n x^{1+m} \cosh ^3\left (a+b \log \left (c x^n\right )\right ) \sinh \left (a+b \log \left (c x^n\right )\right )}{(1+m)^2-16 b^2 n^2}-\frac {\left (12 b^2 n^2\right ) \int x^m \cosh ^2\left (a+b \log \left (c x^n\right )\right ) \, dx}{(1+m)^2-16 b^2 n^2} \\ & = -\frac {12 b^2 (1+m) n^2 x^{1+m} \cosh ^2\left (a+b \log \left (c x^n\right )\right )}{\left ((1+m)^2-16 b^2 n^2\right ) \left (1+2 m+m^2-4 b^2 n^2\right )}+\frac {(1+m) x^{1+m} \cosh ^4\left (a+b \log \left (c x^n\right )\right )}{1+2 m+m^2-16 b^2 n^2}+\frac {24 b^3 n^3 x^{1+m} \cosh \left (a+b \log \left (c x^n\right )\right ) \sinh \left (a+b \log \left (c x^n\right )\right )}{(1+m)^4-20 b^2 (1+m)^2 n^2+64 b^4 n^4}-\frac {4 b n x^{1+m} \cosh ^3\left (a+b \log \left (c x^n\right )\right ) \sinh \left (a+b \log \left (c x^n\right )\right )}{(1+m)^2-16 b^2 n^2}+\frac {\left (24 b^4 n^4\right ) \int x^m \, dx}{(1+m)^4-20 b^2 (1+m)^2 n^2+64 b^4 n^4} \\ & = \frac {24 b^4 n^4 x^{1+m}}{(1+m) \left ((1+m)^4-20 b^2 (1+m)^2 n^2+64 b^4 n^4\right )}-\frac {12 b^2 (1+m) n^2 x^{1+m} \cosh ^2\left (a+b \log \left (c x^n\right )\right )}{\left ((1+m)^2-16 b^2 n^2\right ) \left (1+2 m+m^2-4 b^2 n^2\right )}+\frac {(1+m) x^{1+m} \cosh ^4\left (a+b \log \left (c x^n\right )\right )}{1+2 m+m^2-16 b^2 n^2}+\frac {24 b^3 n^3 x^{1+m} \cosh \left (a+b \log \left (c x^n\right )\right ) \sinh \left (a+b \log \left (c x^n\right )\right )}{(1+m)^4-20 b^2 (1+m)^2 n^2+64 b^4 n^4}-\frac {4 b n x^{1+m} \cosh ^3\left (a+b \log \left (c x^n\right )\right ) \sinh \left (a+b \log \left (c x^n\right )\right )}{(1+m)^2-16 b^2 n^2} \\ \end{align*}

Mathematica [A] (verified)

Time = 2.52 (sec) , antiderivative size = 311, normalized size of antiderivative = 1.17 \[ \int x^m \cosh ^4\left (a+b \log \left (c x^n\right )\right ) \, dx=\frac {1}{8} x^{1+m} \left (\frac {3}{1+m}+\frac {4 \sinh (2 b n \log (x)) \left (-2 b n \cosh \left (2 \left (a-b n \log (x)+b \log \left (c x^n\right )\right )\right )+(1+m) \sinh \left (2 \left (a-b n \log (x)+b \log \left (c x^n\right )\right )\right )\right )}{(1+m-2 b n) (1+m+2 b n)}+\frac {4 \cosh (2 b n \log (x)) \left ((1+m) \cosh \left (2 \left (a-b n \log (x)+b \log \left (c x^n\right )\right )\right )-2 b n \sinh \left (2 \left (a-b n \log (x)+b \log \left (c x^n\right )\right )\right )\right )}{(1+m-2 b n) (1+m+2 b n)}+\frac {\sinh (4 b n \log (x)) \left (-4 b n \cosh \left (4 \left (a-b n \log (x)+b \log \left (c x^n\right )\right )\right )+(1+m) \sinh \left (4 \left (a-b n \log (x)+b \log \left (c x^n\right )\right )\right )\right )}{(1+m-4 b n) (1+m+4 b n)}+\frac {\cosh (4 b n \log (x)) \left ((1+m) \cosh \left (4 \left (a-b n \log (x)+b \log \left (c x^n\right )\right )\right )-4 b n \sinh \left (4 \left (a-b n \log (x)+b \log \left (c x^n\right )\right )\right )\right )}{(1+m-4 b n) (1+m+4 b n)}\right ) \]

[In]

Integrate[x^m*Cosh[a + b*Log[c*x^n]]^4,x]

[Out]

(x^(1 + m)*(3/(1 + m) + (4*Sinh[2*b*n*Log[x]]*(-2*b*n*Cosh[2*(a - b*n*Log[x] + b*Log[c*x^n])] + (1 + m)*Sinh[2
*(a - b*n*Log[x] + b*Log[c*x^n])]))/((1 + m - 2*b*n)*(1 + m + 2*b*n)) + (4*Cosh[2*b*n*Log[x]]*((1 + m)*Cosh[2*
(a - b*n*Log[x] + b*Log[c*x^n])] - 2*b*n*Sinh[2*(a - b*n*Log[x] + b*Log[c*x^n])]))/((1 + m - 2*b*n)*(1 + m + 2
*b*n)) + (Sinh[4*b*n*Log[x]]*(-4*b*n*Cosh[4*(a - b*n*Log[x] + b*Log[c*x^n])] + (1 + m)*Sinh[4*(a - b*n*Log[x]
+ b*Log[c*x^n])]))/((1 + m - 4*b*n)*(1 + m + 4*b*n)) + (Cosh[4*b*n*Log[x]]*((1 + m)*Cosh[4*(a - b*n*Log[x] + b
*Log[c*x^n])] - 4*b*n*Sinh[4*(a - b*n*Log[x] + b*Log[c*x^n])]))/((1 + m - 4*b*n)*(1 + m + 4*b*n))))/8

Maple [A] (verified)

Time = 236.93 (sec) , antiderivative size = 214, normalized size of antiderivative = 0.80

method result size
parallelrisch \(-\frac {8 x^{1+m} \left (-3 b^{4} n^{4}-2 \left (1+m \right ) n^{3} \left (\sinh \left (2 b \ln \left (c \,x^{n}\right )+2 a \right )+\frac {\sinh \left (4 b \ln \left (c \,x^{n}\right )+4 a \right )}{8}\right ) b^{3}+\left (1+m \right )^{2} n^{2} \left (\cosh \left (2 b \ln \left (c \,x^{n}\right )+2 a \right )+\frac {\cosh \left (4 b \ln \left (c \,x^{n}\right )+4 a \right )}{16}+\frac {15}{16}\right ) b^{2}+\frac {\left (1+m \right )^{3} \left (\sinh \left (2 b \ln \left (c \,x^{n}\right )+2 a \right )+\frac {\sinh \left (4 b \ln \left (c \,x^{n}\right )+4 a \right )}{2}\right ) n b}{8}-\frac {\left (1+m \right )^{4} \left (\cosh \left (2 b \ln \left (c \,x^{n}\right )+2 a \right )+\frac {\cosh \left (4 b \ln \left (c \,x^{n}\right )+4 a \right )}{4}+\frac {3}{4}\right )}{16}\right )}{64 n^{4} \left (1+m \right ) b^{4}-20 n^{2} \left (1+m \right )^{3} b^{2}+\left (1+m \right )^{5}}\) \(214\)

[In]

int(x^m*cosh(a+b*ln(c*x^n))^4,x,method=_RETURNVERBOSE)

[Out]

-8*x^(1+m)*(-3*b^4*n^4-2*(1+m)*n^3*(sinh(2*b*ln(c*x^n)+2*a)+1/8*sinh(4*b*ln(c*x^n)+4*a))*b^3+(1+m)^2*n^2*(cosh
(2*b*ln(c*x^n)+2*a)+1/16*cosh(4*b*ln(c*x^n)+4*a)+15/16)*b^2+1/8*(1+m)^3*(sinh(2*b*ln(c*x^n)+2*a)+1/2*sinh(4*b*
ln(c*x^n)+4*a))*n*b-1/16*(1+m)^4*(cosh(2*b*ln(c*x^n)+2*a)+1/4*cosh(4*b*ln(c*x^n)+4*a)+3/4))/(64*n^4*(1+m)*b^4-
20*n^2*(1+m)^3*b^2+(1+m)^5)

Fricas [B] (verification not implemented)

Leaf count of result is larger than twice the leaf count of optimal. 1123 vs. \(2 (283) = 566\).

Time = 0.32 (sec) , antiderivative size = 1123, normalized size of antiderivative = 4.22 \[ \int x^m \cosh ^4\left (a+b \log \left (c x^n\right )\right ) \, dx=\text {Too large to display} \]

[In]

integrate(x^m*cosh(a+b*log(c*x^n))^4,x, algorithm="fricas")

[Out]

1/8*((m^4 + 4*m^3 - 4*(b^2*m^2 + 2*b^2*m + b^2)*n^2 + 6*m^2 + 4*m + 1)*x*cosh(b*n*log(x) + b*log(c) + a)^4*cos
h(m*log(x)) + 4*(m^4 + 4*m^3 - 16*(b^2*m^2 + 2*b^2*m + b^2)*n^2 + 6*m^2 + 4*m + 1)*x*cosh(b*n*log(x) + b*log(c
) + a)^2*cosh(m*log(x)) + ((m^4 + 4*m^3 - 4*(b^2*m^2 + 2*b^2*m + b^2)*n^2 + 6*m^2 + 4*m + 1)*x*cosh(m*log(x))
+ (m^4 + 4*m^3 - 4*(b^2*m^2 + 2*b^2*m + b^2)*n^2 + 6*m^2 + 4*m + 1)*x*sinh(m*log(x)))*sinh(b*n*log(x) + b*log(
c) + a)^4 + 16*((4*(b^3*m + b^3)*n^3 - (b*m^3 + 3*b*m^2 + 3*b*m + b)*n)*x*cosh(b*n*log(x) + b*log(c) + a)*cosh
(m*log(x)) + (4*(b^3*m + b^3)*n^3 - (b*m^3 + 3*b*m^2 + 3*b*m + b)*n)*x*cosh(b*n*log(x) + b*log(c) + a)*sinh(m*
log(x)))*sinh(b*n*log(x) + b*log(c) + a)^3 + 3*(64*b^4*n^4 + m^4 + 4*m^3 - 20*(b^2*m^2 + 2*b^2*m + b^2)*n^2 +
6*m^2 + 4*m + 1)*x*cosh(m*log(x)) + 2*(3*(m^4 + 4*m^3 - 4*(b^2*m^2 + 2*b^2*m + b^2)*n^2 + 6*m^2 + 4*m + 1)*x*c
osh(b*n*log(x) + b*log(c) + a)^2*cosh(m*log(x)) + 2*(m^4 + 4*m^3 - 16*(b^2*m^2 + 2*b^2*m + b^2)*n^2 + 6*m^2 +
4*m + 1)*x*cosh(m*log(x)) + (3*(m^4 + 4*m^3 - 4*(b^2*m^2 + 2*b^2*m + b^2)*n^2 + 6*m^2 + 4*m + 1)*x*cosh(b*n*lo
g(x) + b*log(c) + a)^2 + 2*(m^4 + 4*m^3 - 16*(b^2*m^2 + 2*b^2*m + b^2)*n^2 + 6*m^2 + 4*m + 1)*x)*sinh(m*log(x)
))*sinh(b*n*log(x) + b*log(c) + a)^2 + 16*((4*(b^3*m + b^3)*n^3 - (b*m^3 + 3*b*m^2 + 3*b*m + b)*n)*x*cosh(b*n*
log(x) + b*log(c) + a)^3*cosh(m*log(x)) + (16*(b^3*m + b^3)*n^3 - (b*m^3 + 3*b*m^2 + 3*b*m + b)*n)*x*cosh(b*n*
log(x) + b*log(c) + a)*cosh(m*log(x)) + ((4*(b^3*m + b^3)*n^3 - (b*m^3 + 3*b*m^2 + 3*b*m + b)*n)*x*cosh(b*n*lo
g(x) + b*log(c) + a)^3 + (16*(b^3*m + b^3)*n^3 - (b*m^3 + 3*b*m^2 + 3*b*m + b)*n)*x*cosh(b*n*log(x) + b*log(c)
 + a))*sinh(m*log(x)))*sinh(b*n*log(x) + b*log(c) + a) + ((m^4 + 4*m^3 - 4*(b^2*m^2 + 2*b^2*m + b^2)*n^2 + 6*m
^2 + 4*m + 1)*x*cosh(b*n*log(x) + b*log(c) + a)^4 + 4*(m^4 + 4*m^3 - 16*(b^2*m^2 + 2*b^2*m + b^2)*n^2 + 6*m^2
+ 4*m + 1)*x*cosh(b*n*log(x) + b*log(c) + a)^2 + 3*(64*b^4*n^4 + m^4 + 4*m^3 - 20*(b^2*m^2 + 2*b^2*m + b^2)*n^
2 + 6*m^2 + 4*m + 1)*x)*sinh(m*log(x)))/(m^5 + 64*(b^4*m + b^4)*n^4 + 5*m^4 + 10*m^3 - 20*(b^2*m^3 + 3*b^2*m^2
 + 3*b^2*m + b^2)*n^2 + 10*m^2 + 5*m + 1)

Sympy [F]

\[ \int x^m \cosh ^4\left (a+b \log \left (c x^n\right )\right ) \, dx=\text {Too large to display} \]

[In]

integrate(x**m*cosh(a+b*ln(c*x**n))**4,x)

[Out]

Piecewise((log(x)*cosh(a)**4, Eq(b, 0) & Eq(m, -1)), (Integral(x**m*cosh(-a + m*log(c*x**n)/(4*n) + log(c*x**n
)/(4*n))**4, x), Eq(b, (-m - 1)/(4*n))), (Integral(x**m*cosh(-a + m*log(c*x**n)/(2*n) + log(c*x**n)/(2*n))**4,
 x), Eq(b, (-m - 1)/(2*n))), (Integral(x**m*cosh(a + m*log(c*x**n)/(4*n) + log(c*x**n)/(4*n))**4, x), Eq(b, (m
 + 1)/(4*n))), (Integral(x**m*cosh(a + m*log(c*x**n)/(2*n) + log(c*x**n)/(2*n))**4, x), Eq(b, (m + 1)/(2*n))),
 (Integral(cosh(a + b*log(c*x**n))**4/x, x), Eq(m, -1)), (24*b**4*n**4*x*x**m*sinh(a + b*log(c*x**n))**4/(64*b
**4*m*n**4 + 64*b**4*n**4 - 20*b**2*m**3*n**2 - 60*b**2*m**2*n**2 - 60*b**2*m*n**2 - 20*b**2*n**2 + m**5 + 5*m
**4 + 10*m**3 + 10*m**2 + 5*m + 1) - 48*b**4*n**4*x*x**m*sinh(a + b*log(c*x**n))**2*cosh(a + b*log(c*x**n))**2
/(64*b**4*m*n**4 + 64*b**4*n**4 - 20*b**2*m**3*n**2 - 60*b**2*m**2*n**2 - 60*b**2*m*n**2 - 20*b**2*n**2 + m**5
 + 5*m**4 + 10*m**3 + 10*m**2 + 5*m + 1) + 24*b**4*n**4*x*x**m*cosh(a + b*log(c*x**n))**4/(64*b**4*m*n**4 + 64
*b**4*n**4 - 20*b**2*m**3*n**2 - 60*b**2*m**2*n**2 - 60*b**2*m*n**2 - 20*b**2*n**2 + m**5 + 5*m**4 + 10*m**3 +
 10*m**2 + 5*m + 1) - 24*b**3*m*n**3*x*x**m*sinh(a + b*log(c*x**n))**3*cosh(a + b*log(c*x**n))/(64*b**4*m*n**4
 + 64*b**4*n**4 - 20*b**2*m**3*n**2 - 60*b**2*m**2*n**2 - 60*b**2*m*n**2 - 20*b**2*n**2 + m**5 + 5*m**4 + 10*m
**3 + 10*m**2 + 5*m + 1) + 40*b**3*m*n**3*x*x**m*sinh(a + b*log(c*x**n))*cosh(a + b*log(c*x**n))**3/(64*b**4*m
*n**4 + 64*b**4*n**4 - 20*b**2*m**3*n**2 - 60*b**2*m**2*n**2 - 60*b**2*m*n**2 - 20*b**2*n**2 + m**5 + 5*m**4 +
 10*m**3 + 10*m**2 + 5*m + 1) - 24*b**3*n**3*x*x**m*sinh(a + b*log(c*x**n))**3*cosh(a + b*log(c*x**n))/(64*b**
4*m*n**4 + 64*b**4*n**4 - 20*b**2*m**3*n**2 - 60*b**2*m**2*n**2 - 60*b**2*m*n**2 - 20*b**2*n**2 + m**5 + 5*m**
4 + 10*m**3 + 10*m**2 + 5*m + 1) + 40*b**3*n**3*x*x**m*sinh(a + b*log(c*x**n))*cosh(a + b*log(c*x**n))**3/(64*
b**4*m*n**4 + 64*b**4*n**4 - 20*b**2*m**3*n**2 - 60*b**2*m**2*n**2 - 60*b**2*m*n**2 - 20*b**2*n**2 + m**5 + 5*
m**4 + 10*m**3 + 10*m**2 + 5*m + 1) + 12*b**2*m**2*n**2*x*x**m*sinh(a + b*log(c*x**n))**2*cosh(a + b*log(c*x**
n))**2/(64*b**4*m*n**4 + 64*b**4*n**4 - 20*b**2*m**3*n**2 - 60*b**2*m**2*n**2 - 60*b**2*m*n**2 - 20*b**2*n**2
+ m**5 + 5*m**4 + 10*m**3 + 10*m**2 + 5*m + 1) - 16*b**2*m**2*n**2*x*x**m*cosh(a + b*log(c*x**n))**4/(64*b**4*
m*n**4 + 64*b**4*n**4 - 20*b**2*m**3*n**2 - 60*b**2*m**2*n**2 - 60*b**2*m*n**2 - 20*b**2*n**2 + m**5 + 5*m**4
+ 10*m**3 + 10*m**2 + 5*m + 1) + 24*b**2*m*n**2*x*x**m*sinh(a + b*log(c*x**n))**2*cosh(a + b*log(c*x**n))**2/(
64*b**4*m*n**4 + 64*b**4*n**4 - 20*b**2*m**3*n**2 - 60*b**2*m**2*n**2 - 60*b**2*m*n**2 - 20*b**2*n**2 + m**5 +
 5*m**4 + 10*m**3 + 10*m**2 + 5*m + 1) - 32*b**2*m*n**2*x*x**m*cosh(a + b*log(c*x**n))**4/(64*b**4*m*n**4 + 64
*b**4*n**4 - 20*b**2*m**3*n**2 - 60*b**2*m**2*n**2 - 60*b**2*m*n**2 - 20*b**2*n**2 + m**5 + 5*m**4 + 10*m**3 +
 10*m**2 + 5*m + 1) + 12*b**2*n**2*x*x**m*sinh(a + b*log(c*x**n))**2*cosh(a + b*log(c*x**n))**2/(64*b**4*m*n**
4 + 64*b**4*n**4 - 20*b**2*m**3*n**2 - 60*b**2*m**2*n**2 - 60*b**2*m*n**2 - 20*b**2*n**2 + m**5 + 5*m**4 + 10*
m**3 + 10*m**2 + 5*m + 1) - 16*b**2*n**2*x*x**m*cosh(a + b*log(c*x**n))**4/(64*b**4*m*n**4 + 64*b**4*n**4 - 20
*b**2*m**3*n**2 - 60*b**2*m**2*n**2 - 60*b**2*m*n**2 - 20*b**2*n**2 + m**5 + 5*m**4 + 10*m**3 + 10*m**2 + 5*m
+ 1) - 4*b*m**3*n*x*x**m*sinh(a + b*log(c*x**n))*cosh(a + b*log(c*x**n))**3/(64*b**4*m*n**4 + 64*b**4*n**4 - 2
0*b**2*m**3*n**2 - 60*b**2*m**2*n**2 - 60*b**2*m*n**2 - 20*b**2*n**2 + m**5 + 5*m**4 + 10*m**3 + 10*m**2 + 5*m
 + 1) - 12*b*m**2*n*x*x**m*sinh(a + b*log(c*x**n))*cosh(a + b*log(c*x**n))**3/(64*b**4*m*n**4 + 64*b**4*n**4 -
 20*b**2*m**3*n**2 - 60*b**2*m**2*n**2 - 60*b**2*m*n**2 - 20*b**2*n**2 + m**5 + 5*m**4 + 10*m**3 + 10*m**2 + 5
*m + 1) - 12*b*m*n*x*x**m*sinh(a + b*log(c*x**n))*cosh(a + b*log(c*x**n))**3/(64*b**4*m*n**4 + 64*b**4*n**4 -
20*b**2*m**3*n**2 - 60*b**2*m**2*n**2 - 60*b**2*m*n**2 - 20*b**2*n**2 + m**5 + 5*m**4 + 10*m**3 + 10*m**2 + 5*
m + 1) - 4*b*n*x*x**m*sinh(a + b*log(c*x**n))*cosh(a + b*log(c*x**n))**3/(64*b**4*m*n**4 + 64*b**4*n**4 - 20*b
**2*m**3*n**2 - 60*b**2*m**2*n**2 - 60*b**2*m*n**2 - 20*b**2*n**2 + m**5 + 5*m**4 + 10*m**3 + 10*m**2 + 5*m +
1) + m**4*x*x**m*cosh(a + b*log(c*x**n))**4/(64*b**4*m*n**4 + 64*b**4*n**4 - 20*b**2*m**3*n**2 - 60*b**2*m**2*
n**2 - 60*b**2*m*n**2 - 20*b**2*n**2 + m**5 + 5*m**4 + 10*m**3 + 10*m**2 + 5*m + 1) + 4*m**3*x*x**m*cosh(a + b
*log(c*x**n))**4/(64*b**4*m*n**4 + 64*b**4*n**4 - 20*b**2*m**3*n**2 - 60*b**2*m**2*n**2 - 60*b**2*m*n**2 - 20*
b**2*n**2 + m**5 + 5*m**4 + 10*m**3 + 10*m**2 + 5*m + 1) + 6*m**2*x*x**m*cosh(a + b*log(c*x**n))**4/(64*b**4*m
*n**4 + 64*b**4*n**4 - 20*b**2*m**3*n**2 - 60*b**2*m**2*n**2 - 60*b**2*m*n**2 - 20*b**2*n**2 + m**5 + 5*m**4 +
 10*m**3 + 10*m**2 + 5*m + 1) + 4*m*x*x**m*cosh(a + b*log(c*x**n))**4/(64*b**4*m*n**4 + 64*b**4*n**4 - 20*b**2
*m**3*n**2 - 60*b**2*m**2*n**2 - 60*b**2*m*n**2 - 20*b**2*n**2 + m**5 + 5*m**4 + 10*m**3 + 10*m**2 + 5*m + 1)
+ x*x**m*cosh(a + b*log(c*x**n))**4/(64*b**4*m*n**4 + 64*b**4*n**4 - 20*b**2*m**3*n**2 - 60*b**2*m**2*n**2 - 6
0*b**2*m*n**2 - 20*b**2*n**2 + m**5 + 5*m**4 + 10*m**3 + 10*m**2 + 5*m + 1), True))

Maxima [A] (verification not implemented)

none

Time = 0.23 (sec) , antiderivative size = 161, normalized size of antiderivative = 0.61 \[ \int x^m \cosh ^4\left (a+b \log \left (c x^n\right )\right ) \, dx=\frac {c^{4 \, b} x e^{\left (4 \, b \log \left (x^{n}\right ) + m \log \left (x\right ) + 4 \, a\right )}}{16 \, {\left (4 \, b n + m + 1\right )}} + \frac {c^{2 \, b} x e^{\left (2 \, b \log \left (x^{n}\right ) + m \log \left (x\right ) + 2 \, a\right )}}{4 \, {\left (2 \, b n + m + 1\right )}} - \frac {x e^{\left (-2 \, b \log \left (x^{n}\right ) + m \log \left (x\right ) - 2 \, a\right )}}{4 \, {\left (2 \, b c^{2 \, b} n - c^{2 \, b} {\left (m + 1\right )}\right )}} - \frac {x e^{\left (-4 \, b \log \left (x^{n}\right ) + m \log \left (x\right ) - 4 \, a\right )}}{16 \, {\left (4 \, b c^{4 \, b} n - c^{4 \, b} {\left (m + 1\right )}\right )}} + \frac {3 \, x^{m + 1}}{8 \, {\left (m + 1\right )}} \]

[In]

integrate(x^m*cosh(a+b*log(c*x^n))^4,x, algorithm="maxima")

[Out]

1/16*c^(4*b)*x*e^(4*b*log(x^n) + m*log(x) + 4*a)/(4*b*n + m + 1) + 1/4*c^(2*b)*x*e^(2*b*log(x^n) + m*log(x) +
2*a)/(2*b*n + m + 1) - 1/4*x*e^(-2*b*log(x^n) + m*log(x) - 2*a)/(2*b*c^(2*b)*n - c^(2*b)*(m + 1)) - 1/16*x*e^(
-4*b*log(x^n) + m*log(x) - 4*a)/(4*b*c^(4*b)*n - c^(4*b)*(m + 1)) + 3/8*x^(m + 1)/(m + 1)

Giac [B] (verification not implemented)

Leaf count of result is larger than twice the leaf count of optimal. 6880 vs. \(2 (283) = 566\).

Time = 0.40 (sec) , antiderivative size = 6880, normalized size of antiderivative = 25.86 \[ \int x^m \cosh ^4\left (a+b \log \left (c x^n\right )\right ) \, dx=\text {Too large to display} \]

[In]

integrate(x^m*cosh(a+b*log(c*x^n))^4,x, algorithm="giac")

[Out]

b^3*c^(4*b)*m*n^3*x*x^(4*b*n)*x^m*e^(4*a)/(64*b^4*m*n^4 + 64*b^4*n^4 - 20*b^2*m^3*n^2 - 60*b^2*m^2*n^2 + m^5 -
 60*b^2*m*n^2 + 5*m^4 - 20*b^2*n^2 + 10*m^3 + 10*m^2 + 5*m + 1) + 8*b^3*c^(2*b)*m*n^3*x*x^(2*b*n)*x^m*e^(2*a)/
(64*b^4*m*n^4 + 64*b^4*n^4 - 20*b^2*m^3*n^2 - 60*b^2*m^2*n^2 + m^5 - 60*b^2*m*n^2 + 5*m^4 - 20*b^2*n^2 + 10*m^
3 + 10*m^2 + 5*m + 1) - 1/4*b^2*c^(4*b)*m^2*n^2*x*x^(4*b*n)*x^m*e^(4*a)/(64*b^4*m*n^4 + 64*b^4*n^4 - 20*b^2*m^
3*n^2 - 60*b^2*m^2*n^2 + m^5 - 60*b^2*m*n^2 + 5*m^4 - 20*b^2*n^2 + 10*m^3 + 10*m^2 + 5*m + 1) + b^3*c^(4*b)*n^
3*x*x^(4*b*n)*x^m*e^(4*a)/(64*b^4*m*n^4 + 64*b^4*n^4 - 20*b^2*m^3*n^2 - 60*b^2*m^2*n^2 + m^5 - 60*b^2*m*n^2 +
5*m^4 - 20*b^2*n^2 + 10*m^3 + 10*m^2 + 5*m + 1) - 4*b^2*c^(2*b)*m^2*n^2*x*x^(2*b*n)*x^m*e^(2*a)/(64*b^4*m*n^4
+ 64*b^4*n^4 - 20*b^2*m^3*n^2 - 60*b^2*m^2*n^2 + m^5 - 60*b^2*m*n^2 + 5*m^4 - 20*b^2*n^2 + 10*m^3 + 10*m^2 + 5
*m + 1) + 8*b^3*c^(2*b)*n^3*x*x^(2*b*n)*x^m*e^(2*a)/(64*b^4*m*n^4 + 64*b^4*n^4 - 20*b^2*m^3*n^2 - 60*b^2*m^2*n
^2 + m^5 - 60*b^2*m*n^2 + 5*m^4 - 20*b^2*n^2 + 10*m^3 + 10*m^2 + 5*m + 1) + 24*b^4*n^4*x*x^m/(64*b^4*m*n^4 + 6
4*b^4*n^4 - 20*b^2*m^3*n^2 - 60*b^2*m^2*n^2 + m^5 - 60*b^2*m*n^2 + 5*m^4 - 20*b^2*n^2 + 10*m^3 + 10*m^2 + 5*m
+ 1) - 1/4*b*c^(4*b)*m^3*n*x*x^(4*b*n)*x^m*e^(4*a)/(64*b^4*m*n^4 + 64*b^4*n^4 - 20*b^2*m^3*n^2 - 60*b^2*m^2*n^
2 + m^5 - 60*b^2*m*n^2 + 5*m^4 - 20*b^2*n^2 + 10*m^3 + 10*m^2 + 5*m + 1) - 1/2*b^2*c^(4*b)*m*n^2*x*x^(4*b*n)*x
^m*e^(4*a)/(64*b^4*m*n^4 + 64*b^4*n^4 - 20*b^2*m^3*n^2 - 60*b^2*m^2*n^2 + m^5 - 60*b^2*m*n^2 + 5*m^4 - 20*b^2*
n^2 + 10*m^3 + 10*m^2 + 5*m + 1) - 1/2*b*c^(2*b)*m^3*n*x*x^(2*b*n)*x^m*e^(2*a)/(64*b^4*m*n^4 + 64*b^4*n^4 - 20
*b^2*m^3*n^2 - 60*b^2*m^2*n^2 + m^5 - 60*b^2*m*n^2 + 5*m^4 - 20*b^2*n^2 + 10*m^3 + 10*m^2 + 5*m + 1) - 8*b^2*c
^(2*b)*m*n^2*x*x^(2*b*n)*x^m*e^(2*a)/(64*b^4*m*n^4 + 64*b^4*n^4 - 20*b^2*m^3*n^2 - 60*b^2*m^2*n^2 + m^5 - 60*b
^2*m*n^2 + 5*m^4 - 20*b^2*n^2 + 10*m^3 + 10*m^2 + 5*m + 1) + 1/16*c^(4*b)*m^4*x*x^(4*b*n)*x^m*e^(4*a)/(64*b^4*
m*n^4 + 64*b^4*n^4 - 20*b^2*m^3*n^2 - 60*b^2*m^2*n^2 + m^5 - 60*b^2*m*n^2 + 5*m^4 - 20*b^2*n^2 + 10*m^3 + 10*m
^2 + 5*m + 1) - 3/4*b*c^(4*b)*m^2*n*x*x^(4*b*n)*x^m*e^(4*a)/(64*b^4*m*n^4 + 64*b^4*n^4 - 20*b^2*m^3*n^2 - 60*b
^2*m^2*n^2 + m^5 - 60*b^2*m*n^2 + 5*m^4 - 20*b^2*n^2 + 10*m^3 + 10*m^2 + 5*m + 1) - 1/4*b^2*c^(4*b)*n^2*x*x^(4
*b*n)*x^m*e^(4*a)/(64*b^4*m*n^4 + 64*b^4*n^4 - 20*b^2*m^3*n^2 - 60*b^2*m^2*n^2 + m^5 - 60*b^2*m*n^2 + 5*m^4 -
20*b^2*n^2 + 10*m^3 + 10*m^2 + 5*m + 1) + 1/4*c^(2*b)*m^4*x*x^(2*b*n)*x^m*e^(2*a)/(64*b^4*m*n^4 + 64*b^4*n^4 -
 20*b^2*m^3*n^2 - 60*b^2*m^2*n^2 + m^5 - 60*b^2*m*n^2 + 5*m^4 - 20*b^2*n^2 + 10*m^3 + 10*m^2 + 5*m + 1) - 3/2*
b*c^(2*b)*m^2*n*x*x^(2*b*n)*x^m*e^(2*a)/(64*b^4*m*n^4 + 64*b^4*n^4 - 20*b^2*m^3*n^2 - 60*b^2*m^2*n^2 + m^5 - 6
0*b^2*m*n^2 + 5*m^4 - 20*b^2*n^2 + 10*m^3 + 10*m^2 + 5*m + 1) - 4*b^2*c^(2*b)*n^2*x*x^(2*b*n)*x^m*e^(2*a)/(64*
b^4*m*n^4 + 64*b^4*n^4 - 20*b^2*m^3*n^2 - 60*b^2*m^2*n^2 + m^5 - 60*b^2*m*n^2 + 5*m^4 - 20*b^2*n^2 + 10*m^3 +
10*m^2 + 5*m + 1) - 15/2*b^2*m^2*n^2*x*x^m/(64*b^4*m*n^4 + 64*b^4*n^4 - 20*b^2*m^3*n^2 - 60*b^2*m^2*n^2 + m^5
- 60*b^2*m*n^2 + 5*m^4 - 20*b^2*n^2 + 10*m^3 + 10*m^2 + 5*m + 1) + 1/4*c^(4*b)*m^3*x*x^(4*b*n)*x^m*e^(4*a)/(64
*b^4*m*n^4 + 64*b^4*n^4 - 20*b^2*m^3*n^2 - 60*b^2*m^2*n^2 + m^5 - 60*b^2*m*n^2 + 5*m^4 - 20*b^2*n^2 + 10*m^3 +
 10*m^2 + 5*m + 1) - 3/4*b*c^(4*b)*m*n*x*x^(4*b*n)*x^m*e^(4*a)/(64*b^4*m*n^4 + 64*b^4*n^4 - 20*b^2*m^3*n^2 - 6
0*b^2*m^2*n^2 + m^5 - 60*b^2*m*n^2 + 5*m^4 - 20*b^2*n^2 + 10*m^3 + 10*m^2 + 5*m + 1) + c^(2*b)*m^3*x*x^(2*b*n)
*x^m*e^(2*a)/(64*b^4*m*n^4 + 64*b^4*n^4 - 20*b^2*m^3*n^2 - 60*b^2*m^2*n^2 + m^5 - 60*b^2*m*n^2 + 5*m^4 - 20*b^
2*n^2 + 10*m^3 + 10*m^2 + 5*m + 1) - 3/2*b*c^(2*b)*m*n*x*x^(2*b*n)*x^m*e^(2*a)/(64*b^4*m*n^4 + 64*b^4*n^4 - 20
*b^2*m^3*n^2 - 60*b^2*m^2*n^2 + m^5 - 60*b^2*m*n^2 + 5*m^4 - 20*b^2*n^2 + 10*m^3 + 10*m^2 + 5*m + 1) - 8*b^3*m
*n^3*x*x^m*e^(-2*a)/((64*b^4*m*n^4 + 64*b^4*n^4 - 20*b^2*m^3*n^2 - 60*b^2*m^2*n^2 + m^5 - 60*b^2*m*n^2 + 5*m^4
 - 20*b^2*n^2 + 10*m^3 + 10*m^2 + 5*m + 1)*c^(2*b)*x^(2*b*n)) - b^3*m*n^3*x*x^m*e^(-4*a)/((64*b^4*m*n^4 + 64*b
^4*n^4 - 20*b^2*m^3*n^2 - 60*b^2*m^2*n^2 + m^5 - 60*b^2*m*n^2 + 5*m^4 - 20*b^2*n^2 + 10*m^3 + 10*m^2 + 5*m + 1
)*c^(4*b)*x^(4*b*n)) - 15*b^2*m*n^2*x*x^m/(64*b^4*m*n^4 + 64*b^4*n^4 - 20*b^2*m^3*n^2 - 60*b^2*m^2*n^2 + m^5 -
 60*b^2*m*n^2 + 5*m^4 - 20*b^2*n^2 + 10*m^3 + 10*m^2 + 5*m + 1) + 3/8*c^(4*b)*m^2*x*x^(4*b*n)*x^m*e^(4*a)/(64*
b^4*m*n^4 + 64*b^4*n^4 - 20*b^2*m^3*n^2 - 60*b^2*m^2*n^2 + m^5 - 60*b^2*m*n^2 + 5*m^4 - 20*b^2*n^2 + 10*m^3 +
10*m^2 + 5*m + 1) - 1/4*b*c^(4*b)*n*x*x^(4*b*n)*x^m*e^(4*a)/(64*b^4*m*n^4 + 64*b^4*n^4 - 20*b^2*m^3*n^2 - 60*b
^2*m^2*n^2 + m^5 - 60*b^2*m*n^2 + 5*m^4 - 20*b^2*n^2 + 10*m^3 + 10*m^2 + 5*m + 1) + 3/2*c^(2*b)*m^2*x*x^(2*b*n
)*x^m*e^(2*a)/(64*b^4*m*n^4 + 64*b^4*n^4 - 20*b^2*m^3*n^2 - 60*b^2*m^2*n^2 + m^5 - 60*b^2*m*n^2 + 5*m^4 - 20*b
^2*n^2 + 10*m^3 + 10*m^2 + 5*m + 1) - 1/2*b*c^(2*b)*n*x*x^(2*b*n)*x^m*e^(2*a)/(64*b^4*m*n^4 + 64*b^4*n^4 - 20*
b^2*m^3*n^2 - 60*b^2*m^2*n^2 + m^5 - 60*b^2*m*n^2 + 5*m^4 - 20*b^2*n^2 + 10*m^3 + 10*m^2 + 5*m + 1) - 4*b^2*m^
2*n^2*x*x^m*e^(-2*a)/((64*b^4*m*n^4 + 64*b^4*n^4 - 20*b^2*m^3*n^2 - 60*b^2*m^2*n^2 + m^5 - 60*b^2*m*n^2 + 5*m^
4 - 20*b^2*n^2 + 10*m^3 + 10*m^2 + 5*m + 1)*c^(2*b)*x^(2*b*n)) - 8*b^3*n^3*x*x^m*e^(-2*a)/((64*b^4*m*n^4 + 64*
b^4*n^4 - 20*b^2*m^3*n^2 - 60*b^2*m^2*n^2 + m^5 - 60*b^2*m*n^2 + 5*m^4 - 20*b^2*n^2 + 10*m^3 + 10*m^2 + 5*m +
1)*c^(2*b)*x^(2*b*n)) - 1/4*b^2*m^2*n^2*x*x^m*e^(-4*a)/((64*b^4*m*n^4 + 64*b^4*n^4 - 20*b^2*m^3*n^2 - 60*b^2*m
^2*n^2 + m^5 - 60*b^2*m*n^2 + 5*m^4 - 20*b^2*n^2 + 10*m^3 + 10*m^2 + 5*m + 1)*c^(4*b)*x^(4*b*n)) - b^3*n^3*x*x
^m*e^(-4*a)/((64*b^4*m*n^4 + 64*b^4*n^4 - 20*b^2*m^3*n^2 - 60*b^2*m^2*n^2 + m^5 - 60*b^2*m*n^2 + 5*m^4 - 20*b^
2*n^2 + 10*m^3 + 10*m^2 + 5*m + 1)*c^(4*b)*x^(4*b*n)) + 3/8*m^4*x*x^m/(64*b^4*m*n^4 + 64*b^4*n^4 - 20*b^2*m^3*
n^2 - 60*b^2*m^2*n^2 + m^5 - 60*b^2*m*n^2 + 5*m^4 - 20*b^2*n^2 + 10*m^3 + 10*m^2 + 5*m + 1) - 15/2*b^2*n^2*x*x
^m/(64*b^4*m*n^4 + 64*b^4*n^4 - 20*b^2*m^3*n^2 - 60*b^2*m^2*n^2 + m^5 - 60*b^2*m*n^2 + 5*m^4 - 20*b^2*n^2 + 10
*m^3 + 10*m^2 + 5*m + 1) + 1/4*c^(4*b)*m*x*x^(4*b*n)*x^m*e^(4*a)/(64*b^4*m*n^4 + 64*b^4*n^4 - 20*b^2*m^3*n^2 -
 60*b^2*m^2*n^2 + m^5 - 60*b^2*m*n^2 + 5*m^4 - 20*b^2*n^2 + 10*m^3 + 10*m^2 + 5*m + 1) + c^(2*b)*m*x*x^(2*b*n)
*x^m*e^(2*a)/(64*b^4*m*n^4 + 64*b^4*n^4 - 20*b^2*m^3*n^2 - 60*b^2*m^2*n^2 + m^5 - 60*b^2*m*n^2 + 5*m^4 - 20*b^
2*n^2 + 10*m^3 + 10*m^2 + 5*m + 1) + 1/2*b*m^3*n*x*x^m*e^(-2*a)/((64*b^4*m*n^4 + 64*b^4*n^4 - 20*b^2*m^3*n^2 -
 60*b^2*m^2*n^2 + m^5 - 60*b^2*m*n^2 + 5*m^4 - 20*b^2*n^2 + 10*m^3 + 10*m^2 + 5*m + 1)*c^(2*b)*x^(2*b*n)) - 8*
b^2*m*n^2*x*x^m*e^(-2*a)/((64*b^4*m*n^4 + 64*b^4*n^4 - 20*b^2*m^3*n^2 - 60*b^2*m^2*n^2 + m^5 - 60*b^2*m*n^2 +
5*m^4 - 20*b^2*n^2 + 10*m^3 + 10*m^2 + 5*m + 1)*c^(2*b)*x^(2*b*n)) + 1/4*b*m^3*n*x*x^m*e^(-4*a)/((64*b^4*m*n^4
 + 64*b^4*n^4 - 20*b^2*m^3*n^2 - 60*b^2*m^2*n^2 + m^5 - 60*b^2*m*n^2 + 5*m^4 - 20*b^2*n^2 + 10*m^3 + 10*m^2 +
5*m + 1)*c^(4*b)*x^(4*b*n)) - 1/2*b^2*m*n^2*x*x^m*e^(-4*a)/((64*b^4*m*n^4 + 64*b^4*n^4 - 20*b^2*m^3*n^2 - 60*b
^2*m^2*n^2 + m^5 - 60*b^2*m*n^2 + 5*m^4 - 20*b^2*n^2 + 10*m^3 + 10*m^2 + 5*m + 1)*c^(4*b)*x^(4*b*n)) + 3/2*m^3
*x*x^m/(64*b^4*m*n^4 + 64*b^4*n^4 - 20*b^2*m^3*n^2 - 60*b^2*m^2*n^2 + m^5 - 60*b^2*m*n^2 + 5*m^4 - 20*b^2*n^2
+ 10*m^3 + 10*m^2 + 5*m + 1) + 1/16*c^(4*b)*x*x^(4*b*n)*x^m*e^(4*a)/(64*b^4*m*n^4 + 64*b^4*n^4 - 20*b^2*m^3*n^
2 - 60*b^2*m^2*n^2 + m^5 - 60*b^2*m*n^2 + 5*m^4 - 20*b^2*n^2 + 10*m^3 + 10*m^2 + 5*m + 1) + 1/4*c^(2*b)*x*x^(2
*b*n)*x^m*e^(2*a)/(64*b^4*m*n^4 + 64*b^4*n^4 - 20*b^2*m^3*n^2 - 60*b^2*m^2*n^2 + m^5 - 60*b^2*m*n^2 + 5*m^4 -
20*b^2*n^2 + 10*m^3 + 10*m^2 + 5*m + 1) + 1/4*m^4*x*x^m*e^(-2*a)/((64*b^4*m*n^4 + 64*b^4*n^4 - 20*b^2*m^3*n^2
- 60*b^2*m^2*n^2 + m^5 - 60*b^2*m*n^2 + 5*m^4 - 20*b^2*n^2 + 10*m^3 + 10*m^2 + 5*m + 1)*c^(2*b)*x^(2*b*n)) + 3
/2*b*m^2*n*x*x^m*e^(-2*a)/((64*b^4*m*n^4 + 64*b^4*n^4 - 20*b^2*m^3*n^2 - 60*b^2*m^2*n^2 + m^5 - 60*b^2*m*n^2 +
 5*m^4 - 20*b^2*n^2 + 10*m^3 + 10*m^2 + 5*m + 1)*c^(2*b)*x^(2*b*n)) - 4*b^2*n^2*x*x^m*e^(-2*a)/((64*b^4*m*n^4
+ 64*b^4*n^4 - 20*b^2*m^3*n^2 - 60*b^2*m^2*n^2 + m^5 - 60*b^2*m*n^2 + 5*m^4 - 20*b^2*n^2 + 10*m^3 + 10*m^2 + 5
*m + 1)*c^(2*b)*x^(2*b*n)) + 1/16*m^4*x*x^m*e^(-4*a)/((64*b^4*m*n^4 + 64*b^4*n^4 - 20*b^2*m^3*n^2 - 60*b^2*m^2
*n^2 + m^5 - 60*b^2*m*n^2 + 5*m^4 - 20*b^2*n^2 + 10*m^3 + 10*m^2 + 5*m + 1)*c^(4*b)*x^(4*b*n)) + 3/4*b*m^2*n*x
*x^m*e^(-4*a)/((64*b^4*m*n^4 + 64*b^4*n^4 - 20*b^2*m^3*n^2 - 60*b^2*m^2*n^2 + m^5 - 60*b^2*m*n^2 + 5*m^4 - 20*
b^2*n^2 + 10*m^3 + 10*m^2 + 5*m + 1)*c^(4*b)*x^(4*b*n)) - 1/4*b^2*n^2*x*x^m*e^(-4*a)/((64*b^4*m*n^4 + 64*b^4*n
^4 - 20*b^2*m^3*n^2 - 60*b^2*m^2*n^2 + m^5 - 60*b^2*m*n^2 + 5*m^4 - 20*b^2*n^2 + 10*m^3 + 10*m^2 + 5*m + 1)*c^
(4*b)*x^(4*b*n)) + 9/4*m^2*x*x^m/(64*b^4*m*n^4 + 64*b^4*n^4 - 20*b^2*m^3*n^2 - 60*b^2*m^2*n^2 + m^5 - 60*b^2*m
*n^2 + 5*m^4 - 20*b^2*n^2 + 10*m^3 + 10*m^2 + 5*m + 1) + m^3*x*x^m*e^(-2*a)/((64*b^4*m*n^4 + 64*b^4*n^4 - 20*b
^2*m^3*n^2 - 60*b^2*m^2*n^2 + m^5 - 60*b^2*m*n^2 + 5*m^4 - 20*b^2*n^2 + 10*m^3 + 10*m^2 + 5*m + 1)*c^(2*b)*x^(
2*b*n)) + 3/2*b*m*n*x*x^m*e^(-2*a)/((64*b^4*m*n^4 + 64*b^4*n^4 - 20*b^2*m^3*n^2 - 60*b^2*m^2*n^2 + m^5 - 60*b^
2*m*n^2 + 5*m^4 - 20*b^2*n^2 + 10*m^3 + 10*m^2 + 5*m + 1)*c^(2*b)*x^(2*b*n)) + 1/4*m^3*x*x^m*e^(-4*a)/((64*b^4
*m*n^4 + 64*b^4*n^4 - 20*b^2*m^3*n^2 - 60*b^2*m^2*n^2 + m^5 - 60*b^2*m*n^2 + 5*m^4 - 20*b^2*n^2 + 10*m^3 + 10*
m^2 + 5*m + 1)*c^(4*b)*x^(4*b*n)) + 3/4*b*m*n*x*x^m*e^(-4*a)/((64*b^4*m*n^4 + 64*b^4*n^4 - 20*b^2*m^3*n^2 - 60
*b^2*m^2*n^2 + m^5 - 60*b^2*m*n^2 + 5*m^4 - 20*b^2*n^2 + 10*m^3 + 10*m^2 + 5*m + 1)*c^(4*b)*x^(4*b*n)) + 3/2*m
*x*x^m/(64*b^4*m*n^4 + 64*b^4*n^4 - 20*b^2*m^3*n^2 - 60*b^2*m^2*n^2 + m^5 - 60*b^2*m*n^2 + 5*m^4 - 20*b^2*n^2
+ 10*m^3 + 10*m^2 + 5*m + 1) + 3/2*m^2*x*x^m*e^(-2*a)/((64*b^4*m*n^4 + 64*b^4*n^4 - 20*b^2*m^3*n^2 - 60*b^2*m^
2*n^2 + m^5 - 60*b^2*m*n^2 + 5*m^4 - 20*b^2*n^2 + 10*m^3 + 10*m^2 + 5*m + 1)*c^(2*b)*x^(2*b*n)) + 1/2*b*n*x*x^
m*e^(-2*a)/((64*b^4*m*n^4 + 64*b^4*n^4 - 20*b^2*m^3*n^2 - 60*b^2*m^2*n^2 + m^5 - 60*b^2*m*n^2 + 5*m^4 - 20*b^2
*n^2 + 10*m^3 + 10*m^2 + 5*m + 1)*c^(2*b)*x^(2*b*n)) + 3/8*m^2*x*x^m*e^(-4*a)/((64*b^4*m*n^4 + 64*b^4*n^4 - 20
*b^2*m^3*n^2 - 60*b^2*m^2*n^2 + m^5 - 60*b^2*m*n^2 + 5*m^4 - 20*b^2*n^2 + 10*m^3 + 10*m^2 + 5*m + 1)*c^(4*b)*x
^(4*b*n)) + 1/4*b*n*x*x^m*e^(-4*a)/((64*b^4*m*n^4 + 64*b^4*n^4 - 20*b^2*m^3*n^2 - 60*b^2*m^2*n^2 + m^5 - 60*b^
2*m*n^2 + 5*m^4 - 20*b^2*n^2 + 10*m^3 + 10*m^2 + 5*m + 1)*c^(4*b)*x^(4*b*n)) + 3/8*x*x^m/(64*b^4*m*n^4 + 64*b^
4*n^4 - 20*b^2*m^3*n^2 - 60*b^2*m^2*n^2 + m^5 - 60*b^2*m*n^2 + 5*m^4 - 20*b^2*n^2 + 10*m^3 + 10*m^2 + 5*m + 1)
 + m*x*x^m*e^(-2*a)/((64*b^4*m*n^4 + 64*b^4*n^4 - 20*b^2*m^3*n^2 - 60*b^2*m^2*n^2 + m^5 - 60*b^2*m*n^2 + 5*m^4
 - 20*b^2*n^2 + 10*m^3 + 10*m^2 + 5*m + 1)*c^(2*b)*x^(2*b*n)) + 1/4*m*x*x^m*e^(-4*a)/((64*b^4*m*n^4 + 64*b^4*n
^4 - 20*b^2*m^3*n^2 - 60*b^2*m^2*n^2 + m^5 - 60*b^2*m*n^2 + 5*m^4 - 20*b^2*n^2 + 10*m^3 + 10*m^2 + 5*m + 1)*c^
(4*b)*x^(4*b*n)) + 1/4*x*x^m*e^(-2*a)/((64*b^4*m*n^4 + 64*b^4*n^4 - 20*b^2*m^3*n^2 - 60*b^2*m^2*n^2 + m^5 - 60
*b^2*m*n^2 + 5*m^4 - 20*b^2*n^2 + 10*m^3 + 10*m^2 + 5*m + 1)*c^(2*b)*x^(2*b*n)) + 1/16*x*x^m*e^(-4*a)/((64*b^4
*m*n^4 + 64*b^4*n^4 - 20*b^2*m^3*n^2 - 60*b^2*m^2*n^2 + m^5 - 60*b^2*m*n^2 + 5*m^4 - 20*b^2*n^2 + 10*m^3 + 10*
m^2 + 5*m + 1)*c^(4*b)*x^(4*b*n))

Mupad [B] (verification not implemented)

Time = 1.97 (sec) , antiderivative size = 134, normalized size of antiderivative = 0.50 \[ \int x^m \cosh ^4\left (a+b \log \left (c x^n\right )\right ) \, dx=\frac {3\,x\,x^m}{8\,m+8}+\frac {x\,x^m\,{\mathrm {e}}^{-2\,a}}{{\left (c\,x^n\right )}^{2\,b}\,\left (4\,m-8\,b\,n+4\right )}+\frac {x\,x^m\,{\mathrm {e}}^{2\,a}\,{\left (c\,x^n\right )}^{2\,b}}{4\,m+8\,b\,n+4}+\frac {x\,x^m\,{\mathrm {e}}^{-4\,a}}{{\left (c\,x^n\right )}^{4\,b}\,\left (16\,m-64\,b\,n+16\right )}+\frac {x\,x^m\,{\mathrm {e}}^{4\,a}\,{\left (c\,x^n\right )}^{4\,b}}{16\,m+64\,b\,n+16} \]

[In]

int(x^m*cosh(a + b*log(c*x^n))^4,x)

[Out]

(3*x*x^m)/(8*m + 8) + (x*x^m*exp(-2*a))/((c*x^n)^(2*b)*(4*m - 8*b*n + 4)) + (x*x^m*exp(2*a)*(c*x^n)^(2*b))/(4*
m + 8*b*n + 4) + (x*x^m*exp(-4*a))/((c*x^n)^(4*b)*(16*m - 64*b*n + 16)) + (x*x^m*exp(4*a)*(c*x^n)^(4*b))/(16*m
 + 64*b*n + 16)