Integrand size = 14, antiderivative size = 94 \[ \int x^3 \sqrt {W\left (\frac {a}{x}\right )} \, dx=-\frac {256}{105} a^4 \sqrt {\pi } \text {erf}\left (2 \sqrt {W\left (\frac {a}{x}\right )}\right )+\frac {2}{7} x^4 \sqrt {W\left (\frac {a}{x}\right )}-\frac {2}{35} x^4 W\left (\frac {a}{x}\right )^{3/2}+\frac {16}{105} x^4 W\left (\frac {a}{x}\right )^{5/2}-\frac {128}{105} x^4 W\left (\frac {a}{x}\right )^{7/2} \] Output:
-256/105*a^4*Pi^(1/2)*erf(2*LambertW(a/x)^(1/2))+2/7*x^4*LambertW(a/x)^(1/ 2)-2/35*x^4*LambertW(a/x)^(3/2)+16/105*x^4*LambertW(a/x)^(5/2)-128/105*x^4 *LambertW(a/x)^(7/2)
Time = 0.01 (sec) , antiderivative size = 88, normalized size of antiderivative = 0.94 \[ \int x^3 \sqrt {W\left (\frac {a}{x}\right )} \, dx=-\frac {2}{105} \left (128 a^4 \sqrt {\pi } \text {erf}\left (2 \sqrt {W\left (\frac {a}{x}\right )}\right )-15 x^4 \sqrt {W\left (\frac {a}{x}\right )}+3 x^4 W\left (\frac {a}{x}\right )^{3/2}-8 x^4 W\left (\frac {a}{x}\right )^{5/2}+64 x^4 W\left (\frac {a}{x}\right )^{7/2}\right ) \] Input:
Integrate[x^3*Sqrt[ProductLog[a/x]],x]
Output:
(-2*(128*a^4*Sqrt[Pi]*Erf[2*Sqrt[ProductLog[a/x]]] - 15*x^4*Sqrt[ProductLo g[a/x]] + 3*x^4*ProductLog[a/x]^(3/2) - 8*x^4*ProductLog[a/x]^(5/2) + 64*x ^4*ProductLog[a/x]^(7/2)))/105
Time = 0.89 (sec) , antiderivative size = 105, normalized size of antiderivative = 1.12, number of steps used = 5, number of rules used = 5, \(\frac {\text {number of rules}}{\text {integrand size}}\) = 0.357, Rules used = {7173, 7206, 7206, 7206, 7203}
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^3 \sqrt {W\left (\frac {a}{x}\right )} \, dx\) |
\(\Big \downarrow \) 7173 |
\(\displaystyle \frac {2}{7} x^4 \sqrt {W\left (\frac {a}{x}\right )}-\frac {1}{7} \int \frac {x^3 W\left (\frac {a}{x}\right )^{3/2}}{W\left (\frac {a}{x}\right )+1}dx\) |
\(\Big \downarrow \) 7206 |
\(\displaystyle \frac {1}{7} \left (\frac {8}{5} \int \frac {x^3 W\left (\frac {a}{x}\right )^{5/2}}{W\left (\frac {a}{x}\right )+1}dx-\frac {2}{5} x^4 W\left (\frac {a}{x}\right )^{3/2}\right )+\frac {2}{7} x^4 \sqrt {W\left (\frac {a}{x}\right )}\) |
\(\Big \downarrow \) 7206 |
\(\displaystyle \frac {1}{7} \left (\frac {8}{5} \left (\frac {2}{3} x^4 W\left (\frac {a}{x}\right )^{5/2}-\frac {8}{3} \int \frac {x^3 W\left (\frac {a}{x}\right )^{7/2}}{W\left (\frac {a}{x}\right )+1}dx\right )-\frac {2}{5} x^4 W\left (\frac {a}{x}\right )^{3/2}\right )+\frac {2}{7} x^4 \sqrt {W\left (\frac {a}{x}\right )}\) |
\(\Big \downarrow \) 7206 |
\(\displaystyle \frac {1}{7} \left (\frac {8}{5} \left (\frac {2}{3} x^4 W\left (\frac {a}{x}\right )^{5/2}-\frac {8}{3} \left (2 x^4 W\left (\frac {a}{x}\right )^{7/2}-8 \int \frac {x^3 W\left (\frac {a}{x}\right )^{9/2}}{W\left (\frac {a}{x}\right )+1}dx\right )\right )-\frac {2}{5} x^4 W\left (\frac {a}{x}\right )^{3/2}\right )+\frac {2}{7} x^4 \sqrt {W\left (\frac {a}{x}\right )}\) |
\(\Big \downarrow \) 7203 |
\(\displaystyle \frac {1}{7} \left (\frac {8}{5} \left (\frac {2}{3} x^4 W\left (\frac {a}{x}\right )^{5/2}-\frac {8}{3} \left (4 \sqrt {\pi } a^4 \text {erf}\left (2 \sqrt {W\left (\frac {a}{x}\right )}\right )+2 x^4 W\left (\frac {a}{x}\right )^{7/2}\right )\right )-\frac {2}{5} x^4 W\left (\frac {a}{x}\right )^{3/2}\right )+\frac {2}{7} x^4 \sqrt {W\left (\frac {a}{x}\right )}\) |
Input:
Int[x^3*Sqrt[ProductLog[a/x]],x]
Output:
(2*x^4*Sqrt[ProductLog[a/x]])/7 + ((-2*x^4*ProductLog[a/x]^(3/2))/5 + (8*( (2*x^4*ProductLog[a/x]^(5/2))/3 - (8*(4*a^4*Sqrt[Pi]*Erf[2*Sqrt[ProductLog [a/x]]] + 2*x^4*ProductLog[a/x]^(7/2)))/3))/5)/7
Int[(x_)^(m_.)*((c_.)*ProductLog[(a_.)*(x_)^(n_.)])^(p_.), x_Symbol] :> Sim p[x^(m + 1)*((c*ProductLog[a*x^n])^p/(m + n*p + 1)), x] + Simp[n*(p/(c*(m + n*p + 1))) Int[x^m*((c*ProductLog[a*x^n])^(p + 1)/(1 + ProductLog[a*x^n] )), x], x] /; FreeQ[{a, c, m, n, p}, x] && (EqQ[m, -1] || (IntegerQ[p - 1/2 ] && ILtQ[Simplify[p + (m + 1)/n] - 1/2, 0]) || ( !IntegerQ[p - 1/2] && ILt Q[Simplify[p + (m + 1)/n], 0]))
Int[((x_)^(m_.)*((c_.)*ProductLog[(a_.)*(x_)^(n_.)])^(p_))/((d_) + (d_.)*Pr oductLog[(a_.)*(x_)^(n_.)]), x_Symbol] :> Simp[a^(p - 1/2)*c^(p - 1/2)*Rt[P i*(c/(p - 1/2)), 2]*(Erf[Sqrt[c*ProductLog[a*x^n]]/Rt[c/(p - 1/2), 2]]/(d*n )), x] /; FreeQ[{a, c, d, m, n}, x] && NeQ[m, -1] && IntegerQ[p - 1/2] && E qQ[m + n*(p - 1/2), -1] && PosQ[c/(p - 1/2)]
Int[((x_)^(m_.)*((c_.)*ProductLog[(a_.)*(x_)^(n_.)])^(p_.))/((d_) + (d_.)*P roductLog[(a_.)*(x_)^(n_.)]), x_Symbol] :> Simp[x^(m + 1)*((c*ProductLog[a* x^n])^p/(d*(m + n*p + 1))), x] - Simp[(m + 1)/(c*(m + n*p + 1)) Int[x^m*( (c*ProductLog[a*x^n])^(p + 1)/(d + d*ProductLog[a*x^n])), x], x] /; FreeQ[{ a, c, d, m, n, p}, x] && NeQ[m, -1] && LtQ[Simplify[p + (m + 1)/n], 0]
\[\int x^{3} \sqrt {\operatorname {LambertW}\left (\frac {a}{x}\right )}d x\]
Input:
int(x^3*LambertW(a/x)^(1/2),x)
Output:
int(x^3*LambertW(a/x)^(1/2),x)
\[ \int x^3 \sqrt {W\left (\frac {a}{x}\right )} \, dx=\int { x^{3} \sqrt {\operatorname {W}({\frac {a}{x}})} \,d x } \] Input:
integrate(x^3*lambert_w(a/x)^(1/2),x, algorithm="fricas")
Output:
integral(x^3*sqrt(lambert_w(a/x)), x)
\[ \int x^3 \sqrt {W\left (\frac {a}{x}\right )} \, dx=\int x^{3} \sqrt {W\left (\frac {a}{x}\right )}\, dx \] Input:
integrate(x**3*LambertW(a/x)**(1/2),x)
Output:
Integral(x**3*sqrt(LambertW(a/x)), x)
\[ \int x^3 \sqrt {W\left (\frac {a}{x}\right )} \, dx=\int { x^{3} \sqrt {\operatorname {W}({\frac {a}{x}})} \,d x } \] Input:
integrate(x^3*lambert_w(a/x)^(1/2),x, algorithm="maxima")
Output:
integrate(x^3*sqrt(lambert_w(a/x)), x)
\[ \int x^3 \sqrt {W\left (\frac {a}{x}\right )} \, dx=\int { x^{3} \sqrt {\operatorname {W}({\frac {a}{x}})} \,d x } \] Input:
integrate(x^3*lambert_w(a/x)^(1/2),x, algorithm="giac")
Output:
integrate(x^3*sqrt(lambert_w(a/x)), x)
Timed out. \[ \int x^3 \sqrt {W\left (\frac {a}{x}\right )} \, dx=\int x^3\,\sqrt {\mathrm {LambertW}\left (\frac {a}{x}\right )} \,d x \] Input:
int(x^3*LambertW(a/x)^(1/2),x)
Output:
int(x^3*LambertW(a/x)^(1/2), x)
\[ \int x^3 \sqrt {W\left (\frac {a}{x}\right )} \, dx=\frac {\sqrt {\textit {lambert\_w} \left (\frac {a}{x}\right )}\, x^{4}}{4}+\frac {\left (\int \frac {\sqrt {\textit {lambert\_w} \left (\frac {a}{x}\right )}\, x^{2}}{e^{\textit {lambert\_w} \left (\frac {a}{x}\right )} \textit {lambert\_w} \left (\frac {a}{x}\right )^{2}+e^{\textit {lambert\_w} \left (\frac {a}{x}\right )} \textit {lambert\_w} \left (\frac {a}{x}\right )}d x \right ) a}{8} \] Input:
int(x^3*Lambert_W(a/x)^(1/2),x)
Output:
(2*sqrt(lambert_w(a/x))*x**4 + int((sqrt(lambert_w(a/x))*x**2)/(e**lambert _w(a/x)*lambert_w(a/x)**2 + e**lambert_w(a/x)*lambert_w(a/x)),x)*a)/8