3.7.58 \(\int \frac {\sqrt {x+x^4} (-b+a x^6)}{x^6} \, dx\) [658]

Optimal. Leaf size=52 \[ \frac {\sqrt {x+x^4} \left (2 b+2 b x^3+3 a x^6\right )}{9 x^5}+\frac {1}{3} a \tanh ^{-1}\left (\frac {x^2}{\sqrt {x+x^4}}\right ) \]

[Out]

1/9*(x^4+x)^(1/2)*(3*a*x^6+2*b*x^3+2*b)/x^5+1/3*a*arctanh(x^2/(x^4+x)^(1/2))

________________________________________________________________________________________

Rubi [C] Result contains higher order function than in optimal. Order 4 vs. order 3 in optimal.
time = 0.18, antiderivative size = 163, normalized size of antiderivative = 3.13, number of steps used = 8, number of rules used = 7, integrand size = 22, \(\frac {\text {number of rules}}{\text {integrand size}}\) = 0.318, Rules used = {2073, 2077, 2039, 2045, 2036, 335, 231} \begin {gather*} \frac {3^{3/4} a x (x+1) \sqrt {\frac {x^2-x+1}{\left (\left (1+\sqrt {3}\right ) x+1\right )^2}} F\left (\text {ArcCos}\left (\frac {\left (1-\sqrt {3}\right ) x+1}{\left (1+\sqrt {3}\right ) x+1}\right )|\frac {1}{4} \left (2+\sqrt {3}\right )\right )}{5 \sqrt {\frac {x (x+1)}{\left (\left (1+\sqrt {3}\right ) x+1\right )^2}} \sqrt {x^4+x}}+\frac {a \left (x^4+x\right )^{3/2}}{3 x^3}-\frac {2 a \sqrt {x^4+x}}{5 x^3}+\frac {2 b \left (x^4+x\right )^{3/2}}{9 x^6} \end {gather*}

Warning: Unable to verify antiderivative.

[In]

Int[(Sqrt[x + x^4]*(-b + a*x^6))/x^6,x]

[Out]

(-2*a*Sqrt[x + x^4])/(5*x^3) + (2*b*(x + x^4)^(3/2))/(9*x^6) + (a*(x + x^4)^(3/2))/(3*x^3) + (3^(3/4)*a*x*(1 +
x)*Sqrt[(1 - x + x^2)/(1 + (1 + Sqrt[3])*x)^2]*EllipticF[ArcCos[(1 + (1 - Sqrt[3])*x)/(1 + (1 + Sqrt[3])*x)],
(2 + Sqrt[3])/4])/(5*Sqrt[(x*(1 + x))/(1 + (1 + Sqrt[3])*x)^2]*Sqrt[x + x^4])

Rule 231

Int[1/Sqrt[(a_) + (b_.)*(x_)^6], x_Symbol] :> With[{r = Numer[Rt[b/a, 3]], s = Denom[Rt[b/a, 3]]}, Simp[x*(s +
r*x^2)*(Sqrt[(s^2 - r*s*x^2 + r^2*x^4)/(s + (1 + Sqrt[3])*r*x^2)^2]/(2*3^(1/4)*s*Sqrt[a + b*x^6]*Sqrt[r*x^2*(
(s + r*x^2)/(s + (1 + Sqrt[3])*r*x^2)^2)]))*EllipticF[ArcCos[(s + (1 - Sqrt[3])*r*x^2)/(s + (1 + Sqrt[3])*r*x^
2)], (2 + Sqrt[3])/4], x]] /; FreeQ[{a, b}, x]

Rule 335

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

Rule 2036

Int[((a_.)*(x_)^(j_.) + (b_.)*(x_)^(n_.))^(p_), x_Symbol] :> Dist[(a*x^j + b*x^n)^FracPart[p]/(x^(j*FracPart[p
])*(a + b*x^(n - j))^FracPart[p]), Int[x^(j*p)*(a + b*x^(n - j))^p, x], x] /; FreeQ[{a, b, j, n, p}, x] &&  !I
ntegerQ[p] && NeQ[n, j] && PosQ[n - j]

Rule 2039

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

Rule 2045

Int[((c_.)*(x_))^(m_)*((a_.)*(x_)^(j_.) + (b_.)*(x_)^(n_.))^(p_), x_Symbol] :> Simp[(c*x)^(m + 1)*((a*x^j + b*
x^n)^p/(c*(m + j*p + 1))), x] - Dist[b*p*((n - j)/(c^n*(m + j*p + 1))), Int[(c*x)^(m + n)*(a*x^j + b*x^n)^(p -
1), x], x] /; FreeQ[{a, b, c}, x] &&  !IntegerQ[p] && LtQ[0, j, n] && (IntegersQ[j, n] || GtQ[c, 0]) && GtQ[p
, 0] && LtQ[m + j*p + 1, 0]

Rule 2073

Int[(Pq_)*((c_.)*(x_))^(m_.)*((a_.)*(x_)^(j_.) + (b_.)*(x_)^(n_))^(p_), x_Symbol] :> With[{q = Expon[Pq, x]},
With[{Pqq = Coeff[Pq, x, q]}, Int[(c*x)^m*ExpandToSum[Pq - Pqq*x^q - a*Pqq*(m + q - n + 1)*(x^(q - n)/(b*(m +
q + n*p + 1))), x]*(a*x^j + b*x^n)^p, x] + Simp[Pqq*(c*x)^(m + q - n + 1)*((a*x^j + b*x^n)^(p + 1)/(b*c^(q - n
+ 1)*(m + q + n*p + 1))), x]] /; GtQ[q, n - 1] && NeQ[m + q + n*p + 1, 0] && (IntegerQ[2*p] || IntegerQ[p + (
q + 1)/(2*n)])] /; FreeQ[{a, b, c, m, p}, x] && PolyQ[Pq, x] &&  !IntegerQ[p] && IGtQ[j, 0] && IGtQ[n, 0] && L
tQ[j, n]

Rule 2077

Int[(Pq_)*((c_.)*(x_))^(m_.)*((a_.)*(x_)^(j_.) + (b_.)*(x_)^(n_))^(p_), x_Symbol] :> Int[ExpandIntegrand[(c*x)
^m*Pq*(a*x^j + b*x^n)^p, x], x] /; FreeQ[{a, b, c, j, m, n, p}, x] && (PolyQ[Pq, x] || PolyQ[Pq, x^n]) &&  !In
tegerQ[p] && NeQ[n, j]

Rubi steps

\begin {align*} \int \frac {\sqrt {x+x^4} \left (-b+a x^6\right )}{x^6} \, dx &=\frac {a \left (x+x^4\right )^{3/2}}{3 x^3}+\int \frac {\left (-b+a x^2\right ) \sqrt {x+x^4}}{x^6} \, dx\\ &=\frac {a \left (x+x^4\right )^{3/2}}{3 x^3}+\int \left (-\frac {b \sqrt {x+x^4}}{x^6}+\frac {a \sqrt {x+x^4}}{x^4}\right ) \, dx\\ &=\frac {a \left (x+x^4\right )^{3/2}}{3 x^3}+a \int \frac {\sqrt {x+x^4}}{x^4} \, dx-b \int \frac {\sqrt {x+x^4}}{x^6} \, dx\\ &=-\frac {2 a \sqrt {x+x^4}}{5 x^3}+\frac {2 b \left (x+x^4\right )^{3/2}}{9 x^6}+\frac {a \left (x+x^4\right )^{3/2}}{3 x^3}+\frac {1}{5} (3 a) \int \frac {1}{\sqrt {x+x^4}} \, dx\\ &=-\frac {2 a \sqrt {x+x^4}}{5 x^3}+\frac {2 b \left (x+x^4\right )^{3/2}}{9 x^6}+\frac {a \left (x+x^4\right )^{3/2}}{3 x^3}+\frac {\left (3 a \sqrt {x} \sqrt {1+x^3}\right ) \int \frac {1}{\sqrt {x} \sqrt {1+x^3}} \, dx}{5 \sqrt {x+x^4}}\\ &=-\frac {2 a \sqrt {x+x^4}}{5 x^3}+\frac {2 b \left (x+x^4\right )^{3/2}}{9 x^6}+\frac {a \left (x+x^4\right )^{3/2}}{3 x^3}+\frac {\left (6 a \sqrt {x} \sqrt {1+x^3}\right ) \text {Subst}\left (\int \frac {1}{\sqrt {1+x^6}} \, dx,x,\sqrt {x}\right )}{5 \sqrt {x+x^4}}\\ &=-\frac {2 a \sqrt {x+x^4}}{5 x^3}+\frac {2 b \left (x+x^4\right )^{3/2}}{9 x^6}+\frac {a \left (x+x^4\right )^{3/2}}{3 x^3}+\frac {3^{3/4} a x (1+x) \sqrt {\frac {1-x+x^2}{\left (1+\left (1+\sqrt {3}\right ) x\right )^2}} F\left (\cos ^{-1}\left (\frac {1+\left (1-\sqrt {3}\right ) x}{1+\left (1+\sqrt {3}\right ) x}\right )|\frac {1}{4} \left (2+\sqrt {3}\right )\right )}{5 \sqrt {\frac {x (1+x)}{\left (1+\left (1+\sqrt {3}\right ) x\right )^2}} \sqrt {x+x^4}}\\ \end {align*}

________________________________________________________________________________________

Mathematica [A]
time = 0.18, size = 77, normalized size = 1.48 \begin {gather*} \frac {\sqrt {x+x^4} \left (2 b+2 b x^3+3 a x^6\right )}{9 x^5}+\frac {a \sqrt {x+x^4} \tanh ^{-1}\left (\frac {x^{3/2}}{\sqrt {1+x^3}}\right )}{3 \sqrt {x} \sqrt {1+x^3}} \end {gather*}

Antiderivative was successfully verified.

[In]

Integrate[(Sqrt[x + x^4]*(-b + a*x^6))/x^6,x]

[Out]

(Sqrt[x + x^4]*(2*b + 2*b*x^3 + 3*a*x^6))/(9*x^5) + (a*Sqrt[x + x^4]*ArcTanh[x^(3/2)/Sqrt[1 + x^3]])/(3*Sqrt[x
]*Sqrt[1 + x^3])

________________________________________________________________________________________

Maple [C] Result contains higher order function than in optimal. Order 4 vs. order 3.
time = 1.54, size = 332, normalized size = 6.38 Too large to display

Verification of antiderivative is not currently implemented for this CAS.

[In]

int((x^4+x)^(1/2)*(a*x^6-b)/x^6,x,method=_RETURNVERBOSE)

[Out]

a*(1/3*x*(x^4+x)^(1/2)-(-1/2-1/2*I*3^(1/2))*((3/2+1/2*I*3^(1/2))*x/(1/2+1/2*I*3^(1/2))/(1+x))^(1/2)*(1+x)^2*(-
(x-1/2+1/2*I*3^(1/2))/(1/2-1/2*I*3^(1/2))/(1+x))^(1/2)*(-(x-1/2-1/2*I*3^(1/2))/(1/2+1/2*I*3^(1/2))/(1+x))^(1/2
)/(3/2+1/2*I*3^(1/2))/(x*(1+x)*(x-1/2+1/2*I*3^(1/2))*(x-1/2-1/2*I*3^(1/2)))^(1/2)*(-EllipticF(((3/2+1/2*I*3^(1
/2))*x/(1/2+1/2*I*3^(1/2))/(1+x))^(1/2),((-3/2+1/2*I*3^(1/2))*(-1/2-1/2*I*3^(1/2))/(-1/2+1/2*I*3^(1/2))/(-3/2-
1/2*I*3^(1/2)))^(1/2))+EllipticPi(((3/2+1/2*I*3^(1/2))*x/(1/2+1/2*I*3^(1/2))/(1+x))^(1/2),(1/2+1/2*I*3^(1/2))/
(3/2+1/2*I*3^(1/2)),((-3/2+1/2*I*3^(1/2))*(-1/2-1/2*I*3^(1/2))/(-1/2+1/2*I*3^(1/2))/(-3/2-1/2*I*3^(1/2)))^(1/2
))))-b*(-2/9*(x^4+x)^(1/2)/x^5-2/9*(x^4+x)^(1/2)/x^2)

________________________________________________________________________________________

Maxima [F]
time = 0.00, size = 0, normalized size = 0.00 \begin {gather*} \text {Failed to integrate} \end {gather*}

Verification of antiderivative is not currently implemented for this CAS.

[In]

integrate((x^4+x)^(1/2)*(a*x^6-b)/x^6,x, algorithm="maxima")

[Out]

integrate((a*x^6 - b)*sqrt(x^4 + x)/x^6, x)

________________________________________________________________________________________

Fricas [A]
time = 0.39, size = 55, normalized size = 1.06 \begin {gather*} \frac {3 \, a x^{5} \log \left (-2 \, x^{3} - 2 \, \sqrt {x^{4} + x} x - 1\right ) + 2 \, {\left (3 \, a x^{6} + 2 \, b x^{3} + 2 \, b\right )} \sqrt {x^{4} + x}}{18 \, x^{5}} \end {gather*}

Verification of antiderivative is not currently implemented for this CAS.

[In]

integrate((x^4+x)^(1/2)*(a*x^6-b)/x^6,x, algorithm="fricas")

[Out]

1/18*(3*a*x^5*log(-2*x^3 - 2*sqrt(x^4 + x)*x - 1) + 2*(3*a*x^6 + 2*b*x^3 + 2*b)*sqrt(x^4 + x))/x^5

________________________________________________________________________________________

Sympy [F]
time = 0.00, size = 0, normalized size = 0.00 \begin {gather*} \int \frac {\sqrt {x \left (x + 1\right ) \left (x^{2} - x + 1\right )} \left (a x^{6} - b\right )}{x^{6}}\, dx \end {gather*}

Verification of antiderivative is not currently implemented for this CAS.

[In]

integrate((x**4+x)**(1/2)*(a*x**6-b)/x**6,x)

[Out]

Integral(sqrt(x*(x + 1)*(x**2 - x + 1))*(a*x**6 - b)/x**6, x)

________________________________________________________________________________________

Giac [A]
time = 0.43, size = 49, normalized size = 0.94 \begin {gather*} \frac {1}{3} \, \sqrt {x^{4} + x} a x + \frac {2}{9} \, b {\left (\frac {1}{x^{3}} + 1\right )}^{\frac {3}{2}} + \frac {1}{6} \, a \log \left (\sqrt {\frac {1}{x^{3}} + 1} + 1\right ) - \frac {1}{6} \, a \log \left ({\left | \sqrt {\frac {1}{x^{3}} + 1} - 1 \right |}\right ) \end {gather*}

Verification of antiderivative is not currently implemented for this CAS.

[In]

integrate((x^4+x)^(1/2)*(a*x^6-b)/x^6,x, algorithm="giac")

[Out]

1/3*sqrt(x^4 + x)*a*x + 2/9*b*(1/x^3 + 1)^(3/2) + 1/6*a*log(sqrt(1/x^3 + 1) + 1) - 1/6*a*log(abs(sqrt(1/x^3 +
1) - 1))

________________________________________________________________________________________

Mupad [B]
time = 0.99, size = 47, normalized size = 0.90 \begin {gather*} \frac {2\,b\,\left (x^3+1\right )\,\sqrt {x^4+x}}{9\,x^5}+\frac {2\,a\,x\,\sqrt {x^4+x}\,{{}}_2{\mathrm {F}}_1\left (-\frac {1}{2},\frac {1}{2};\ \frac {3}{2};\ -x^3\right )}{3\,\sqrt {x^3+1}} \end {gather*}

Verification of antiderivative is not currently implemented for this CAS.

[In]

int(-((b - a*x^6)*(x + x^4)^(1/2))/x^6,x)

[Out]

(2*b*(x^3 + 1)*(x + x^4)^(1/2))/(9*x^5) + (2*a*x*(x + x^4)^(1/2)*hypergeom([-1/2, 1/2], 3/2, -x^3))/(3*(x^3 +
1)^(1/2))

________________________________________________________________________________________