\(\int x^2 \text {erf}(d (a+b \log (c x^n))) \, dx\) [40]

Optimal result

Integrand size = 17, antiderivative size = 102 \[ \int x^2 \text {erf}\left (d \left (a+b \log \left (c x^n\right )\right )\right ) \, dx=\frac {1}{3} x^3 \text {erf}\left (d \left (a+b \log \left (c x^n\right )\right )\right )-\frac {1}{3} e^{\frac {9-12 a b d^2 n}{4 b^2 d^2 n^2}} x^3 \left (c x^n\right )^{-3/n} \text {erf}\left (\frac {2 a b d^2-\frac {3}{n}+2 b^2 d^2 \log \left (c x^n\right )}{2 b d}\right ) \] Output:

1/3*x^3*erf(d*(a+b*ln(c*x^n)))-1/3*exp(1/4*(-12*a*b*d^2*n+9)/b^2/d^2/n^2)* 
x^3*erf(1/2*(2*a*b*d^2-3/n+2*b^2*d^2*ln(c*x^n))/b/d)/((c*x^n)^(3/n))
 

Mathematica [A] (verified)

Time = 0.26 (sec) , antiderivative size = 88, normalized size of antiderivative = 0.86 \[ \int x^2 \text {erf}\left (d \left (a+b \log \left (c x^n\right )\right )\right ) \, dx=\frac {1}{3} \left (x^3 \text {erf}\left (d \left (a+b \log \left (c x^n\right )\right )\right )-e^{-\frac {3 \left (\frac {-\frac {3}{d^2}+4 a b n}{b^2}+4 n \log \left (c x^n\right )\right )}{4 n^2}} x^3 \text {erf}\left (a d-\frac {3}{2 b d n}+b d \log \left (c x^n\right )\right )\right ) \] Input:

Integrate[x^2*Erf[d*(a + b*Log[c*x^n])],x]
 

Output:

(x^3*Erf[d*(a + b*Log[c*x^n])] - (x^3*Erf[a*d - 3/(2*b*d*n) + b*d*Log[c*x^ 
n]])/E^((3*((-3/d^2 + 4*a*b*n)/b^2 + 4*n*Log[c*x^n]))/(4*n^2)))/3
 

Rubi [A] (verified)

Time = 0.55 (sec) , antiderivative size = 102, normalized size of antiderivative = 1.00, number of steps used = 6, number of rules used = 5, \(\frac {\text {number of rules}}{\text {integrand size}}\) = 0.294, Rules used = {6955, 2712, 2706, 2664, 2634}

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^2 \text {erf}\left (d \left (a+b \log \left (c x^n\right )\right )\right ) \, dx\)

\(\Big \downarrow \) 6955

\(\displaystyle \frac {1}{3} x^3 \text {erf}\left (d \left (a+b \log \left (c x^n\right )\right )\right )-\frac {2 b d n \int e^{-d^2 \left (a+b \log \left (c x^n\right )\right )^2} x^2dx}{3 \sqrt {\pi }}\)

\(\Big \downarrow \) 2712

\(\displaystyle \frac {1}{3} x^3 \text {erf}\left (d \left (a+b \log \left (c x^n\right )\right )\right )-\frac {2 b d n x^{2 a b d^2 n} \left (c x^n\right )^{-2 a b d^2} \int e^{-a^2 d^2-b^2 \log ^2\left (c x^n\right ) d^2} x^{2-2 a b d^2 n}dx}{3 \sqrt {\pi }}\)

\(\Big \downarrow \) 2706

\(\displaystyle \frac {1}{3} x^3 \text {erf}\left (d \left (a+b \log \left (c x^n\right )\right )\right )-\frac {2 b d x^3 \left (c x^n\right )^{-3/n} \int \exp \left (-a^2 d^2-b^2 \log ^2\left (c x^n\right ) d^2+\frac {\left (3-2 a b d^2 n\right ) \log \left (c x^n\right )}{n}\right )d\log \left (c x^n\right )}{3 \sqrt {\pi }}\)

\(\Big \downarrow \) 2664

\(\displaystyle \frac {1}{3} x^3 \text {erf}\left (d \left (a+b \log \left (c x^n\right )\right )\right )-\frac {2 b d x^3 \left (c x^n\right )^{-3/n} e^{\frac {9-12 a b d^2 n}{4 b^2 d^2 n^2}} \int \exp \left (-\frac {\left (2 a b d^2+2 b^2 \log \left (c x^n\right ) d^2-\frac {3}{n}\right )^2}{4 b^2 d^2}\right )d\log \left (c x^n\right )}{3 \sqrt {\pi }}\)

\(\Big \downarrow \) 2634

\(\displaystyle \frac {1}{3} x^3 \text {erf}\left (d \left (a+b \log \left (c x^n\right )\right )\right )-\frac {1}{3} x^3 \left (c x^n\right )^{-3/n} e^{\frac {9-12 a b d^2 n}{4 b^2 d^2 n^2}} \text {erf}\left (\frac {2 a b d^2+2 b^2 d^2 \log \left (c x^n\right )-\frac {3}{n}}{2 b d}\right )\)

Input:

Int[x^2*Erf[d*(a + b*Log[c*x^n])],x]
 

Output:

(x^3*Erf[d*(a + b*Log[c*x^n])])/3 - (E^((9 - 12*a*b*d^2*n)/(4*b^2*d^2*n^2) 
)*x^3*Erf[(2*a*b*d^2 - 3/n + 2*b^2*d^2*Log[c*x^n])/(2*b*d)])/(3*(c*x^n)^(3 
/n))
 

Defintions of rubi rules used

rule 2634

Int[(F_)^((a_.) + (b_.)*((c_.) + (d_.)*(x_))^2), x_Symbol] :> Simp[F^a*Sqrt 
[Pi]*(Erf[(c + d*x)*Rt[(-b)*Log[F], 2]]/(2*d*Rt[(-b)*Log[F], 2])), x] /; Fr 
eeQ[{F, a, b, c, d}, x] && NegQ[b]
 

rule 2664

Int[(F_)^((a_.) + (b_.)*(x_) + (c_.)*(x_)^2), x_Symbol] :> Simp[F^(a - b^2/ 
(4*c))   Int[F^((b + 2*c*x)^2/(4*c)), x], x] /; FreeQ[{F, a, b, c}, x]
 

rule 2706

Int[(F_)^(((a_.) + Log[(c_.)*((d_.) + (e_.)*(x_))^(n_.)]^2*(b_.))*(f_.))*(( 
g_.) + (h_.)*(x_))^(m_.), x_Symbol] :> Simp[(g + h*x)^(m + 1)/(h*n*(c*(d + 
e*x)^n)^((m + 1)/n))   Subst[Int[E^(a*f*Log[F] + ((m + 1)*x)/n + b*f*Log[F] 
*x^2), x], x, Log[c*(d + e*x)^n]], x] /; FreeQ[{F, a, b, c, d, e, f, g, h, 
m, n}, x] && EqQ[e*g - d*h, 0]
 

rule 2712

Int[(F_)^(((a_.) + Log[(c_.)*((d_.) + (e_.)*(x_))^(n_.)]*(b_.))^2*(f_.))*(( 
g_.) + (h_.)*(x_))^(m_.), x_Symbol] :> Simp[(g + h*x)^m*((c*(d + e*x)^n)^(2 
*a*b*f*Log[F])/(d + e*x)^(m + 2*a*b*f*n*Log[F]))*Int[(d + e*x)^(m + 2*a*b*f 
*n*Log[F])*F^(a^2*f + b^2*f*Log[c*(d + e*x)^n]^2), x], x] /; FreeQ[{F, a, b 
, c, d, e, f, g, h, m, n}, x] && EqQ[e*g - d*h, 0]
 

rule 6955

Int[Erf[((a_.) + Log[(c_.)*(x_)^(n_.)]*(b_.))*(d_.)]*((e_.)*(x_))^(m_.), x_ 
Symbol] :> Simp[(e*x)^(m + 1)*(Erf[d*(a + b*Log[c*x^n])]/(e*(m + 1))), x] - 
 Simp[2*b*d*(n/(Sqrt[Pi]*(m + 1)))   Int[(e*x)^m/E^(d*(a + b*Log[c*x^n]))^2 
, x], x] /; FreeQ[{a, b, c, d, e, m, n}, x] && NeQ[m, -1]
 

Maple [F]

\[\int x^{2} \operatorname {erf}\left (d \left (a +b \ln \left (c \,x^{n}\right )\right )\right )d x\]

Input:

int(x^2*erf(d*(a+b*ln(c*x^n))),x)
 

Output:

int(x^2*erf(d*(a+b*ln(c*x^n))),x)
 

Fricas [A] (verification not implemented)

Time = 0.09 (sec) , antiderivative size = 125, normalized size of antiderivative = 1.23 \[ \int x^2 \text {erf}\left (d \left (a+b \log \left (c x^n\right )\right )\right ) \, dx=\frac {1}{3} \, x^{3} \operatorname {erf}\left (b d \log \left (c x^{n}\right ) + a d\right ) - \frac {1}{3} \, \sqrt {b^{2} d^{2} n^{2}} \operatorname {erf}\left (\frac {{\left (2 \, b^{2} d^{2} n^{2} \log \left (x\right ) + 2 \, b^{2} d^{2} n \log \left (c\right ) + 2 \, a b d^{2} n - 3\right )} \sqrt {b^{2} d^{2} n^{2}}}{2 \, b^{2} d^{2} n^{2}}\right ) e^{\left (-\frac {3 \, {\left (4 \, b^{2} d^{2} n \log \left (c\right ) + 4 \, a b d^{2} n - 3\right )}}{4 \, b^{2} d^{2} n^{2}}\right )} \] Input:

integrate(x^2*erf(d*(a+b*log(c*x^n))),x, algorithm="fricas")
 

Output:

1/3*x^3*erf(b*d*log(c*x^n) + a*d) - 1/3*sqrt(b^2*d^2*n^2)*erf(1/2*(2*b^2*d 
^2*n^2*log(x) + 2*b^2*d^2*n*log(c) + 2*a*b*d^2*n - 3)*sqrt(b^2*d^2*n^2)/(b 
^2*d^2*n^2))*e^(-3/4*(4*b^2*d^2*n*log(c) + 4*a*b*d^2*n - 3)/(b^2*d^2*n^2))
 

Sympy [F]

\[ \int x^2 \text {erf}\left (d \left (a+b \log \left (c x^n\right )\right )\right ) \, dx=\int x^{2} \operatorname {erf}{\left (a d + b d \log {\left (c x^{n} \right )} \right )}\, dx \] Input:

integrate(x**2*erf(d*(a+b*ln(c*x**n))),x)
 

Output:

Integral(x**2*erf(a*d + b*d*log(c*x**n)), x)
 

Maxima [F]

\[ \int x^2 \text {erf}\left (d \left (a+b \log \left (c x^n\right )\right )\right ) \, dx=\int { x^{2} \operatorname {erf}\left ({\left (b \log \left (c x^{n}\right ) + a\right )} d\right ) \,d x } \] Input:

integrate(x^2*erf(d*(a+b*log(c*x^n))),x, algorithm="maxima")
 

Output:

1/3*x^3*erf(b*d*log(x^n) + (b*log(c) + a)*d) - 2/3*b*d*n*integrate(x^2*e^( 
-b^2*d^2*log(c)^2 - 2*b^2*d^2*log(c)*log(x^n) - b^2*d^2*log(x^n)^2 - 2*a*b 
*d^2*log(x^n) - a^2*d^2), x)/(sqrt(pi)*c^(2*a*b*d^2))
 

Giac [A] (verification not implemented)

Time = 0.28 (sec) , antiderivative size = 85, normalized size of antiderivative = 0.83 \[ \int x^2 \text {erf}\left (d \left (a+b \log \left (c x^n\right )\right )\right ) \, dx=\frac {1}{3} \, x^{3} \operatorname {erf}\left (b d n \log \left (x\right ) + b d \log \left (c\right ) + a d\right ) + \frac {\operatorname {erf}\left (-b d n \log \left (x\right ) - b d \log \left (c\right ) - a d + \frac {3}{2 \, b d n}\right ) e^{\left (-\frac {3 \, a}{b n} + \frac {9}{4 \, b^{2} d^{2} n^{2}}\right )}}{3 \, c^{\frac {3}{n}}} \] Input:

integrate(x^2*erf(d*(a+b*log(c*x^n))),x, algorithm="giac")
 

Output:

1/3*x^3*erf(b*d*n*log(x) + b*d*log(c) + a*d) + 1/3*erf(-b*d*n*log(x) - b*d 
*log(c) - a*d + 3/2/(b*d*n))*e^(-3*a/(b*n) + 9/4/(b^2*d^2*n^2))/c^(3/n)
 

Mupad [F(-1)]

Timed out. \[ \int x^2 \text {erf}\left (d \left (a+b \log \left (c x^n\right )\right )\right ) \, dx=\int x^2\,\mathrm {erf}\left (d\,\left (a+b\,\ln \left (c\,x^n\right )\right )\right ) \,d x \] Input:

int(x^2*erf(d*(a + b*log(c*x^n))),x)
 

Output:

int(x^2*erf(d*(a + b*log(c*x^n))), x)
 

Reduce [F]

\[ \int x^2 \text {erf}\left (d \left (a+b \log \left (c x^n\right )\right )\right ) \, dx=\int \mathrm {erf}\left (\mathrm {log}\left (x^{n} c \right ) b d +a d \right ) x^{2}d x \] Input:

int(x^2*erf(d*(a+b*log(c*x^n))),x)
 

Output:

int(erf(log(x**n*c)*b*d + a*d)*x**2,x)