3.4.46 \(\int f^{a+b x} \sinh ^2(d+e x+f x^2) \, dx\) [346]

3.4.46.1 Optimal result
3.4.46.2 Mathematica [A] (verified)
3.4.46.3 Rubi [A] (verified)
3.4.46.4 Maple [A] (verified)
3.4.46.5 Fricas [B] (verification not implemented)
3.4.46.6 Sympy [F]
3.4.46.7 Maxima [A] (verification not implemented)
3.4.46.8 Giac [C] (verification not implemented)
3.4.46.9 Mupad [F(-1)]

3.4.46.1 Optimal result

Integrand size = 21, antiderivative size = 161 \[ \int f^{a+b x} \sinh ^2\left (d+e x+f x^2\right ) \, dx=\frac {1}{8} e^{-2 d+\frac {(2 e-b \log (f))^2}{8 f}} f^{-\frac {1}{2}+a} \sqrt {\frac {\pi }{2}} \text {erf}\left (\frac {2 e+4 f x-b \log (f)}{2 \sqrt {2} \sqrt {f}}\right )+\frac {1}{8} e^{2 d-\frac {(2 e+b \log (f))^2}{8 f}} f^{-\frac {1}{2}+a} \sqrt {\frac {\pi }{2}} \text {erfi}\left (\frac {2 e+4 f x+b \log (f)}{2 \sqrt {2} \sqrt {f}}\right )-\frac {f^{a+b x}}{2 b \log (f)} \]

output
-1/2*f^(b*x+a)/b/ln(f)+1/16*exp(-2*d+1/8*(2*e-b*ln(f))^2/f)*f^(-1/2+a)*erf 
(1/4*(2*e+4*f*x-b*ln(f))*2^(1/2)/f^(1/2))*2^(1/2)*Pi^(1/2)+1/16*exp(2*d-1/ 
8*(2*e+b*ln(f))^2/f)*f^(-1/2+a)*erfi(1/4*(2*e+4*f*x+b*ln(f))*2^(1/2)/f^(1/ 
2))*2^(1/2)*Pi^(1/2)
 
3.4.46.2 Mathematica [A] (verified)

Time = 0.45 (sec) , antiderivative size = 220, normalized size of antiderivative = 1.37 \[ \int f^{a+b x} \sinh ^2\left (d+e x+f x^2\right ) \, dx=\frac {e^{-\frac {4 e^2+b^2 \log ^2(f)}{8 f}} f^{a-\frac {b e+f}{2 f}} \left (-4 \sqrt {2} e^{\frac {4 e^2+b^2 \log ^2(f)}{8 f}} f^{\frac {1}{2}+b \left (\frac {e}{2 f}+x\right )}+b e^{\frac {4 e^2+b^2 \log ^2(f)}{4 f}} \sqrt {\pi } \text {erf}\left (\frac {2 e+4 f x-b \log (f)}{2 \sqrt {2} \sqrt {f}}\right ) \log (f) (\cosh (2 d)-\sinh (2 d))+b \sqrt {\pi } \text {erfi}\left (\frac {2 e+4 f x+b \log (f)}{2 \sqrt {2} \sqrt {f}}\right ) \log (f) (\cosh (2 d)+\sinh (2 d))\right )}{8 \sqrt {2} b \log (f)} \]

input
Integrate[f^(a + b*x)*Sinh[d + e*x + f*x^2]^2,x]
 
output
(f^(a - (b*e + f)/(2*f))*(-4*Sqrt[2]*E^((4*e^2 + b^2*Log[f]^2)/(8*f))*f^(1 
/2 + b*(e/(2*f) + x)) + b*E^((4*e^2 + b^2*Log[f]^2)/(4*f))*Sqrt[Pi]*Erf[(2 
*e + 4*f*x - b*Log[f])/(2*Sqrt[2]*Sqrt[f])]*Log[f]*(Cosh[2*d] - Sinh[2*d]) 
 + b*Sqrt[Pi]*Erfi[(2*e + 4*f*x + b*Log[f])/(2*Sqrt[2]*Sqrt[f])]*Log[f]*(C 
osh[2*d] + Sinh[2*d])))/(8*Sqrt[2]*b*E^((4*e^2 + b^2*Log[f]^2)/(8*f))*Log[ 
f])
 
3.4.46.3 Rubi [A] (verified)

Time = 0.51 (sec) , antiderivative size = 161, 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.095, Rules used = {6038, 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 f^{a+b x} \sinh ^2\left (d+e x+f x^2\right ) \, dx\)

\(\Big \downarrow \) 6038

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

\(\Big \downarrow \) 2009

\(\displaystyle \frac {1}{8} \sqrt {\frac {\pi }{2}} f^{a-\frac {1}{2}} e^{\frac {(2 e-b \log (f))^2}{8 f}-2 d} \text {erf}\left (\frac {-b \log (f)+2 e+4 f x}{2 \sqrt {2} \sqrt {f}}\right )+\frac {1}{8} \sqrt {\frac {\pi }{2}} f^{a-\frac {1}{2}} e^{2 d-\frac {(b \log (f)+2 e)^2}{8 f}} \text {erfi}\left (\frac {b \log (f)+2 e+4 f x}{2 \sqrt {2} \sqrt {f}}\right )-\frac {f^{a+b x}}{2 b \log (f)}\)

input
Int[f^(a + b*x)*Sinh[d + e*x + f*x^2]^2,x]
 
output
(E^(-2*d + (2*e - b*Log[f])^2/(8*f))*f^(-1/2 + a)*Sqrt[Pi/2]*Erf[(2*e + 4* 
f*x - b*Log[f])/(2*Sqrt[2]*Sqrt[f])])/8 + (E^(2*d - (2*e + b*Log[f])^2/(8* 
f))*f^(-1/2 + a)*Sqrt[Pi/2]*Erfi[(2*e + 4*f*x + b*Log[f])/(2*Sqrt[2]*Sqrt[ 
f])])/8 - f^(a + b*x)/(2*b*Log[f])
 

3.4.46.3.1 Defintions of rubi rules used

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

rule 6038
Int[(F_)^(u_)*Sinh[v_]^(n_.), x_Symbol] :> Int[ExpandTrigToExp[F^u, Sinh[v] 
^n, x], x] /; FreeQ[F, x] && (LinearQ[u, x] || PolyQ[u, x, 2]) && (LinearQ[ 
v, x] || PolyQ[v, x, 2]) && IGtQ[n, 0]
 
3.4.46.4 Maple [A] (verified)

Time = 0.89 (sec) , antiderivative size = 158, normalized size of antiderivative = 0.98

method result size
risch \(-\frac {\operatorname {erf}\left (-\sqrt {2}\, \sqrt {f}\, x +\frac {\left (b \ln \left (f \right )-2 e \right ) \sqrt {2}}{4 \sqrt {f}}\right ) \sqrt {2}\, \sqrt {\pi }\, f^{a} {\mathrm e}^{\frac {b^{2} \ln \left (f \right )^{2}-4 \ln \left (f \right ) b e -16 d f +4 e^{2}}{8 f}}}{16 \sqrt {f}}-\frac {\operatorname {erf}\left (-\sqrt {-2 f}\, x +\frac {2 e +b \ln \left (f \right )}{2 \sqrt {-2 f}}\right ) \sqrt {\pi }\, f^{a} {\mathrm e}^{-\frac {b^{2} \ln \left (f \right )^{2}+4 \ln \left (f \right ) b e -16 d f +4 e^{2}}{8 f}}}{8 \sqrt {-2 f}}-\frac {f^{a} f^{b x}}{2 b \ln \left (f \right )}\) \(158\)

input
int(f^(b*x+a)*sinh(f*x^2+e*x+d)^2,x,method=_RETURNVERBOSE)
 
output
-1/16*erf(-2^(1/2)*f^(1/2)*x+1/4*(b*ln(f)-2*e)*2^(1/2)/f^(1/2))/f^(1/2)*2^ 
(1/2)*Pi^(1/2)*f^a*exp(1/8*(b^2*ln(f)^2-4*ln(f)*b*e-16*d*f+4*e^2)/f)-1/8*e 
rf(-(-2*f)^(1/2)*x+1/2*(2*e+b*ln(f))/(-2*f)^(1/2))/(-2*f)^(1/2)*Pi^(1/2)*f 
^a*exp(-1/8*(b^2*ln(f)^2+4*ln(f)*b*e-16*d*f+4*e^2)/f)-1/2*f^a*f^(b*x)/b/ln 
(f)
 
3.4.46.5 Fricas [B] (verification not implemented)

Leaf count of result is larger than twice the leaf count of optimal. 334 vs. \(2 (126) = 252\).

Time = 0.30 (sec) , antiderivative size = 334, normalized size of antiderivative = 2.07 \[ \int f^{a+b x} \sinh ^2\left (d+e x+f x^2\right ) \, dx=-\frac {\sqrt {2} \sqrt {\pi } b \sqrt {-f} \cosh \left (\frac {b^{2} \log \left (f\right )^{2} + 4 \, e^{2} - 16 \, d f + 4 \, {\left (b e - 2 \, a f\right )} \log \left (f\right )}{8 \, f}\right ) \operatorname {erf}\left (\frac {\sqrt {2} {\left (4 \, f x + b \log \left (f\right ) + 2 \, e\right )} \sqrt {-f}}{4 \, f}\right ) \log \left (f\right ) + \sqrt {2} \sqrt {\pi } b \sqrt {f} \cosh \left (\frac {b^{2} \log \left (f\right )^{2} + 4 \, e^{2} - 16 \, d f - 4 \, {\left (b e - 2 \, a f\right )} \log \left (f\right )}{8 \, f}\right ) \operatorname {erf}\left (-\frac {\sqrt {2} {\left (4 \, f x - b \log \left (f\right ) + 2 \, e\right )}}{4 \, \sqrt {f}}\right ) \log \left (f\right ) - \sqrt {2} \sqrt {\pi } b \sqrt {-f} \operatorname {erf}\left (\frac {\sqrt {2} {\left (4 \, f x + b \log \left (f\right ) + 2 \, e\right )} \sqrt {-f}}{4 \, f}\right ) \log \left (f\right ) \sinh \left (\frac {b^{2} \log \left (f\right )^{2} + 4 \, e^{2} - 16 \, d f + 4 \, {\left (b e - 2 \, a f\right )} \log \left (f\right )}{8 \, f}\right ) + \sqrt {2} \sqrt {\pi } b \sqrt {f} \operatorname {erf}\left (-\frac {\sqrt {2} {\left (4 \, f x - b \log \left (f\right ) + 2 \, e\right )}}{4 \, \sqrt {f}}\right ) \log \left (f\right ) \sinh \left (\frac {b^{2} \log \left (f\right )^{2} + 4 \, e^{2} - 16 \, d f - 4 \, {\left (b e - 2 \, a f\right )} \log \left (f\right )}{8 \, f}\right ) + 8 \, f \cosh \left ({\left (b x + a\right )} \log \left (f\right )\right ) + 8 \, f \sinh \left ({\left (b x + a\right )} \log \left (f\right )\right )}{16 \, b f \log \left (f\right )} \]

input
integrate(f^(b*x+a)*sinh(f*x^2+e*x+d)^2,x, algorithm="fricas")
 
output
-1/16*(sqrt(2)*sqrt(pi)*b*sqrt(-f)*cosh(1/8*(b^2*log(f)^2 + 4*e^2 - 16*d*f 
 + 4*(b*e - 2*a*f)*log(f))/f)*erf(1/4*sqrt(2)*(4*f*x + b*log(f) + 2*e)*sqr 
t(-f)/f)*log(f) + sqrt(2)*sqrt(pi)*b*sqrt(f)*cosh(1/8*(b^2*log(f)^2 + 4*e^ 
2 - 16*d*f - 4*(b*e - 2*a*f)*log(f))/f)*erf(-1/4*sqrt(2)*(4*f*x - b*log(f) 
 + 2*e)/sqrt(f))*log(f) - sqrt(2)*sqrt(pi)*b*sqrt(-f)*erf(1/4*sqrt(2)*(4*f 
*x + b*log(f) + 2*e)*sqrt(-f)/f)*log(f)*sinh(1/8*(b^2*log(f)^2 + 4*e^2 - 1 
6*d*f + 4*(b*e - 2*a*f)*log(f))/f) + sqrt(2)*sqrt(pi)*b*sqrt(f)*erf(-1/4*s 
qrt(2)*(4*f*x - b*log(f) + 2*e)/sqrt(f))*log(f)*sinh(1/8*(b^2*log(f)^2 + 4 
*e^2 - 16*d*f - 4*(b*e - 2*a*f)*log(f))/f) + 8*f*cosh((b*x + a)*log(f)) + 
8*f*sinh((b*x + a)*log(f)))/(b*f*log(f))
 
3.4.46.6 Sympy [F]

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

input
integrate(f**(b*x+a)*sinh(f*x**2+e*x+d)**2,x)
 
output
Integral(f**(a + b*x)*sinh(d + e*x + f*x**2)**2, x)
 
3.4.46.7 Maxima [A] (verification not implemented)

Time = 0.26 (sec) , antiderivative size = 143, normalized size of antiderivative = 0.89 \[ \int f^{a+b x} \sinh ^2\left (d+e x+f x^2\right ) \, dx=\frac {\sqrt {2} \sqrt {\pi } f^{a} \operatorname {erf}\left (\sqrt {2} \sqrt {-f} x - \frac {\sqrt {2} {\left (b \log \left (f\right ) + 2 \, e\right )}}{4 \, \sqrt {-f}}\right ) e^{\left (2 \, d - \frac {{\left (b \log \left (f\right ) + 2 \, e\right )}^{2}}{8 \, f}\right )}}{16 \, \sqrt {-f}} + \frac {\sqrt {2} \sqrt {\pi } f^{a} \operatorname {erf}\left (\sqrt {2} \sqrt {f} x - \frac {\sqrt {2} {\left (b \log \left (f\right ) - 2 \, e\right )}}{4 \, \sqrt {f}}\right ) e^{\left (-2 \, d + \frac {{\left (b \log \left (f\right ) - 2 \, e\right )}^{2}}{8 \, f}\right )}}{16 \, \sqrt {f}} - \frac {f^{b x + a}}{2 \, b \log \left (f\right )} \]

input
integrate(f^(b*x+a)*sinh(f*x^2+e*x+d)^2,x, algorithm="maxima")
 
output
1/16*sqrt(2)*sqrt(pi)*f^a*erf(sqrt(2)*sqrt(-f)*x - 1/4*sqrt(2)*(b*log(f) + 
 2*e)/sqrt(-f))*e^(2*d - 1/8*(b*log(f) + 2*e)^2/f)/sqrt(-f) + 1/16*sqrt(2) 
*sqrt(pi)*f^a*erf(sqrt(2)*sqrt(f)*x - 1/4*sqrt(2)*(b*log(f) - 2*e)/sqrt(f) 
)*e^(-2*d + 1/8*(b*log(f) - 2*e)^2/f)/sqrt(f) - 1/2*f^(b*x + a)/(b*log(f))
 
3.4.46.8 Giac [C] (verification not implemented)

Result contains complex when optimal does not.

Time = 0.28 (sec) , antiderivative size = 388, normalized size of antiderivative = 2.41 \[ \int f^{a+b x} \sinh ^2\left (d+e x+f x^2\right ) \, dx=-\frac {\sqrt {2} \sqrt {\pi } \operatorname {erf}\left (-\frac {1}{4} \, \sqrt {2} \sqrt {-f} {\left (4 \, x + \frac {b \log \left (f\right ) + 2 \, e}{f}\right )}\right ) e^{\left (-\frac {b^{2} \log \left (f\right )^{2} + 4 \, b e \log \left (f\right ) - 8 \, a f \log \left (f\right ) + 4 \, e^{2} - 16 \, d f}{8 \, f}\right )}}{16 \, \sqrt {-f}} - \frac {\sqrt {2} \sqrt {\pi } \operatorname {erf}\left (-\frac {1}{4} \, \sqrt {2} \sqrt {f} {\left (4 \, x - \frac {b \log \left (f\right ) - 2 \, e}{f}\right )}\right ) e^{\left (\frac {b^{2} \log \left (f\right )^{2} - 4 \, b e \log \left (f\right ) + 8 \, a f \log \left (f\right ) + 4 \, e^{2} - 16 \, d f}{8 \, f}\right )}}{16 \, \sqrt {f}} - {\left (\frac {2 \, b \cos \left (-\frac {1}{2} \, \pi b x \mathrm {sgn}\left (f\right ) + \frac {1}{2} \, \pi b x - \frac {1}{2} \, \pi a \mathrm {sgn}\left (f\right ) + \frac {1}{2} \, \pi a\right ) \log \left ({\left | f \right |}\right )}{4 \, b^{2} \log \left ({\left | f \right |}\right )^{2} + {\left (\pi b \mathrm {sgn}\left (f\right ) - \pi b\right )}^{2}} - \frac {{\left (\pi b \mathrm {sgn}\left (f\right ) - \pi b\right )} \sin \left (-\frac {1}{2} \, \pi b x \mathrm {sgn}\left (f\right ) + \frac {1}{2} \, \pi b x - \frac {1}{2} \, \pi a \mathrm {sgn}\left (f\right ) + \frac {1}{2} \, \pi a\right )}{4 \, b^{2} \log \left ({\left | f \right |}\right )^{2} + {\left (\pi b \mathrm {sgn}\left (f\right ) - \pi b\right )}^{2}}\right )} e^{\left (b x \log \left ({\left | f \right |}\right ) + a \log \left ({\left | f \right |}\right )\right )} + i \, {\left (-\frac {i \, e^{\left (\frac {1}{2} i \, \pi b x \mathrm {sgn}\left (f\right ) - \frac {1}{2} i \, \pi b x + \frac {1}{2} i \, \pi a \mathrm {sgn}\left (f\right ) - \frac {1}{2} i \, \pi a\right )}}{2 i \, \pi b \mathrm {sgn}\left (f\right ) - 2 i \, \pi b + 4 \, b \log \left ({\left | f \right |}\right )} + \frac {i \, e^{\left (-\frac {1}{2} i \, \pi b x \mathrm {sgn}\left (f\right ) + \frac {1}{2} i \, \pi b x - \frac {1}{2} i \, \pi a \mathrm {sgn}\left (f\right ) + \frac {1}{2} i \, \pi a\right )}}{-2 i \, \pi b \mathrm {sgn}\left (f\right ) + 2 i \, \pi b + 4 \, b \log \left ({\left | f \right |}\right )}\right )} e^{\left (b x \log \left ({\left | f \right |}\right ) + a \log \left ({\left | f \right |}\right )\right )} \]

input
integrate(f^(b*x+a)*sinh(f*x^2+e*x+d)^2,x, algorithm="giac")
 
output
-1/16*sqrt(2)*sqrt(pi)*erf(-1/4*sqrt(2)*sqrt(-f)*(4*x + (b*log(f) + 2*e)/f 
))*e^(-1/8*(b^2*log(f)^2 + 4*b*e*log(f) - 8*a*f*log(f) + 4*e^2 - 16*d*f)/f 
)/sqrt(-f) - 1/16*sqrt(2)*sqrt(pi)*erf(-1/4*sqrt(2)*sqrt(f)*(4*x - (b*log( 
f) - 2*e)/f))*e^(1/8*(b^2*log(f)^2 - 4*b*e*log(f) + 8*a*f*log(f) + 4*e^2 - 
 16*d*f)/f)/sqrt(f) - (2*b*cos(-1/2*pi*b*x*sgn(f) + 1/2*pi*b*x - 1/2*pi*a* 
sgn(f) + 1/2*pi*a)*log(abs(f))/(4*b^2*log(abs(f))^2 + (pi*b*sgn(f) - pi*b) 
^2) - (pi*b*sgn(f) - pi*b)*sin(-1/2*pi*b*x*sgn(f) + 1/2*pi*b*x - 1/2*pi*a* 
sgn(f) + 1/2*pi*a)/(4*b^2*log(abs(f))^2 + (pi*b*sgn(f) - pi*b)^2))*e^(b*x* 
log(abs(f)) + a*log(abs(f))) + I*(-I*e^(1/2*I*pi*b*x*sgn(f) - 1/2*I*pi*b*x 
 + 1/2*I*pi*a*sgn(f) - 1/2*I*pi*a)/(2*I*pi*b*sgn(f) - 2*I*pi*b + 4*b*log(a 
bs(f))) + I*e^(-1/2*I*pi*b*x*sgn(f) + 1/2*I*pi*b*x - 1/2*I*pi*a*sgn(f) + 1 
/2*I*pi*a)/(-2*I*pi*b*sgn(f) + 2*I*pi*b + 4*b*log(abs(f))))*e^(b*x*log(abs 
(f)) + a*log(abs(f)))
 
3.4.46.9 Mupad [F(-1)]

Timed out. \[ \int f^{a+b x} \sinh ^2\left (d+e x+f x^2\right ) \, dx=\int f^{a+b\,x}\,{\mathrm {sinh}\left (f\,x^2+e\,x+d\right )}^2 \,d x \]

input
int(f^(a + b*x)*sinh(d + e*x + f*x^2)^2,x)
 
output
int(f^(a + b*x)*sinh(d + e*x + f*x^2)^2, x)