\(\int \frac {x^{20}}{\sqrt {1+x^6}} \, dx\) [615]

Optimal result

Integrand size = 13, antiderivative size = 48 \[ \int \frac {x^{20}}{\sqrt {1+x^6}} \, dx=\frac {1}{144} \sqrt {1+x^6} \left (15 x^3-10 x^9+8 x^{15}\right )-\frac {5}{48} \log \left (x^3+\sqrt {1+x^6}\right ) \]

[Out]

1/144*(x^6+1)^(1/2)*(8*x^15-10*x^9+15*x^3)-5/48*ln(x^3+(x^6+1)^(1/2))

Rubi [A] (verified)

Time = 0.02 (sec) , antiderivative size = 57, normalized size of antiderivative = 1.19, number of steps used = 5, number of rules used = 3, \(\frac {\text {number of rules}}{\text {integrand size}}\) = 0.231, Rules used = {281, 327, 221} \[ \int \frac {x^{20}}{\sqrt {1+x^6}} \, dx=-\frac {5 \text {arcsinh}\left (x^3\right )}{48}+\frac {1}{18} \sqrt {x^6+1} x^{15}-\frac {5}{72} \sqrt {x^6+1} x^9+\frac {5}{48} \sqrt {x^6+1} x^3 \]

[In]

Int[x^20/Sqrt[1 + x^6],x]

[Out]

(5*x^3*Sqrt[1 + x^6])/48 - (5*x^9*Sqrt[1 + x^6])/72 + (x^15*Sqrt[1 + x^6])/18 - (5*ArcSinh[x^3])/48

Rule 221

Int[1/Sqrt[(a_) + (b_.)*(x_)^2], x_Symbol] :> Simp[ArcSinh[Rt[b, 2]*(x/Sqrt[a])]/Rt[b, 2], x] /; FreeQ[{a, b},
x] && GtQ[a, 0] && PosQ[b]

Rule 281

Int[(x_)^(m_.)*((a_) + (b_.)*(x_)^(n_))^(p_), x_Symbol] :> With[{k = GCD[m + 1, n]}, Dist[1/k, Subst[Int[x^((m
+ 1)/k - 1)*(a + b*x^(n/k))^p, x], x, x^k], x] /; k != 1] /; FreeQ[{a, b, p}, x] && IGtQ[n, 0] && IntegerQ[m]

Rule 327

Int[((c_.)*(x_))^(m_)*((a_) + (b_.)*(x_)^(n_))^(p_), x_Symbol] :> Simp[c^(n - 1)*(c*x)^(m - n + 1)*((a + b*x^n
)^(p + 1)/(b*(m + n*p + 1))), x] - Dist[a*c^n*((m - n + 1)/(b*(m + n*p + 1))), Int[(c*x)^(m - n)*(a + b*x^n)^p
, x], x] /; FreeQ[{a, b, c, p}, x] && IGtQ[n, 0] && GtQ[m, n - 1] && NeQ[m + n*p + 1, 0] && IntBinomialQ[a, b,
c, n, m, p, x]

Rubi steps \begin{align*} \text {integral}& = \frac {1}{3} \text {Subst}\left (\int \frac {x^6}{\sqrt {1+x^2}} \, dx,x,x^3\right ) \\ & = \frac {1}{18} x^{15} \sqrt {1+x^6}-\frac {5}{18} \text {Subst}\left (\int \frac {x^4}{\sqrt {1+x^2}} \, dx,x,x^3\right ) \\ & = -\frac {5}{72} x^9 \sqrt {1+x^6}+\frac {1}{18} x^{15} \sqrt {1+x^6}+\frac {5}{24} \text {Subst}\left (\int \frac {x^2}{\sqrt {1+x^2}} \, dx,x,x^3\right ) \\ & = \frac {5}{48} x^3 \sqrt {1+x^6}-\frac {5}{72} x^9 \sqrt {1+x^6}+\frac {1}{18} x^{15} \sqrt {1+x^6}-\frac {5}{48} \text {Subst}\left (\int \frac {1}{\sqrt {1+x^2}} \, dx,x,x^3\right ) \\ & = \frac {5}{48} x^3 \sqrt {1+x^6}-\frac {5}{72} x^9 \sqrt {1+x^6}+\frac {1}{18} x^{15} \sqrt {1+x^6}-\frac {5 \text {arcsinh}\left (x^3\right )}{48} \\ \end{align*}

Mathematica [A] (verified)

Time = 0.14 (sec) , antiderivative size = 47, normalized size of antiderivative = 0.98 \[ \int \frac {x^{20}}{\sqrt {1+x^6}} \, dx=\frac {1}{144} x^3 \sqrt {1+x^6} \left (15-10 x^6+8 x^{12}\right )-\frac {5}{48} \log \left (x^3+\sqrt {1+x^6}\right ) \]

[In]

Integrate[x^20/Sqrt[1 + x^6],x]

[Out]

(x^3*Sqrt[1 + x^6]*(15 - 10*x^6 + 8*x^12))/144 - (5*Log[x^3 + Sqrt[1 + x^6]])/48

Maple [A] (verified)

Time = 1.06 (sec) , antiderivative size = 32, normalized size of antiderivative = 0.67

method	result	size
risch	\(\frac {x^{3} \left (8 x^{12}-10 x^{6}+15\right ) \sqrt {x^{6}+1}}{144}-\frac {5 \,\operatorname {arcsinh}\left (x^{3}\right )}{48}\)	\(32\)
pseudoelliptic	\(\frac {\sqrt {x^{6}+1}\, \left (8 x^{15}-10 x^{9}+15 x^{3}\right )}{144}-\frac {5 \,\operatorname {arcsinh}\left (x^{3}\right )}{48}\)	\(33\)
trager	\(\frac {x^{3} \left (8 x^{12}-10 x^{6}+15\right ) \sqrt {x^{6}+1}}{144}-\frac {5 \ln \left (x^{3}+\sqrt {x^{6}+1}\right )}{48}\)	\(40\)
meijerg	\(\frac {\frac {\sqrt {\pi }\, x^{3} \left (56 x^{12}-70 x^{6}+105\right ) \sqrt {x^{6}+1}}{168}-\frac {5 \sqrt {\pi }\, \operatorname {arcsinh}\left (x^{3}\right )}{8}}{6 \sqrt {\pi }}\)	\(43\)

[In]

int(x^20/(x^6+1)^(1/2),x,method=_RETURNVERBOSE)

[Out]

1/144*x^3*(8*x^12-10*x^6+15)*(x^6+1)^(1/2)-5/48*arcsinh(x^3)

Fricas [A] (verification not implemented)

none

Time = 0.25 (sec) , antiderivative size = 42, normalized size of antiderivative = 0.88 \[ \int \frac {x^{20}}{\sqrt {1+x^6}} \, dx=\frac {1}{144} \, {\left (8 \, x^{15} - 10 \, x^{9} + 15 \, x^{3}\right )} \sqrt {x^{6} + 1} + \frac {5}{48} \, \log \left (-x^{3} + \sqrt {x^{6} + 1}\right ) \]

[In]

integrate(x^20/(x^6+1)^(1/2),x, algorithm="fricas")

[Out]

1/144*(8*x^15 - 10*x^9 + 15*x^3)*sqrt(x^6 + 1) + 5/48*log(-x^3 + sqrt(x^6 + 1))

Sympy [A] (verification not implemented)

Time = 6.77 (sec) , antiderivative size = 65, normalized size of antiderivative = 1.35 \[ \int \frac {x^{20}}{\sqrt {1+x^6}} \, dx=\frac {x^{21}}{18 \sqrt {x^{6} + 1}} - \frac {x^{15}}{72 \sqrt {x^{6} + 1}} + \frac {5 x^{9}}{144 \sqrt {x^{6} + 1}} + \frac {5 x^{3}}{48 \sqrt {x^{6} + 1}} - \frac {5 \operatorname {asinh}{\left (x^{3} \right )}}{48} \]

[In]

integrate(x**20/(x**6+1)**(1/2),x)

[Out]

x**21/(18*sqrt(x**6 + 1)) - x**15/(72*sqrt(x**6 + 1)) + 5*x**9/(144*sqrt(x**6 + 1)) + 5*x**3/(48*sqrt(x**6 + 1
)) - 5*asinh(x**3)/48

Maxima [B] (verification not implemented)

Leaf count of result is larger than twice the leaf count of optimal. 109 vs. \(2 (40) = 80\).

Time = 0.17 (sec) , antiderivative size = 109, normalized size of antiderivative = 2.27 \[ \int \frac {x^{20}}{\sqrt {1+x^6}} \, dx=\frac {\frac {33 \, \sqrt {x^{6} + 1}}{x^{3}} - \frac {40 \, {\left (x^{6} + 1\right )}^{\frac {3}{2}}}{x^{9}} + \frac {15 \, {\left (x^{6} + 1\right )}^{\frac {5}{2}}}{x^{15}}}{144 \, {\left (\frac {3 \, {\left (x^{6} + 1\right )}}{x^{6}} - \frac {3 \, {\left (x^{6} + 1\right )}^{2}}{x^{12}} + \frac {{\left (x^{6} + 1\right )}^{3}}{x^{18}} - 1\right )}} - \frac {5}{96} \, \log \left (\frac {\sqrt {x^{6} + 1}}{x^{3}} + 1\right ) + \frac {5}{96} \, \log \left (\frac {\sqrt {x^{6} + 1}}{x^{3}} - 1\right ) \]

[In]

integrate(x^20/(x^6+1)^(1/2),x, algorithm="maxima")

[Out]

1/144*(33*sqrt(x^6 + 1)/x^3 - 40*(x^6 + 1)^(3/2)/x^9 + 15*(x^6 + 1)^(5/2)/x^15)/(3*(x^6 + 1)/x^6 - 3*(x^6 + 1)
^2/x^12 + (x^6 + 1)^3/x^18 - 1) - 5/96*log(sqrt(x^6 + 1)/x^3 + 1) + 5/96*log(sqrt(x^6 + 1)/x^3 - 1)

Giac [A] (verification not implemented)

none

Time = 0.29 (sec) , antiderivative size = 56, normalized size of antiderivative = 1.17 \[ \int \frac {x^{20}}{\sqrt {1+x^6}} \, dx=\frac {1}{144} \, {\left (2 \, {\left (4 \, x^{6} - 5\right )} x^{6} + 15\right )} \sqrt {x^{6} + 1} x^{3} - \frac {5 \, {\left (\log \left (\sqrt {\frac {1}{x^{6}} + 1} + 1\right ) - \log \left (\sqrt {\frac {1}{x^{6}} + 1} - 1\right )\right )}}{96 \, \mathrm {sgn}\left (x\right )} \]

[In]

integrate(x^20/(x^6+1)^(1/2),x, algorithm="giac")

[Out]

1/144*(2*(4*x^6 - 5)*x^6 + 15)*sqrt(x^6 + 1)*x^3 - 5/96*(log(sqrt(1/x^6 + 1) + 1) - log(sqrt(1/x^6 + 1) - 1))/
sgn(x)

Mupad [F(-1)]

Timed out. \[ \int \frac {x^{20}}{\sqrt {1+x^6}} \, dx=\int \frac {x^{20}}{\sqrt {x^6+1}} \,d x \]

[In]

int(x^20/(x^6 + 1)^(1/2),x)

[Out]

int(x^20/(x^6 + 1)^(1/2), x)