\(\int \frac {e^{c+d x^2} \text {erf}(b x)}{x^5} \, dx\) [58]

Optimal result

Integrand size = 17, antiderivative size = 17 \[ \int \frac {e^{c+d x^2} \text {erf}(b x)}{x^5} \, dx=-\frac {b e^{c-\left (b^2-d\right ) x^2}}{6 \sqrt {\pi } x^3}+\frac {b \left (b^2-d\right ) e^{c-\left (b^2-d\right ) x^2}}{3 \sqrt {\pi } x}-\frac {b d e^{c-\left (b^2-d\right ) x^2}}{2 \sqrt {\pi } x}-\frac {e^{c+d x^2} \text {erf}(b x)}{4 x^4}-\frac {d e^{c+d x^2} \text {erf}(b x)}{4 x^2}+\frac {1}{3} b \left (b^2-d\right )^{3/2} e^c \text {erf}\left (\sqrt {b^2-d} x\right )-\frac {1}{2} b \sqrt {b^2-d} d e^c \text {erf}\left (\sqrt {b^2-d} x\right )+\frac {1}{2} d^2 \text {Int}\left (\frac {e^{c+d x^2} \text {erf}(b x)}{x},x\right ) \]

[Out]

-1/4*exp(d*x^2+c)*erf(b*x)/x^4-1/4*d*exp(d*x^2+c)*erf(b*x)/x^2+1/3*b*(b^2-d)^(3/2)*exp(c)*erf(x*(b^2-d)^(1/2))
-1/2*b*d*exp(c)*erf(x*(b^2-d)^(1/2))*(b^2-d)^(1/2)-1/6*b*exp(c-(b^2-d)*x^2)/x^3/Pi^(1/2)+1/3*b*(b^2-d)*exp(c-(
b^2-d)*x^2)/x/Pi^(1/2)-1/2*b*d*exp(c-(b^2-d)*x^2)/x/Pi^(1/2)+1/2*d^2*Unintegrable(exp(d*x^2+c)*erf(b*x)/x,x)

Rubi [N/A]

Not integrable

Time = 0.27 (sec) , antiderivative size = 17, normalized size of antiderivative = 1.00, number of steps used = 0, number of rules used = 0, \(\frac {\text {number of rules}}{\text {integrand size}}\) = 0.000, Rules used = {} \[ \int \frac {e^{c+d x^2} \text {erf}(b x)}{x^5} \, dx=\int \frac {e^{c+d x^2} \text {erf}(b x)}{x^5} \, dx \]

[In]

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

[Out]

-1/6*(b*E^(c - (b^2 - d)*x^2))/(Sqrt[Pi]*x^3) + (b*(b^2 - d)*E^(c - (b^2 - d)*x^2))/(3*Sqrt[Pi]*x) - (b*d*E^(c
- (b^2 - d)*x^2))/(2*Sqrt[Pi]*x) - (E^(c + d*x^2)*Erf[b*x])/(4*x^4) - (d*E^(c + d*x^2)*Erf[b*x])/(4*x^2) + (b
*(b^2 - d)^(3/2)*E^c*Erf[Sqrt[b^2 - d]*x])/3 - (b*Sqrt[b^2 - d]*d*E^c*Erf[Sqrt[b^2 - d]*x])/2 + (d^2*Defer[Int
][(E^(c + d*x^2)*Erf[b*x])/x, x])/2

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

Mathematica [N/A]

Not integrable

Time = 0.17 (sec) , antiderivative size = 19, normalized size of antiderivative = 1.12 \[ \int \frac {e^{c+d x^2} \text {erf}(b x)}{x^5} \, dx=\int \frac {e^{c+d x^2} \text {erf}(b x)}{x^5} \, dx \]

[In]

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

[Out]

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

Maple [N/A] (verified)

Not integrable

Time = 0.18 (sec) , antiderivative size = 16, normalized size of antiderivative = 0.94

\[\int \frac {{\mathrm e}^{d \,x^{2}+c} \operatorname {erf}\left (b x \right )}{x^{5}}d x\]

[In]

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

[Out]

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

Fricas [N/A]

Not integrable

Time = 0.25 (sec) , antiderivative size = 18, normalized size of antiderivative = 1.06 \[ \int \frac {e^{c+d x^2} \text {erf}(b x)}{x^5} \, dx=\int { \frac {\operatorname {erf}\left (b x\right ) e^{\left (d x^{2} + c\right )}}{x^{5}} \,d x } \]

[In]

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

[Out]

integral(erf(b*x)*e^(d*x^2 + c)/x^5, x)

Sympy [N/A]

Not integrable

Time = 31.90 (sec) , antiderivative size = 19, normalized size of antiderivative = 1.12 \[ \int \frac {e^{c+d x^2} \text {erf}(b x)}{x^5} \, dx=e^{c} \int \frac {e^{d x^{2}} \operatorname {erf}{\left (b x \right )}}{x^{5}}\, dx \]

[In]

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

[Out]

exp(c)*Integral(exp(d*x**2)*erf(b*x)/x**5, x)

Maxima [N/A]

Not integrable

Time = 0.28 (sec) , antiderivative size = 18, normalized size of antiderivative = 1.06 \[ \int \frac {e^{c+d x^2} \text {erf}(b x)}{x^5} \, dx=\int { \frac {\operatorname {erf}\left (b x\right ) e^{\left (d x^{2} + c\right )}}{x^{5}} \,d x } \]

[In]

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

[Out]

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

Giac [N/A]

Not integrable

Time = 0.27 (sec) , antiderivative size = 18, normalized size of antiderivative = 1.06 \[ \int \frac {e^{c+d x^2} \text {erf}(b x)}{x^5} \, dx=\int { \frac {\operatorname {erf}\left (b x\right ) e^{\left (d x^{2} + c\right )}}{x^{5}} \,d x } \]

[In]

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

[Out]

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

Mupad [N/A]

Not integrable

Time = 6.12 (sec) , antiderivative size = 18, normalized size of antiderivative = 1.06 \[ \int \frac {e^{c+d x^2} \text {erf}(b x)}{x^5} \, dx=\int \frac {{\mathrm {e}}^{d\,x^2+c}\,\mathrm {erf}\left (b\,x\right )}{x^5} \,d x \]

[In]

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

[Out]

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