\(\int \frac {x^2 \cosh (c+d x)}{(a+b x)^2} \, dx\) [28]

Optimal result

Integrand size = 17, antiderivative size = 147 \[ \int \frac {x^2 \cosh (c+d x)}{(a+b x)^2} \, dx=-\frac {a^2 \cosh (c+d x)}{b^3 (a+b x)}-\frac {2 a \cosh \left (c-\frac {a d}{b}\right ) \text {Chi}\left (\frac {a d}{b}+d x\right )}{b^3}+\frac {a^2 d \text {Chi}\left (\frac {a d}{b}+d x\right ) \sinh \left (c-\frac {a d}{b}\right )}{b^4}+\frac {\sinh (c+d x)}{b^2 d}+\frac {a^2 d \cosh \left (c-\frac {a d}{b}\right ) \text {Shi}\left (\frac {a d}{b}+d x\right )}{b^4}-\frac {2 a \sinh \left (c-\frac {a d}{b}\right ) \text {Shi}\left (\frac {a d}{b}+d x\right )}{b^3} \] Output:

-a^2*cosh(d*x+c)/b^3/(b*x+a)-2*a*cosh(-c+a*d/b)*Chi(a*d/b+d*x)/b^3-a^2*d*C 
hi(a*d/b+d*x)*sinh(-c+a*d/b)/b^4+sinh(d*x+c)/b^2/d+a^2*d*cosh(-c+a*d/b)*Sh 
i(a*d/b+d*x)/b^4+2*a*sinh(-c+a*d/b)*Shi(a*d/b+d*x)/b^3
 

Mathematica [A] (verified)

Time = 0.43 (sec) , antiderivative size = 115, normalized size of antiderivative = 0.78 \[ \int \frac {x^2 \cosh (c+d x)}{(a+b x)^2} \, dx=\frac {a \text {Chi}\left (d \left (\frac {a}{b}+x\right )\right ) \left (-2 b \cosh \left (c-\frac {a d}{b}\right )+a d \sinh \left (c-\frac {a d}{b}\right )\right )+b \left (-\frac {a^2 \cosh (c+d x)}{a+b x}+\frac {b \sinh (c+d x)}{d}\right )+a \left (a d \cosh \left (c-\frac {a d}{b}\right )-2 b \sinh \left (c-\frac {a d}{b}\right )\right ) \text {Shi}\left (d \left (\frac {a}{b}+x\right )\right )}{b^4} \] Input:

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

Output:

(a*CoshIntegral[d*(a/b + x)]*(-2*b*Cosh[c - (a*d)/b] + a*d*Sinh[c - (a*d)/ 
b]) + b*(-((a^2*Cosh[c + d*x])/(a + b*x)) + (b*Sinh[c + d*x])/d) + a*(a*d* 
Cosh[c - (a*d)/b] - 2*b*Sinh[c - (a*d)/b])*SinhIntegral[d*(a/b + x)])/b^4
 

Rubi [A] (verified)

Time = 0.71 (sec) , antiderivative size = 147, normalized size of antiderivative = 1.00, number of steps used = 2, number of rules used = 2, \(\frac {\text {number of rules}}{\text {integrand size}}\) = 0.118, Rules used = {7293, 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 \frac {x^2 \cosh (c+d x)}{(a+b x)^2} \, dx\)

\(\Big \downarrow \) 7293

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

\(\Big \downarrow \) 2009

\(\displaystyle \frac {a^2 d \sinh \left (c-\frac {a d}{b}\right ) \text {Chi}\left (x d+\frac {a d}{b}\right )}{b^4}+\frac {a^2 d \cosh \left (c-\frac {a d}{b}\right ) \text {Shi}\left (x d+\frac {a d}{b}\right )}{b^4}-\frac {a^2 \cosh (c+d x)}{b^3 (a+b x)}-\frac {2 a \cosh \left (c-\frac {a d}{b}\right ) \text {Chi}\left (x d+\frac {a d}{b}\right )}{b^3}-\frac {2 a \sinh \left (c-\frac {a d}{b}\right ) \text {Shi}\left (x d+\frac {a d}{b}\right )}{b^3}+\frac {\sinh (c+d x)}{b^2 d}\)

Input:

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

Output:

-((a^2*Cosh[c + d*x])/(b^3*(a + b*x))) - (2*a*Cosh[c - (a*d)/b]*CoshIntegr 
al[(a*d)/b + d*x])/b^3 + (a^2*d*CoshIntegral[(a*d)/b + d*x]*Sinh[c - (a*d) 
/b])/b^4 + Sinh[c + d*x]/(b^2*d) + (a^2*d*Cosh[c - (a*d)/b]*SinhIntegral[( 
a*d)/b + d*x])/b^4 - (2*a*Sinh[c - (a*d)/b]*SinhIntegral[(a*d)/b + d*x])/b 
^3
 

Defintions of rubi rules used

rule 2009

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

rule 7293

Int[u_, x_Symbol] :> With[{v = ExpandIntegrand[u, x]}, Int[v, x] /; SumQ[v] 
]
 

Maple [B] (verified)

Leaf count of result is larger than twice the leaf count of optimal. \(435\) vs. \(2(152)=304\).

Time = 0.54 (sec) , antiderivative size = 436, normalized size of antiderivative = 2.97

method	result	size
risch	\(-\frac {{\mathrm e}^{-\frac {a d -c b}{b}} \operatorname {expIntegral}_{1}\left (-d x -c -\frac {a d -c b}{b}\right ) a^{2} b \,d^{2} x -{\mathrm e}^{\frac {a d -c b}{b}} \operatorname {expIntegral}_{1}\left (d x +c +\frac {a d -c b}{b}\right ) a^{2} b \,d^{2} x +{\mathrm e}^{-\frac {a d -c b}{b}} \operatorname {expIntegral}_{1}\left (-d x -c -\frac {a d -c b}{b}\right ) a^{3} d^{2}-2 \,{\mathrm e}^{-\frac {a d -c b}{b}} \operatorname {expIntegral}_{1}\left (-d x -c -\frac {a d -c b}{b}\right ) a \,b^{2} d x -{\mathrm e}^{\frac {a d -c b}{b}} \operatorname {expIntegral}_{1}\left (d x +c +\frac {a d -c b}{b}\right ) a^{3} d^{2}-2 \,{\mathrm e}^{\frac {a d -c b}{b}} \operatorname {expIntegral}_{1}\left (d x +c +\frac {a d -c b}{b}\right ) a \,b^{2} d x -2 \,{\mathrm e}^{-\frac {a d -c b}{b}} \operatorname {expIntegral}_{1}\left (-d x -c -\frac {a d -c b}{b}\right ) a^{2} b d -2 \,{\mathrm e}^{\frac {a d -c b}{b}} \operatorname {expIntegral}_{1}\left (d x +c +\frac {a d -c b}{b}\right ) a^{2} b d +{\mathrm e}^{-d x -c} a^{2} b d +{\mathrm e}^{-d x -c} b^{3} x +{\mathrm e}^{d x +c} a^{2} b d -{\mathrm e}^{d x +c} b^{3} x +{\mathrm e}^{-d x -c} a \,b^{2}-{\mathrm e}^{d x +c} a \,b^{2}}{2 d \,b^{4} \left (b x +a \right )}\)	\(436\)

Input:

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

Output:

-1/2/d*(exp(-(a*d-b*c)/b)*Ei(1,-d*x-c-(a*d-b*c)/b)*a^2*b*d^2*x-exp((a*d-b* 
c)/b)*Ei(1,d*x+c+(a*d-b*c)/b)*a^2*b*d^2*x+exp(-(a*d-b*c)/b)*Ei(1,-d*x-c-(a 
*d-b*c)/b)*a^3*d^2-2*exp(-(a*d-b*c)/b)*Ei(1,-d*x-c-(a*d-b*c)/b)*a*b^2*d*x- 
exp((a*d-b*c)/b)*Ei(1,d*x+c+(a*d-b*c)/b)*a^3*d^2-2*exp((a*d-b*c)/b)*Ei(1,d 
*x+c+(a*d-b*c)/b)*a*b^2*d*x-2*exp(-(a*d-b*c)/b)*Ei(1,-d*x-c-(a*d-b*c)/b)*a 
^2*b*d-2*exp((a*d-b*c)/b)*Ei(1,d*x+c+(a*d-b*c)/b)*a^2*b*d+exp(-d*x-c)*a^2* 
b*d+exp(-d*x-c)*b^3*x+exp(d*x+c)*a^2*b*d-exp(d*x+c)*b^3*x+exp(-d*x-c)*a*b^ 
2-exp(d*x+c)*a*b^2)/b^4/(b*x+a)
 

Fricas [A] (verification not implemented)

Time = 0.08 (sec) , antiderivative size = 274, normalized size of antiderivative = 1.86 \[ \int \frac {x^2 \cosh (c+d x)}{(a+b x)^2} \, dx=-\frac {2 \, a^{2} b d \cosh \left (d x + c\right ) - {\left ({\left (a^{3} d^{2} - 2 \, a^{2} b d + {\left (a^{2} b d^{2} - 2 \, a b^{2} d\right )} x\right )} {\rm Ei}\left (\frac {b d x + a d}{b}\right ) - {\left (a^{3} d^{2} + 2 \, a^{2} b d + {\left (a^{2} b d^{2} + 2 \, a b^{2} d\right )} x\right )} {\rm Ei}\left (-\frac {b d x + a d}{b}\right )\right )} \cosh \left (-\frac {b c - a d}{b}\right ) - 2 \, {\left (b^{3} x + a b^{2}\right )} \sinh \left (d x + c\right ) + {\left ({\left (a^{3} d^{2} - 2 \, a^{2} b d + {\left (a^{2} b d^{2} - 2 \, a b^{2} d\right )} x\right )} {\rm Ei}\left (\frac {b d x + a d}{b}\right ) + {\left (a^{3} d^{2} + 2 \, a^{2} b d + {\left (a^{2} b d^{2} + 2 \, a b^{2} d\right )} x\right )} {\rm Ei}\left (-\frac {b d x + a d}{b}\right )\right )} \sinh \left (-\frac {b c - a d}{b}\right )}{2 \, {\left (b^{5} d x + a b^{4} d\right )}} \] Input:

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

Output:

-1/2*(2*a^2*b*d*cosh(d*x + c) - ((a^3*d^2 - 2*a^2*b*d + (a^2*b*d^2 - 2*a*b 
^2*d)*x)*Ei((b*d*x + a*d)/b) - (a^3*d^2 + 2*a^2*b*d + (a^2*b*d^2 + 2*a*b^2 
*d)*x)*Ei(-(b*d*x + a*d)/b))*cosh(-(b*c - a*d)/b) - 2*(b^3*x + a*b^2)*sinh 
(d*x + c) + ((a^3*d^2 - 2*a^2*b*d + (a^2*b*d^2 - 2*a*b^2*d)*x)*Ei((b*d*x + 
 a*d)/b) + (a^3*d^2 + 2*a^2*b*d + (a^2*b*d^2 + 2*a*b^2*d)*x)*Ei(-(b*d*x + 
a*d)/b))*sinh(-(b*c - a*d)/b))/(b^5*d*x + a*b^4*d)
 

Sympy [F]

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

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

Output:

Integral(x**2*cosh(c + d*x)/(a + b*x)**2, x)
 

Maxima [A] (verification not implemented)

Time = 0.13 (sec) , antiderivative size = 236, normalized size of antiderivative = 1.61 \[ \int \frac {x^2 \cosh (c+d x)}{(a+b x)^2} \, dx=\frac {1}{2} \, {\left (a^{2} {\left (\frac {e^{\left (-c + \frac {a d}{b}\right )} E_{1}\left (\frac {{\left (b x + a\right )} d}{b}\right )}{b^{4}} - \frac {e^{\left (c - \frac {a d}{b}\right )} E_{1}\left (-\frac {{\left (b x + a\right )} d}{b}\right )}{b^{4}}\right )} + \frac {2 \, a {\left (\frac {e^{\left (-c + \frac {a d}{b}\right )} E_{1}\left (\frac {{\left (b x + a\right )} d}{b}\right )}{b} + \frac {e^{\left (c - \frac {a d}{b}\right )} E_{1}\left (-\frac {{\left (b x + a\right )} d}{b}\right )}{b}\right )}}{b^{2} d} - \frac {\frac {{\left (d x e^{c} - e^{c}\right )} e^{\left (d x\right )}}{d^{2}} + \frac {{\left (d x + 1\right )} e^{\left (-d x - c\right )}}{d^{2}}}{b^{2}} + \frac {4 \, a \cosh \left (d x + c\right ) \log \left (b x + a\right )}{b^{3} d}\right )} d - {\left (\frac {a^{2}}{b^{4} x + a b^{3}} - \frac {x}{b^{2}} + \frac {2 \, a \log \left (b x + a\right )}{b^{3}}\right )} \cosh \left (d x + c\right ) \] Input:

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

Output:

1/2*(a^2*(e^(-c + a*d/b)*exp_integral_e(1, (b*x + a)*d/b)/b^4 - e^(c - a*d 
/b)*exp_integral_e(1, -(b*x + a)*d/b)/b^4) + 2*a*(e^(-c + a*d/b)*exp_integ 
ral_e(1, (b*x + a)*d/b)/b + e^(c - a*d/b)*exp_integral_e(1, -(b*x + a)*d/b 
)/b)/(b^2*d) - ((d*x*e^c - e^c)*e^(d*x)/d^2 + (d*x + 1)*e^(-d*x - c)/d^2)/ 
b^2 + 4*a*cosh(d*x + c)*log(b*x + a)/(b^3*d))*d - (a^2/(b^4*x + a*b^3) - x 
/b^2 + 2*a*log(b*x + a)/b^3)*cosh(d*x + c)
 

Giac [B] (verification not implemented)

Leaf count of result is larger than twice the leaf count of optimal. 1308 vs. \(2 (152) = 304\).

Time = 0.15 (sec) , antiderivative size = 1308, normalized size of antiderivative = 8.90 \[ \int \frac {x^2 \cosh (c+d x)}{(a+b x)^2} \, dx=\text {Too large to display} \] Input:

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

Output:

1/2*((b*x + a)*a^2*(b*c/(b*x + a) - a*d/(b*x + a) + d)*d^2*Ei(((b*x + a)*( 
b*c/(b*x + a) - a*d/(b*x + a) + d) - b*c + a*d)/b)*e^((b*c - a*d)/b) - a^2 
*b*c*d^2*Ei(((b*x + a)*(b*c/(b*x + a) - a*d/(b*x + a) + d) - b*c + a*d)/b) 
*e^((b*c - a*d)/b) + a^3*d^3*Ei(((b*x + a)*(b*c/(b*x + a) - a*d/(b*x + a) 
+ d) - b*c + a*d)/b)*e^((b*c - a*d)/b) - (b*x + a)*a^2*(b*c/(b*x + a) - a* 
d/(b*x + a) + d)*d^2*Ei(-((b*x + a)*(b*c/(b*x + a) - a*d/(b*x + a) + d) - 
b*c + a*d)/b)*e^(-(b*c - a*d)/b) + a^2*b*c*d^2*Ei(-((b*x + a)*(b*c/(b*x + 
a) - a*d/(b*x + a) + d) - b*c + a*d)/b)*e^(-(b*c - a*d)/b) - a^3*d^3*Ei(-( 
(b*x + a)*(b*c/(b*x + a) - a*d/(b*x + a) + d) - b*c + a*d)/b)*e^(-(b*c - a 
*d)/b) - 2*(b*x + a)*a*b*(b*c/(b*x + a) - a*d/(b*x + a) + d)*d*Ei(((b*x + 
a)*(b*c/(b*x + a) - a*d/(b*x + a) + d) - b*c + a*d)/b)*e^((b*c - a*d)/b) + 
 2*a*b^2*c*d*Ei(((b*x + a)*(b*c/(b*x + a) - a*d/(b*x + a) + d) - b*c + a*d 
)/b)*e^((b*c - a*d)/b) - 2*a^2*b*d^2*Ei(((b*x + a)*(b*c/(b*x + a) - a*d/(b 
*x + a) + d) - b*c + a*d)/b)*e^((b*c - a*d)/b) - 2*(b*x + a)*a*b*(b*c/(b*x 
 + a) - a*d/(b*x + a) + d)*d*Ei(-((b*x + a)*(b*c/(b*x + a) - a*d/(b*x + a) 
 + d) - b*c + a*d)/b)*e^(-(b*c - a*d)/b) + 2*a*b^2*c*d*Ei(-((b*x + a)*(b*c 
/(b*x + a) - a*d/(b*x + a) + d) - b*c + a*d)/b)*e^(-(b*c - a*d)/b) - 2*a^2 
*b*d^2*Ei(-((b*x + a)*(b*c/(b*x + a) - a*d/(b*x + a) + d) - b*c + a*d)/b)* 
e^(-(b*c - a*d)/b) - a^2*b*d^2*e^((b*x + a)*(b*c/(b*x + a) - a*d/(b*x + a) 
 + d)/b) - a^2*b*d^2*e^(-(b*x + a)*(b*c/(b*x + a) - a*d/(b*x + a) + d)/...
 

Mupad [F(-1)]

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

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

Output:

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

Reduce [F]

\[ \int \frac {x^2 \cosh (c+d x)}{(a+b x)^2} \, dx=\text {too large to display} \] Input:

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

Output:

(e**(2*c + 2*d*x)*a**2*d**2*x - e**(2*c + 2*d*x)*a**2*d - e**(2*c + 2*d*x) 
*a*b - e**(2*c + 2*d*x)*b**2*x - e**(2*c + d*x)*int((e**(d*x)*x)/(a**4*d** 
2 + 2*a**3*b*d**2*x + a**2*b**2*d**2*x**2 - a**2*b**2 - 2*a*b**3*x - b**4* 
x**2),x)*a**6*d**5 - e**(2*c + d*x)*int((e**(d*x)*x)/(a**4*d**2 + 2*a**3*b 
*d**2*x + a**2*b**2*d**2*x**2 - a**2*b**2 - 2*a*b**3*x - b**4*x**2),x)*a** 
5*b*d**5*x + e**(2*c + d*x)*int((e**(d*x)*x)/(a**4*d**2 + 2*a**3*b*d**2*x 
+ a**2*b**2*d**2*x**2 - a**2*b**2 - 2*a*b**3*x - b**4*x**2),x)*a**5*b*d**4 
 + e**(2*c + d*x)*int((e**(d*x)*x)/(a**4*d**2 + 2*a**3*b*d**2*x + a**2*b** 
2*d**2*x**2 - a**2*b**2 - 2*a*b**3*x - b**4*x**2),x)*a**4*b**2*d**4*x + 3* 
e**(2*c + d*x)*int((e**(d*x)*x)/(a**4*d**2 + 2*a**3*b*d**2*x + a**2*b**2*d 
**2*x**2 - a**2*b**2 - 2*a*b**3*x - b**4*x**2),x)*a**4*b**2*d**3 + 3*e**(2 
*c + d*x)*int((e**(d*x)*x)/(a**4*d**2 + 2*a**3*b*d**2*x + a**2*b**2*d**2*x 
**2 - a**2*b**2 - 2*a*b**3*x - b**4*x**2),x)*a**3*b**3*d**3*x - e**(2*c + 
d*x)*int((e**(d*x)*x)/(a**4*d**2 + 2*a**3*b*d**2*x + a**2*b**2*d**2*x**2 - 
 a**2*b**2 - 2*a*b**3*x - b**4*x**2),x)*a**3*b**3*d**2 - e**(2*c + d*x)*in 
t((e**(d*x)*x)/(a**4*d**2 + 2*a**3*b*d**2*x + a**2*b**2*d**2*x**2 - a**2*b 
**2 - 2*a*b**3*x - b**4*x**2),x)*a**2*b**4*d**2*x - 2*e**(2*c + d*x)*int(( 
e**(d*x)*x)/(a**4*d**2 + 2*a**3*b*d**2*x + a**2*b**2*d**2*x**2 - a**2*b**2 
 - 2*a*b**3*x - b**4*x**2),x)*a**2*b**4*d - 2*e**(2*c + d*x)*int((e**(d*x) 
*x)/(a**4*d**2 + 2*a**3*b*d**2*x + a**2*b**2*d**2*x**2 - a**2*b**2 - 2*...