\(\int x^{3+m} \cosh ^2(a+b x) \, dx\) [88]

Optimal result

Integrand size = 14, antiderivative size = 86 \[ \int x^{3+m} \cosh ^2(a+b x) \, dx=\frac {x^{4+m}}{2 (4+m)}-\frac {2^{-6-m} e^{2 a} x^m (-b x)^{-m} \Gamma (4+m,-2 b x)}{b^4}-\frac {2^{-6-m} e^{-2 a} x^m (b x)^{-m} \Gamma (4+m,2 b x)}{b^4} \] Output:

x^(4+m)/(8+2*m)-2^(-6-m)*exp(2*a)*x^m*GAMMA(4+m,-2*b*x)/b^4/((-b*x)^m)-2^( 
-6-m)*x^m*GAMMA(4+m,2*b*x)/b^4/exp(2*a)/((b*x)^m)
 

Mathematica [A] (verified)

Time = 0.07 (sec) , antiderivative size = 79, normalized size of antiderivative = 0.92 \[ \int x^{3+m} \cosh ^2(a+b x) \, dx=\frac {1}{64} x^m \left (\frac {32 x^4}{4+m}-\frac {2^{-m} e^{2 a} (-b x)^{-m} \Gamma (4+m,-2 b x)}{b^4}-\frac {2^{-m} e^{-2 a} (b x)^{-m} \Gamma (4+m,2 b x)}{b^4}\right ) \] Input:

Integrate[x^(3 + m)*Cosh[a + b*x]^2,x]
 

Output:

(x^m*((32*x^4)/(4 + m) - (E^(2*a)*Gamma[4 + m, -2*b*x])/(2^m*b^4*(-(b*x))^ 
m) - Gamma[4 + m, 2*b*x]/(2^m*b^4*E^(2*a)*(b*x)^m)))/64
 

Rubi [A] (verified)

Time = 0.37 (sec) , antiderivative size = 86, normalized size of antiderivative = 1.00, number of steps used = 3, number of rules used = 3, \(\frac {\text {number of rules}}{\text {integrand size}}\) = 0.214, Rules used = {3042, 3793, 2009}

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 x^{m+3} \cosh ^2(a+b x) \, dx\)

\(\Big \downarrow \) 3042

\(\displaystyle \int x^{m+3} \sin \left (i a+i b x+\frac {\pi }{2}\right )^2dx\)

\(\Big \downarrow \) 3793

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

\(\Big \downarrow \) 2009

\(\displaystyle -\frac {e^{2 a} 2^{-m-6} x^m (-b x)^{-m} \Gamma (m+4,-2 b x)}{b^4}-\frac {e^{-2 a} 2^{-m-6} x^m (b x)^{-m} \Gamma (m+4,2 b x)}{b^4}+\frac {x^{m+4}}{2 (m+4)}\)

Input:

Int[x^(3 + m)*Cosh[a + b*x]^2,x]
 

Output:

x^(4 + m)/(2*(4 + m)) - (2^(-6 - m)*E^(2*a)*x^m*Gamma[4 + m, -2*b*x])/(b^4 
*(-(b*x))^m) - (2^(-6 - m)*x^m*Gamma[4 + m, 2*b*x])/(b^4*E^(2*a)*(b*x)^m)
 

Defintions of rubi rules used

rule 2009

Int[u_, x_Symbol] :> Simp[IntSum[u, x], x] /; SumQ[u]
 

rule 3042

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

rule 3793

Int[((c_.) + (d_.)*(x_))^(m_)*sin[(e_.) + (f_.)*(x_)]^(n_), x_Symbol] :> In 
t[ExpandTrigReduce[(c + d*x)^m, Sin[e + f*x]^n, x], x] /; FreeQ[{c, d, e, f 
, m}, x] && IGtQ[n, 1] && ( !RationalQ[m] || (GeQ[m, -1] && LtQ[m, 1]))
 

Maple [F]

\[\int x^{3+m} \cosh \left (b x +a \right )^{2}d x\]

Input:

int(x^(3+m)*cosh(b*x+a)^2,x)
 

Output:

int(x^(3+m)*cosh(b*x+a)^2,x)
 

Fricas [A] (verification not implemented)

Time = 0.10 (sec) , antiderivative size = 136, normalized size of antiderivative = 1.58 \[ \int x^{3+m} \cosh ^2(a+b x) \, dx=\frac {4 \, b x \cosh \left ({\left (m + 3\right )} \log \left (x\right )\right ) - {\left (m + 4\right )} \cosh \left ({\left (m + 3\right )} \log \left (2 \, b\right ) + 2 \, a\right ) \Gamma \left (m + 4, 2 \, b x\right ) + {\left (m + 4\right )} \cosh \left ({\left (m + 3\right )} \log \left (-2 \, b\right ) - 2 \, a\right ) \Gamma \left (m + 4, -2 \, b x\right ) + {\left (m + 4\right )} \Gamma \left (m + 4, 2 \, b x\right ) \sinh \left ({\left (m + 3\right )} \log \left (2 \, b\right ) + 2 \, a\right ) - {\left (m + 4\right )} \Gamma \left (m + 4, -2 \, b x\right ) \sinh \left ({\left (m + 3\right )} \log \left (-2 \, b\right ) - 2 \, a\right ) + 4 \, b x \sinh \left ({\left (m + 3\right )} \log \left (x\right )\right )}{8 \, {\left (b m + 4 \, b\right )}} \] Input:

integrate(x^(3+m)*cosh(b*x+a)^2,x, algorithm="fricas")
 

Output:

1/8*(4*b*x*cosh((m + 3)*log(x)) - (m + 4)*cosh((m + 3)*log(2*b) + 2*a)*gam 
ma(m + 4, 2*b*x) + (m + 4)*cosh((m + 3)*log(-2*b) - 2*a)*gamma(m + 4, -2*b 
*x) + (m + 4)*gamma(m + 4, 2*b*x)*sinh((m + 3)*log(2*b) + 2*a) - (m + 4)*g 
amma(m + 4, -2*b*x)*sinh((m + 3)*log(-2*b) - 2*a) + 4*b*x*sinh((m + 3)*log 
(x)))/(b*m + 4*b)
 

Sympy [F]

\[ \int x^{3+m} \cosh ^2(a+b x) \, dx=\int x^{m + 3} \cosh ^{2}{\left (a + b x \right )}\, dx \] Input:

integrate(x**(3+m)*cosh(b*x+a)**2,x)
 

Output:

Integral(x**(m + 3)*cosh(a + b*x)**2, x)
 

Maxima [A] (verification not implemented)

Time = 0.11 (sec) , antiderivative size = 71, normalized size of antiderivative = 0.83 \[ \int x^{3+m} \cosh ^2(a+b x) \, dx=-\frac {1}{4} \, \left (2 \, b x\right )^{-m - 4} x^{m + 4} e^{\left (-2 \, a\right )} \Gamma \left (m + 4, 2 \, b x\right ) - \frac {1}{4} \, \left (-2 \, b x\right )^{-m - 4} x^{m + 4} e^{\left (2 \, a\right )} \Gamma \left (m + 4, -2 \, b x\right ) + \frac {x^{m + 4}}{2 \, {\left (m + 4\right )}} \] Input:

integrate(x^(3+m)*cosh(b*x+a)^2,x, algorithm="maxima")
 

Output:

-1/4*(2*b*x)^(-m - 4)*x^(m + 4)*e^(-2*a)*gamma(m + 4, 2*b*x) - 1/4*(-2*b*x 
)^(-m - 4)*x^(m + 4)*e^(2*a)*gamma(m + 4, -2*b*x) + 1/2*x^(m + 4)/(m + 4)
 

Giac [F]

\[ \int x^{3+m} \cosh ^2(a+b x) \, dx=\int { x^{m + 3} \cosh \left (b x + a\right )^{2} \,d x } \] Input:

integrate(x^(3+m)*cosh(b*x+a)^2,x, algorithm="giac")
 

Output:

integrate(x^(m + 3)*cosh(b*x + a)^2, x)
 

Mupad [F(-1)]

Timed out. \[ \int x^{3+m} \cosh ^2(a+b x) \, dx=\int x^{m+3}\,{\mathrm {cosh}\left (a+b\,x\right )}^2 \,d x \] Input:

int(x^(m + 3)*cosh(a + b*x)^2,x)
 

Output:

int(x^(m + 3)*cosh(a + b*x)^2, x)
 

Reduce [F]

\[ \int x^{3+m} \cosh ^2(a+b x) \, dx=\frac {10 e^{2 b x +4 a} \left (\int \frac {x^{m} e^{2 b x}}{x}d x \right ) m^{4}+35 e^{2 b x +4 a} \left (\int \frac {x^{m} e^{2 b x}}{x}d x \right ) m^{3}+50 e^{2 b x +4 a} \left (\int \frac {x^{m} e^{2 b x}}{x}d x \right ) m^{2}+24 e^{2 b x +4 a} \left (\int \frac {x^{m} e^{2 b x}}{x}d x \right ) m -x^{m} m^{4}-24 x^{m} e^{4 b x +4 a}-10 x^{m} e^{4 b x +4 a} m^{3}-35 x^{m} e^{4 b x +4 a} m^{2}-50 x^{m} e^{4 b x +4 a} m -32 x^{m} b^{3} x^{3}-48 x^{m} b x +32 x^{m} e^{4 b x +4 a} b^{3} x^{3}-48 x^{m} e^{4 b x +4 a} b^{2} x^{2}+48 x^{m} e^{4 b x +4 a} b x +32 x^{m} e^{2 b x +2 a} b^{4} x^{4}-8 x^{m} b^{3} m \,x^{3}-4 x^{m} b^{2} m^{2} x^{2}-28 x^{m} b^{2} m \,x^{2}-2 x^{m} b \,m^{3} x -18 x^{m} b \,m^{2} x +10 e^{2 b x} \left (\int \frac {x^{m}}{e^{2 b x} x}d x \right ) m^{4}+35 e^{2 b x} \left (\int \frac {x^{m}}{e^{2 b x} x}d x \right ) m^{3}+50 e^{2 b x} \left (\int \frac {x^{m}}{e^{2 b x} x}d x \right ) m^{2}+24 e^{2 b x} \left (\int \frac {x^{m}}{e^{2 b x} x}d x \right ) m +e^{2 b x} \left (\int \frac {x^{m}}{e^{2 b x} x}d x \right ) m^{5}+e^{2 b x +4 a} \left (\int \frac {x^{m} e^{2 b x}}{x}d x \right ) m^{5}+8 x^{m} e^{4 b x +4 a} b^{3} m \,x^{3}-4 x^{m} e^{4 b x +4 a} b^{2} m^{2} x^{2}-28 x^{m} e^{4 b x +4 a} b^{2} m \,x^{2}+2 x^{m} e^{4 b x +4 a} b \,m^{3} x +18 x^{m} e^{4 b x +4 a} b \,m^{2} x +52 x^{m} e^{4 b x +4 a} b m x -10 x^{m} m^{3}-24 x^{m}-x^{m} e^{4 b x +4 a} m^{4}-35 x^{m} m^{2}-50 x^{m} m -48 x^{m} b^{2} x^{2}-52 x^{m} b m x}{64 e^{2 b x +2 a} b^{4} \left (m +4\right )} \] Input:

int(x^(3+m)*cosh(b*x+a)^2,x)
 

Output:

(8*x**m*e**(4*a + 4*b*x)*b**3*m*x**3 + 32*x**m*e**(4*a + 4*b*x)*b**3*x**3 
- 4*x**m*e**(4*a + 4*b*x)*b**2*m**2*x**2 - 28*x**m*e**(4*a + 4*b*x)*b**2*m 
*x**2 - 48*x**m*e**(4*a + 4*b*x)*b**2*x**2 + 2*x**m*e**(4*a + 4*b*x)*b*m** 
3*x + 18*x**m*e**(4*a + 4*b*x)*b*m**2*x + 52*x**m*e**(4*a + 4*b*x)*b*m*x + 
 48*x**m*e**(4*a + 4*b*x)*b*x - x**m*e**(4*a + 4*b*x)*m**4 - 10*x**m*e**(4 
*a + 4*b*x)*m**3 - 35*x**m*e**(4*a + 4*b*x)*m**2 - 50*x**m*e**(4*a + 4*b*x 
)*m - 24*x**m*e**(4*a + 4*b*x) + e**(4*a + 2*b*x)*int((x**m*e**(2*b*x))/x, 
x)*m**5 + 10*e**(4*a + 2*b*x)*int((x**m*e**(2*b*x))/x,x)*m**4 + 35*e**(4*a 
 + 2*b*x)*int((x**m*e**(2*b*x))/x,x)*m**3 + 50*e**(4*a + 2*b*x)*int((x**m* 
e**(2*b*x))/x,x)*m**2 + 24*e**(4*a + 2*b*x)*int((x**m*e**(2*b*x))/x,x)*m + 
 32*x**m*e**(2*a + 2*b*x)*b**4*x**4 + e**(2*b*x)*int(x**m/(e**(2*b*x)*x),x 
)*m**5 + 10*e**(2*b*x)*int(x**m/(e**(2*b*x)*x),x)*m**4 + 35*e**(2*b*x)*int 
(x**m/(e**(2*b*x)*x),x)*m**3 + 50*e**(2*b*x)*int(x**m/(e**(2*b*x)*x),x)*m* 
*2 + 24*e**(2*b*x)*int(x**m/(e**(2*b*x)*x),x)*m - 8*x**m*b**3*m*x**3 - 32* 
x**m*b**3*x**3 - 4*x**m*b**2*m**2*x**2 - 28*x**m*b**2*m*x**2 - 48*x**m*b** 
2*x**2 - 2*x**m*b*m**3*x - 18*x**m*b*m**2*x - 52*x**m*b*m*x - 48*x**m*b*x 
- x**m*m**4 - 10*x**m*m**3 - 35*x**m*m**2 - 50*x**m*m - 24*x**m)/(64*e**(2 
*a + 2*b*x)*b**4*(m + 4))