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

Optimal result

Integrand size = 20, antiderivative size = 132 \[ \int \frac {e^{a+b x} x^2}{c+d x^2} \, dx=\frac {e^{a+b x}}{b d}+\frac {\sqrt {-c} e^{a+\frac {b \sqrt {-c}}{\sqrt {d}}} \operatorname {ExpIntegralEi}\left (-\frac {b \left (\sqrt {-c}-\sqrt {d} x\right )}{\sqrt {d}}\right )}{2 d^{3/2}}-\frac {\sqrt {-c} e^{a-\frac {b \sqrt {-c}}{\sqrt {d}}} \operatorname {ExpIntegralEi}\left (\frac {b \left (\sqrt {-c}+\sqrt {d} x\right )}{\sqrt {d}}\right )}{2 d^{3/2}} \]

[Out]

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

Rubi [A] (verified)

Time = 0.14 (sec) , antiderivative size = 132, normalized size of antiderivative = 1.00, number of steps used = 7, number of rules used = 4, \(\frac {\text {number of rules}}{\text {integrand size}}\) = 0.200, Rules used = {2303, 2225, 2301, 2209} \[ \int \frac {e^{a+b x} x^2}{c+d x^2} \, dx=\frac {\sqrt {-c} e^{a+\frac {b \sqrt {-c}}{\sqrt {d}}} \operatorname {ExpIntegralEi}\left (-\frac {b \left (\sqrt {-c}-\sqrt {d} x\right )}{\sqrt {d}}\right )}{2 d^{3/2}}-\frac {\sqrt {-c} e^{a-\frac {b \sqrt {-c}}{\sqrt {d}}} \operatorname {ExpIntegralEi}\left (\frac {b \left (\sqrt {d} x+\sqrt {-c}\right )}{\sqrt {d}}\right )}{2 d^{3/2}}+\frac {e^{a+b x}}{b d} \]

[In]

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

[Out]

E^(a + b*x)/(b*d) + (Sqrt[-c]*E^(a + (b*Sqrt[-c])/Sqrt[d])*ExpIntegralEi[-((b*(Sqrt[-c] - Sqrt[d]*x))/Sqrt[d])
])/(2*d^(3/2)) - (Sqrt[-c]*E^(a - (b*Sqrt[-c])/Sqrt[d])*ExpIntegralEi[(b*(Sqrt[-c] + Sqrt[d]*x))/Sqrt[d]])/(2*
d^(3/2))

Rule 2209

Int[(F_)^((g_.)*((e_.) + (f_.)*(x_)))/((c_.) + (d_.)*(x_)), x_Symbol] :> Simp[(F^(g*(e - c*(f/d)))/d)*ExpInteg
ralEi[f*g*(c + d*x)*(Log[F]/d)], x] /; FreeQ[{F, c, d, e, f, g}, x] &&  !TrueQ[$UseGamma]

Rule 2225

Int[((F_)^((c_.)*((a_.) + (b_.)*(x_))))^(n_.), x_Symbol] :> Simp[(F^(c*(a + b*x)))^n/(b*c*n*Log[F]), x] /; Fre
eQ[{F, a, b, c, n}, x]

Rule 2301

Int[(F_)^((g_.)*((d_.) + (e_.)*(x_))^(n_.))/((a_) + (c_.)*(x_)^2), x_Symbol] :> Int[ExpandIntegrand[F^(g*(d +
e*x)^n), 1/(a + c*x^2), x], x] /; FreeQ[{F, a, c, d, e, g, n}, x]

Rule 2303

Int[((F_)^((g_.)*((d_.) + (e_.)*(x_))^(n_.))*(u_)^(m_.))/((a_) + (c_)*(x_)^2), x_Symbol] :> Int[ExpandIntegran
d[F^(g*(d + e*x)^n), u^m/(a + c*x^2), x], x] /; FreeQ[{F, a, c, d, e, g, n}, x] && PolynomialQ[u, x] && Intege
rQ[m]

Rubi steps \begin{align*} \text {integral}& = \int \left (\frac {e^{a+b x}}{d}-\frac {c e^{a+b x}}{d \left (c+d x^2\right )}\right ) \, dx \\ & = \frac {\int e^{a+b x} \, dx}{d}-\frac {c \int \frac {e^{a+b x}}{c+d x^2} \, dx}{d} \\ & = \frac {e^{a+b x}}{b d}-\frac {c \int \left (\frac {\sqrt {-c} e^{a+b x}}{2 c \left (\sqrt {-c}-\sqrt {d} x\right )}+\frac {\sqrt {-c} e^{a+b x}}{2 c \left (\sqrt {-c}+\sqrt {d} x\right )}\right ) \, dx}{d} \\ & = \frac {e^{a+b x}}{b d}-\frac {\sqrt {-c} \int \frac {e^{a+b x}}{\sqrt {-c}-\sqrt {d} x} \, dx}{2 d}-\frac {\sqrt {-c} \int \frac {e^{a+b x}}{\sqrt {-c}+\sqrt {d} x} \, dx}{2 d} \\ & = \frac {e^{a+b x}}{b d}+\frac {\sqrt {-c} e^{a+\frac {b \sqrt {-c}}{\sqrt {d}}} \text {Ei}\left (-\frac {b \left (\sqrt {-c}-\sqrt {d} x\right )}{\sqrt {d}}\right )}{2 d^{3/2}}-\frac {\sqrt {-c} e^{a-\frac {b \sqrt {-c}}{\sqrt {d}}} \text {Ei}\left (\frac {b \left (\sqrt {-c}+\sqrt {d} x\right )}{\sqrt {d}}\right )}{2 d^{3/2}} \\ \end{align*}

Mathematica [C] (verified)

Result contains complex when optimal does not.

Time = 0.14 (sec) , antiderivative size = 120, normalized size of antiderivative = 0.91 \[ \int \frac {e^{a+b x} x^2}{c+d x^2} \, dx=\frac {e^a \left (2 \sqrt {d} e^{b x}+i b \sqrt {c} e^{\frac {i b \sqrt {c}}{\sqrt {d}}} \operatorname {ExpIntegralEi}\left (b \left (-\frac {i \sqrt {c}}{\sqrt {d}}+x\right )\right )-i b \sqrt {c} e^{-\frac {i b \sqrt {c}}{\sqrt {d}}} \operatorname {ExpIntegralEi}\left (b \left (\frac {i \sqrt {c}}{\sqrt {d}}+x\right )\right )\right )}{2 b d^{3/2}} \]

[In]

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

[Out]

(E^a*(2*Sqrt[d]*E^(b*x) + I*b*Sqrt[c]*E^((I*b*Sqrt[c])/Sqrt[d])*ExpIntegralEi[b*(((-I)*Sqrt[c])/Sqrt[d] + x)]
- (I*b*Sqrt[c]*ExpIntegralEi[b*((I*Sqrt[c])/Sqrt[d] + x)])/E^((I*b*Sqrt[c])/Sqrt[d])))/(2*b*d^(3/2))

Maple [A] (verified)

Time = 0.44 (sec) , antiderivative size = 127, normalized size of antiderivative = 0.96

method	result	size
risch	\(\frac {{\mathrm e}^{\frac {b \sqrt {-c d}+a d}{d}} \operatorname {Ei}_{1}\left (\frac {b \sqrt {-c d}+a d -d \left (b x +a \right )}{d}\right ) c}{2 d \sqrt {-c d}}-\frac {{\mathrm e}^{\frac {-b \sqrt {-c d}+a d}{d}} \operatorname {Ei}_{1}\left (-\frac {b \sqrt {-c d}-a d +d \left (b x +a \right )}{d}\right ) c}{2 d \sqrt {-c d}}+\frac {{\mathrm e}^{b x +a}}{b d}\)	\(127\)
derivativedivides	\(\frac {-\frac {a^{2} b \left ({\mathrm e}^{\frac {b \sqrt {-c d}+a d}{d}} \operatorname {Ei}_{1}\left (\frac {b \sqrt {-c d}+a d -d \left (b x +a \right )}{d}\right )-{\mathrm e}^{-\frac {b \sqrt {-c d}-a d}{d}} \operatorname {Ei}_{1}\left (-\frac {b \sqrt {-c d}-a d +d \left (b x +a \right )}{d}\right )\right )}{2 \sqrt {-c d}}+\frac {b^{2} {\mathrm e}^{b x +a}}{d}-\frac {b \left (2 \,{\mathrm e}^{\frac {b \sqrt {-c d}+a d}{d}} \operatorname {Ei}_{1}\left (\frac {b \sqrt {-c d}+a d -d \left (b x +a \right )}{d}\right ) \sqrt {-c d}\, a b +{\mathrm e}^{\frac {b \sqrt {-c d}+a d}{d}} \operatorname {Ei}_{1}\left (\frac {b \sqrt {-c d}+a d -d \left (b x +a \right )}{d}\right ) a^{2} d -{\mathrm e}^{\frac {b \sqrt {-c d}+a d}{d}} \operatorname {Ei}_{1}\left (\frac {b \sqrt {-c d}+a d -d \left (b x +a \right )}{d}\right ) b^{2} c +2 \,{\mathrm e}^{-\frac {b \sqrt {-c d}-a d}{d}} \operatorname {Ei}_{1}\left (-\frac {b \sqrt {-c d}-a d +d \left (b x +a \right )}{d}\right ) \sqrt {-c d}\, a b -{\mathrm e}^{-\frac {b \sqrt {-c d}-a d}{d}} \operatorname {Ei}_{1}\left (-\frac {b \sqrt {-c d}-a d +d \left (b x +a \right )}{d}\right ) a^{2} d +{\mathrm e}^{-\frac {b \sqrt {-c d}-a d}{d}} \operatorname {Ei}_{1}\left (-\frac {b \sqrt {-c d}-a d +d \left (b x +a \right )}{d}\right ) b^{2} c \right )}{2 d \sqrt {-c d}}+\frac {a b \left (\sqrt {-c d}\, {\mathrm e}^{\frac {b \sqrt {-c d}+a d}{d}} \operatorname {Ei}_{1}\left (\frac {b \sqrt {-c d}+a d -d \left (b x +a \right )}{d}\right ) b +\sqrt {-c d}\, {\mathrm e}^{-\frac {b \sqrt {-c d}-a d}{d}} \operatorname {Ei}_{1}\left (-\frac {b \sqrt {-c d}-a d +d \left (b x +a \right )}{d}\right ) b +{\mathrm e}^{\frac {b \sqrt {-c d}+a d}{d}} \operatorname {Ei}_{1}\left (\frac {b \sqrt {-c d}+a d -d \left (b x +a \right )}{d}\right ) a d -{\mathrm e}^{-\frac {b \sqrt {-c d}-a d}{d}} \operatorname {Ei}_{1}\left (-\frac {b \sqrt {-c d}-a d +d \left (b x +a \right )}{d}\right ) a d \right )}{d \sqrt {-c d}}}{b^{3}}\)	\(660\)
default	\(\frac {-\frac {a^{2} b \left ({\mathrm e}^{\frac {b \sqrt {-c d}+a d}{d}} \operatorname {Ei}_{1}\left (\frac {b \sqrt {-c d}+a d -d \left (b x +a \right )}{d}\right )-{\mathrm e}^{-\frac {b \sqrt {-c d}-a d}{d}} \operatorname {Ei}_{1}\left (-\frac {b \sqrt {-c d}-a d +d \left (b x +a \right )}{d}\right )\right )}{2 \sqrt {-c d}}+\frac {b^{2} {\mathrm e}^{b x +a}}{d}-\frac {b \left (2 \,{\mathrm e}^{\frac {b \sqrt {-c d}+a d}{d}} \operatorname {Ei}_{1}\left (\frac {b \sqrt {-c d}+a d -d \left (b x +a \right )}{d}\right ) \sqrt {-c d}\, a b +{\mathrm e}^{\frac {b \sqrt {-c d}+a d}{d}} \operatorname {Ei}_{1}\left (\frac {b \sqrt {-c d}+a d -d \left (b x +a \right )}{d}\right ) a^{2} d -{\mathrm e}^{\frac {b \sqrt {-c d}+a d}{d}} \operatorname {Ei}_{1}\left (\frac {b \sqrt {-c d}+a d -d \left (b x +a \right )}{d}\right ) b^{2} c +2 \,{\mathrm e}^{-\frac {b \sqrt {-c d}-a d}{d}} \operatorname {Ei}_{1}\left (-\frac {b \sqrt {-c d}-a d +d \left (b x +a \right )}{d}\right ) \sqrt {-c d}\, a b -{\mathrm e}^{-\frac {b \sqrt {-c d}-a d}{d}} \operatorname {Ei}_{1}\left (-\frac {b \sqrt {-c d}-a d +d \left (b x +a \right )}{d}\right ) a^{2} d +{\mathrm e}^{-\frac {b \sqrt {-c d}-a d}{d}} \operatorname {Ei}_{1}\left (-\frac {b \sqrt {-c d}-a d +d \left (b x +a \right )}{d}\right ) b^{2} c \right )}{2 d \sqrt {-c d}}+\frac {a b \left (\sqrt {-c d}\, {\mathrm e}^{\frac {b \sqrt {-c d}+a d}{d}} \operatorname {Ei}_{1}\left (\frac {b \sqrt {-c d}+a d -d \left (b x +a \right )}{d}\right ) b +\sqrt {-c d}\, {\mathrm e}^{-\frac {b \sqrt {-c d}-a d}{d}} \operatorname {Ei}_{1}\left (-\frac {b \sqrt {-c d}-a d +d \left (b x +a \right )}{d}\right ) b +{\mathrm e}^{\frac {b \sqrt {-c d}+a d}{d}} \operatorname {Ei}_{1}\left (\frac {b \sqrt {-c d}+a d -d \left (b x +a \right )}{d}\right ) a d -{\mathrm e}^{-\frac {b \sqrt {-c d}-a d}{d}} \operatorname {Ei}_{1}\left (-\frac {b \sqrt {-c d}-a d +d \left (b x +a \right )}{d}\right ) a d \right )}{d \sqrt {-c d}}}{b^{3}}\)	\(660\)

[In]

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

[Out]

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

Fricas [A] (verification not implemented)

none

Time = 0.34 (sec) , antiderivative size = 106, normalized size of antiderivative = 0.80 \[ \int \frac {e^{a+b x} x^2}{c+d x^2} \, dx=\frac {\sqrt {-\frac {b^{2} c}{d}} {\rm Ei}\left (b x - \sqrt {-\frac {b^{2} c}{d}}\right ) e^{\left (a + \sqrt {-\frac {b^{2} c}{d}}\right )} - \sqrt {-\frac {b^{2} c}{d}} {\rm Ei}\left (b x + \sqrt {-\frac {b^{2} c}{d}}\right ) e^{\left (a - \sqrt {-\frac {b^{2} c}{d}}\right )} + 2 \, e^{\left (b x + a\right )}}{2 \, b d} \]

[In]

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

[Out]

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

Sympy [F]

\[ \int \frac {e^{a+b x} x^2}{c+d x^2} \, dx=e^{a} \int \frac {x^{2} e^{b x}}{c + d x^{2}}\, dx \]

[In]

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

[Out]

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

Maxima [F]

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

[In]

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

[Out]

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

Giac [F]

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

[In]

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

[Out]

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

Mupad [F(-1)]

Timed out. \[ \int \frac {e^{a+b x} x^2}{c+d x^2} \, dx=\int \frac {x^2\,{\mathrm {e}}^{a+b\,x}}{d\,x^2+c} \,d x \]

[In]

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

[Out]

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