\(\int \frac {x^4}{(-1+x^4)^2 \sqrt [4]{x^2+x^4}} \, dx\) [1417]

   Optimal result
   Rubi [C] (verified)
   Mathematica [A] (verified)
   Maple [A] (verified)
   Fricas [C] (verification not implemented)
   Sympy [F]
   Maxima [F]
   Giac [A] (verification not implemented)
   Mupad [F(-1)]

Optimal result

Integrand size = 22, antiderivative size = 101 \[ \int \frac {x^4}{\left (-1+x^4\right )^2 \sqrt [4]{x^2+x^4}} \, dx=\frac {\left (-5-2 x^2-13 x^4\right ) \left (x^2+x^4\right )^{3/4}}{80 x \left (-1+x^2\right ) \left (1+x^2\right )^2}-\frac {\arctan \left (\frac {\sqrt [4]{2} x}{\sqrt [4]{x^2+x^4}}\right )}{32 \sqrt [4]{2}}-\frac {\text {arctanh}\left (\frac {\sqrt [4]{2} x}{\sqrt [4]{x^2+x^4}}\right )}{32 \sqrt [4]{2}} \]

[Out]

1/80*(-13*x^4-2*x^2-5)*(x^4+x^2)^(3/4)/x/(x^2-1)/(x^2+1)^2-1/64*arctan(2^(1/4)*x/(x^4+x^2)^(1/4))*2^(3/4)-1/64
*arctanh(2^(1/4)*x/(x^4+x^2)^(1/4))*2^(3/4)

Rubi [C] (verified)

Result contains higher order function than in optimal. Order 5 vs. order 3 in optimal.

Time = 4.62 (sec) , antiderivative size = 109, normalized size of antiderivative = 1.08, number of steps used = 4, number of rules used = 4, \(\frac {\text {number of rules}}{\text {integrand size}}\) = 0.182, Rules used = {2081, 1268, 477, 524} \[ \int \frac {x^4}{\left (-1+x^4\right )^2 \sqrt [4]{x^2+x^4}} \, dx=\frac {128 x^5 \operatorname {Gamma}\left (\frac {13}{4}\right ) \left (17 \left (-4 x^4-9 x^2+13\right ) \operatorname {Hypergeometric2F1}\left (1,2,\frac {17}{4},-\frac {2 x^2}{1-x^2}\right )-64 x^2 \left (x^2+1\right ) \operatorname {Hypergeometric2F1}\left (2,3,\frac {21}{4},-\frac {2 x^2}{1-x^2}\right )\right )}{89505 \left (1-x^2\right )^3 \left (x^2+1\right ) \sqrt [4]{x^4+x^2} \operatorname {Gamma}\left (\frac {1}{4}\right )} \]

[In]

Int[x^4/((-1 + x^4)^2*(x^2 + x^4)^(1/4)),x]

[Out]

(128*x^5*Gamma[13/4]*(17*(13 - 9*x^2 - 4*x^4)*Hypergeometric2F1[1, 2, 17/4, (-2*x^2)/(1 - x^2)] - 64*x^2*(1 +
x^2)*Hypergeometric2F1[2, 3, 21/4, (-2*x^2)/(1 - x^2)]))/(89505*(1 - x^2)^3*(1 + x^2)*(x^2 + x^4)^(1/4)*Gamma[
1/4])

Rule 477

Int[((e_.)*(x_))^(m_)*((a_) + (b_.)*(x_)^(n_))^(p_)*((c_) + (d_.)*(x_)^(n_))^(q_), x_Symbol] :> With[{k = Deno
minator[m]}, Dist[k/e, Subst[Int[x^(k*(m + 1) - 1)*(a + b*(x^(k*n)/e^n))^p*(c + d*(x^(k*n)/e^n))^q, x], x, (e*
x)^(1/k)], x]] /; FreeQ[{a, b, c, d, e, p, q}, x] && NeQ[b*c - a*d, 0] && IGtQ[n, 0] && FractionQ[m] && Intege
rQ[p]

Rule 524

Int[((e_.)*(x_))^(m_.)*((a_) + (b_.)*(x_)^(n_))^(p_)*((c_) + (d_.)*(x_)^(n_))^(q_), x_Symbol] :> Simp[a^p*c^q*
((e*x)^(m + 1)/(e*(m + 1)))*AppellF1[(m + 1)/n, -p, -q, 1 + (m + 1)/n, (-b)*(x^n/a), (-d)*(x^n/c)], x] /; Free
Q[{a, b, c, d, e, m, n, p, q}, x] && NeQ[b*c - a*d, 0] && NeQ[m, -1] && NeQ[m, n - 1] && (IntegerQ[p] || GtQ[a
, 0]) && (IntegerQ[q] || GtQ[c, 0])

Rule 1268

Int[((f_.)*(x_))^(m_.)*((d_) + (e_.)*(x_)^2)^(q_.)*((a_) + (c_.)*(x_)^4)^(p_.), x_Symbol] :> Int[(f*x)^m*(d +
e*x^2)^(q + p)*(a/d + (c/e)*x^2)^p, x] /; FreeQ[{a, c, d, e, f, q, m, q}, x] && EqQ[c*d^2 + a*e^2, 0] && Integ
erQ[p]

Rule 2081

Int[(u_.)*(P_)^(p_.), x_Symbol] :> With[{m = MinimumMonomialExponent[P, x]}, Dist[P^FracPart[p]/(x^(m*FracPart
[p])*Distrib[1/x^m, P]^FracPart[p]), Int[u*x^(m*p)*Distrib[1/x^m, P]^p, x], x]] /; FreeQ[p, x] &&  !IntegerQ[p
] && SumQ[P] && EveryQ[BinomialQ[#1, x] & , P] &&  !PolyQ[P, x, 2]

Rubi steps \begin{align*} \text {integral}& = \frac {\left (\sqrt {x} \sqrt [4]{1+x^2}\right ) \int \frac {x^{7/2}}{\sqrt [4]{1+x^2} \left (-1+x^4\right )^2} \, dx}{\sqrt [4]{x^2+x^4}} \\ & = \frac {\left (\sqrt {x} \sqrt [4]{1+x^2}\right ) \int \frac {x^{7/2}}{\left (-1+x^2\right )^2 \left (1+x^2\right )^{9/4}} \, dx}{\sqrt [4]{x^2+x^4}} \\ & = \frac {\left (2 \sqrt {x} \sqrt [4]{1+x^2}\right ) \text {Subst}\left (\int \frac {x^8}{\left (-1+x^4\right )^2 \left (1+x^4\right )^{9/4}} \, dx,x,\sqrt {x}\right )}{\sqrt [4]{x^2+x^4}} \\ & = \frac {128 x^5 \operatorname {Gamma}\left (\frac {13}{4}\right ) \left (17 \left (13-9 x^2-4 x^4\right ) \operatorname {Hypergeometric2F1}\left (1,2,\frac {17}{4},-\frac {2 x^2}{1-x^2}\right )-64 x^2 \left (1+x^2\right ) \operatorname {Hypergeometric2F1}\left (2,3,\frac {21}{4},-\frac {2 x^2}{1-x^2}\right )\right )}{89505 \left (1-x^2\right )^3 \left (1+x^2\right ) \sqrt [4]{x^2+x^4} \operatorname {Gamma}\left (\frac {1}{4}\right )} \\ \end{align*}

Mathematica [A] (verified)

Time = 0.48 (sec) , antiderivative size = 131, normalized size of antiderivative = 1.30 \[ \int \frac {x^4}{\left (-1+x^4\right )^2 \sqrt [4]{x^2+x^4}} \, dx=-\frac {\sqrt {x} \left (4 \sqrt {x} \left (5+2 x^2+13 x^4\right )+5\ 2^{3/4} \sqrt [4]{1+x^2} \left (-1+x^4\right ) \arctan \left (\frac {\sqrt [4]{2} \sqrt {x}}{\sqrt [4]{1+x^2}}\right )+5\ 2^{3/4} \sqrt [4]{1+x^2} \left (-1+x^4\right ) \text {arctanh}\left (\frac {\sqrt [4]{2} \sqrt {x}}{\sqrt [4]{1+x^2}}\right )\right )}{320 \left (-1+x^4\right ) \sqrt [4]{x^2+x^4}} \]

[In]

Integrate[x^4/((-1 + x^4)^2*(x^2 + x^4)^(1/4)),x]

[Out]

-1/320*(Sqrt[x]*(4*Sqrt[x]*(5 + 2*x^2 + 13*x^4) + 5*2^(3/4)*(1 + x^2)^(1/4)*(-1 + x^4)*ArcTan[(2^(1/4)*Sqrt[x]
)/(1 + x^2)^(1/4)] + 5*2^(3/4)*(1 + x^2)^(1/4)*(-1 + x^4)*ArcTanh[(2^(1/4)*Sqrt[x])/(1 + x^2)^(1/4)]))/((-1 +
x^4)*(x^2 + x^4)^(1/4))

Maple [A] (verified)

Time = 4.98 (sec) , antiderivative size = 122, normalized size of antiderivative = 1.21

method result size
pseudoelliptic \(\frac {5 \left (x^{2} \left (x^{2}+1\right )\right )^{\frac {1}{4}} \left (x^{4}-1\right ) \left (2 \arctan \left (\frac {\left (x^{2} \left (x^{2}+1\right )\right )^{\frac {1}{4}} 2^{\frac {3}{4}}}{2 x}\right )-\ln \left (\frac {2^{\frac {1}{4}} x +\left (x^{2} \left (x^{2}+1\right )\right )^{\frac {1}{4}}}{-2^{\frac {1}{4}} x +\left (x^{2} \left (x^{2}+1\right )\right )^{\frac {1}{4}}}\right )\right ) 2^{\frac {3}{4}}-104 x^{5}-16 x^{3}-40 x}{\left (x^{2} \left (x^{2}+1\right )\right )^{\frac {1}{4}} \left (640 x^{4}-640\right )}\) \(122\)
trager \(-\frac {\left (13 x^{4}+2 x^{2}+5\right ) \left (x^{4}+x^{2}\right )^{\frac {3}{4}}}{80 \left (x^{2}-1\right ) \left (x^{2}+1\right )^{2} x}-\frac {\operatorname {RootOf}\left (\textit {\_Z}^{2}+\operatorname {RootOf}\left (\textit {\_Z}^{4}-8\right )^{2}\right ) \ln \left (\frac {-\sqrt {x^{4}+x^{2}}\, \operatorname {RootOf}\left (\textit {\_Z}^{2}+\operatorname {RootOf}\left (\textit {\_Z}^{4}-8\right )^{2}\right ) \operatorname {RootOf}\left (\textit {\_Z}^{4}-8\right )^{2} x -2 \left (x^{4}+x^{2}\right )^{\frac {1}{4}} \operatorname {RootOf}\left (\textit {\_Z}^{4}-8\right )^{2} x^{2}+3 \operatorname {RootOf}\left (\textit {\_Z}^{2}+\operatorname {RootOf}\left (\textit {\_Z}^{4}-8\right )^{2}\right ) x^{3}+4 \left (x^{4}+x^{2}\right )^{\frac {3}{4}}+\operatorname {RootOf}\left (\textit {\_Z}^{2}+\operatorname {RootOf}\left (\textit {\_Z}^{4}-8\right )^{2}\right ) x}{\left (1+x \right ) \left (-1+x \right ) x}\right )}{128}-\frac {\operatorname {RootOf}\left (\textit {\_Z}^{4}-8\right ) \ln \left (\frac {\sqrt {x^{4}+x^{2}}\, \operatorname {RootOf}\left (\textit {\_Z}^{4}-8\right )^{3} x +2 \left (x^{4}+x^{2}\right )^{\frac {1}{4}} \operatorname {RootOf}\left (\textit {\_Z}^{4}-8\right )^{2} x^{2}+3 \operatorname {RootOf}\left (\textit {\_Z}^{4}-8\right ) x^{3}+4 \left (x^{4}+x^{2}\right )^{\frac {3}{4}}+\operatorname {RootOf}\left (\textit {\_Z}^{4}-8\right ) x}{\left (1+x \right ) \left (-1+x \right ) x}\right )}{128}\) \(267\)
risch \(-\frac {x \left (13 x^{4}+2 x^{2}+5\right )}{80 \left (x^{2}-1\right ) \left (x^{2} \left (x^{2}+1\right )\right )^{\frac {1}{4}} \left (x^{2}+1\right )}-\frac {\operatorname {RootOf}\left (\textit {\_Z}^{4}-8\right ) \ln \left (\frac {\sqrt {x^{4}+x^{2}}\, \operatorname {RootOf}\left (\textit {\_Z}^{4}-8\right )^{3} x +2 \left (x^{4}+x^{2}\right )^{\frac {1}{4}} \operatorname {RootOf}\left (\textit {\_Z}^{4}-8\right )^{2} x^{2}+3 \operatorname {RootOf}\left (\textit {\_Z}^{4}-8\right ) x^{3}+4 \left (x^{4}+x^{2}\right )^{\frac {3}{4}}+\operatorname {RootOf}\left (\textit {\_Z}^{4}-8\right ) x}{\left (1+x \right ) \left (-1+x \right ) x}\right )}{128}+\frac {\operatorname {RootOf}\left (\textit {\_Z}^{2}+\operatorname {RootOf}\left (\textit {\_Z}^{4}-8\right )^{2}\right ) \ln \left (\frac {\sqrt {x^{4}+x^{2}}\, \operatorname {RootOf}\left (\textit {\_Z}^{2}+\operatorname {RootOf}\left (\textit {\_Z}^{4}-8\right )^{2}\right ) \operatorname {RootOf}\left (\textit {\_Z}^{4}-8\right )^{2} x -2 \left (x^{4}+x^{2}\right )^{\frac {1}{4}} \operatorname {RootOf}\left (\textit {\_Z}^{4}-8\right )^{2} x^{2}-3 \operatorname {RootOf}\left (\textit {\_Z}^{2}+\operatorname {RootOf}\left (\textit {\_Z}^{4}-8\right )^{2}\right ) x^{3}+4 \left (x^{4}+x^{2}\right )^{\frac {3}{4}}-\operatorname {RootOf}\left (\textit {\_Z}^{2}+\operatorname {RootOf}\left (\textit {\_Z}^{4}-8\right )^{2}\right ) x}{\left (1+x \right ) \left (-1+x \right ) x}\right )}{128}\) \(267\)

[In]

int(x^4/(x^4-1)^2/(x^4+x^2)^(1/4),x,method=_RETURNVERBOSE)

[Out]

10/(x^2*(x^2+1))^(1/4)*(1/2*(x^2*(x^2+1))^(1/4)*(x^4-1)*(2*arctan(1/2*(x^2*(x^2+1))^(1/4)/x*2^(3/4))-ln((2^(1/
4)*x+(x^2*(x^2+1))^(1/4))/(-2^(1/4)*x+(x^2*(x^2+1))^(1/4))))*2^(3/4)-52/5*x^5-8/5*x^3-4*x)/(640*x^4-640)

Fricas [C] (verification not implemented)

Result contains complex when optimal does not.

Time = 1.01 (sec) , antiderivative size = 404, normalized size of antiderivative = 4.00 \[ \int \frac {x^4}{\left (-1+x^4\right )^2 \sqrt [4]{x^2+x^4}} \, dx=-\frac {5 \cdot 2^{\frac {3}{4}} {\left (x^{7} + x^{5} - x^{3} - x\right )} \log \left (\frac {4 \, \sqrt {2} {\left (x^{4} + x^{2}\right )}^{\frac {1}{4}} x^{2} + 2^{\frac {3}{4}} {\left (3 \, x^{3} + x\right )} + 4 \cdot 2^{\frac {1}{4}} \sqrt {x^{4} + x^{2}} x + 4 \, {\left (x^{4} + x^{2}\right )}^{\frac {3}{4}}}{x^{3} - x}\right ) - 5 \cdot 2^{\frac {3}{4}} {\left (x^{7} + x^{5} - x^{3} - x\right )} \log \left (\frac {4 \, \sqrt {2} {\left (x^{4} + x^{2}\right )}^{\frac {1}{4}} x^{2} - 2^{\frac {3}{4}} {\left (3 \, x^{3} + x\right )} - 4 \cdot 2^{\frac {1}{4}} \sqrt {x^{4} + x^{2}} x + 4 \, {\left (x^{4} + x^{2}\right )}^{\frac {3}{4}}}{x^{3} - x}\right ) + 5 \cdot 2^{\frac {3}{4}} {\left (i \, x^{7} + i \, x^{5} - i \, x^{3} - i \, x\right )} \log \left (-\frac {4 \, \sqrt {2} {\left (x^{4} + x^{2}\right )}^{\frac {1}{4}} x^{2} - 2^{\frac {3}{4}} {\left (3 i \, x^{3} + i \, x\right )} + 4 i \cdot 2^{\frac {1}{4}} \sqrt {x^{4} + x^{2}} x - 4 \, {\left (x^{4} + x^{2}\right )}^{\frac {3}{4}}}{x^{3} - x}\right ) + 5 \cdot 2^{\frac {3}{4}} {\left (-i \, x^{7} - i \, x^{5} + i \, x^{3} + i \, x\right )} \log \left (-\frac {4 \, \sqrt {2} {\left (x^{4} + x^{2}\right )}^{\frac {1}{4}} x^{2} - 2^{\frac {3}{4}} {\left (-3 i \, x^{3} - i \, x\right )} - 4 i \cdot 2^{\frac {1}{4}} \sqrt {x^{4} + x^{2}} x - 4 \, {\left (x^{4} + x^{2}\right )}^{\frac {3}{4}}}{x^{3} - x}\right ) + 16 \, {\left (13 \, x^{4} + 2 \, x^{2} + 5\right )} {\left (x^{4} + x^{2}\right )}^{\frac {3}{4}}}{1280 \, {\left (x^{7} + x^{5} - x^{3} - x\right )}} \]

[In]

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

[Out]

-1/1280*(5*2^(3/4)*(x^7 + x^5 - x^3 - x)*log((4*sqrt(2)*(x^4 + x^2)^(1/4)*x^2 + 2^(3/4)*(3*x^3 + x) + 4*2^(1/4
)*sqrt(x^4 + x^2)*x + 4*(x^4 + x^2)^(3/4))/(x^3 - x)) - 5*2^(3/4)*(x^7 + x^5 - x^3 - x)*log((4*sqrt(2)*(x^4 +
x^2)^(1/4)*x^2 - 2^(3/4)*(3*x^3 + x) - 4*2^(1/4)*sqrt(x^4 + x^2)*x + 4*(x^4 + x^2)^(3/4))/(x^3 - x)) + 5*2^(3/
4)*(I*x^7 + I*x^5 - I*x^3 - I*x)*log(-(4*sqrt(2)*(x^4 + x^2)^(1/4)*x^2 - 2^(3/4)*(3*I*x^3 + I*x) + 4*I*2^(1/4)
*sqrt(x^4 + x^2)*x - 4*(x^4 + x^2)^(3/4))/(x^3 - x)) + 5*2^(3/4)*(-I*x^7 - I*x^5 + I*x^3 + I*x)*log(-(4*sqrt(2
)*(x^4 + x^2)^(1/4)*x^2 - 2^(3/4)*(-3*I*x^3 - I*x) - 4*I*2^(1/4)*sqrt(x^4 + x^2)*x - 4*(x^4 + x^2)^(3/4))/(x^3
 - x)) + 16*(13*x^4 + 2*x^2 + 5)*(x^4 + x^2)^(3/4))/(x^7 + x^5 - x^3 - x)

Sympy [F]

\[ \int \frac {x^4}{\left (-1+x^4\right )^2 \sqrt [4]{x^2+x^4}} \, dx=\int \frac {x^{4}}{\sqrt [4]{x^{2} \left (x^{2} + 1\right )} \left (x - 1\right )^{2} \left (x + 1\right )^{2} \left (x^{2} + 1\right )^{2}}\, dx \]

[In]

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

[Out]

Integral(x**4/((x**2*(x**2 + 1))**(1/4)*(x - 1)**2*(x + 1)**2*(x**2 + 1)**2), x)

Maxima [F]

\[ \int \frac {x^4}{\left (-1+x^4\right )^2 \sqrt [4]{x^2+x^4}} \, dx=\int { \frac {x^{4}}{{\left (x^{4} + x^{2}\right )}^{\frac {1}{4}} {\left (x^{4} - 1\right )}^{2}} \,d x } \]

[In]

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

[Out]

2/21*(4*x^5 + x^3 - 3*x)*x^(7/2)/((x^8 - 2*x^4 + 1)*(x^2 + 1)^(1/4)) + integrate(16/21*(4*x^4 + x^2 - 3)*x^(7/
2)/((x^12 - 3*x^8 + 3*x^4 - 1)*(x^2 + 1)^(1/4)), x)

Giac [A] (verification not implemented)

none

Time = 0.30 (sec) , antiderivative size = 81, normalized size of antiderivative = 0.80 \[ \int \frac {x^4}{\left (-1+x^4\right )^2 \sqrt [4]{x^2+x^4}} \, dx=\frac {1}{64} \cdot 2^{\frac {3}{4}} \arctan \left (\frac {1}{2} \cdot 2^{\frac {3}{4}} {\left (\frac {1}{x^{2}} + 1\right )}^{\frac {1}{4}}\right ) - \frac {1}{128} \cdot 2^{\frac {3}{4}} \log \left (2^{\frac {1}{4}} + {\left (\frac {1}{x^{2}} + 1\right )}^{\frac {1}{4}}\right ) + \frac {1}{128} \cdot 2^{\frac {3}{4}} \log \left ({\left | -2^{\frac {1}{4}} + {\left (\frac {1}{x^{2}} + 1\right )}^{\frac {1}{4}} \right |}\right ) + \frac {{\left (\frac {1}{x^{2}} + 1\right )}^{\frac {3}{4}}}{16 \, {\left (\frac {1}{x^{2}} - 1\right )}} - \frac {1}{10 \, {\left (\frac {1}{x^{2}} + 1\right )}^{\frac {5}{4}}} \]

[In]

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

[Out]

1/64*2^(3/4)*arctan(1/2*2^(3/4)*(1/x^2 + 1)^(1/4)) - 1/128*2^(3/4)*log(2^(1/4) + (1/x^2 + 1)^(1/4)) + 1/128*2^
(3/4)*log(abs(-2^(1/4) + (1/x^2 + 1)^(1/4))) + 1/16*(1/x^2 + 1)^(3/4)/(1/x^2 - 1) - 1/10/(1/x^2 + 1)^(5/4)

Mupad [F(-1)]

Timed out. \[ \int \frac {x^4}{\left (-1+x^4\right )^2 \sqrt [4]{x^2+x^4}} \, dx=\int \frac {x^4}{{\left (x^4+x^2\right )}^{1/4}\,{\left (x^4-1\right )}^2} \,d x \]

[In]

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

[Out]

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