3.31.22 \(\int \frac {(-q+p x^4) \sqrt {q+p x^4}}{c x^4+b x^2 (q+p x^4)+a (q+p x^4)^2} \, dx\)

Optimal. Leaf size=428 \[ \frac {1}{4} \text {RootSum}\left [\text {$\#$1}^8 a-8 \text {$\#$1}^6 a \sqrt {p} \sqrt {q}+4 \text {$\#$1}^6 b+24 \text {$\#$1}^4 a p q-16 \text {$\#$1}^4 b \sqrt {p} \sqrt {q}+16 \text {$\#$1}^4 c-32 \text {$\#$1}^2 a p^{3/2} q^{3/2}+16 \text {$\#$1}^2 b p q+16 a p^2 q^2\& ,\frac {-\text {$\#$1}^6 \log \left (-\text {$\#$1} x+\sqrt {p x^4+q}-\sqrt {p} x^2-\sqrt {q}\right )+\text {$\#$1}^6 \log (x)+2 \text {$\#$1}^4 \sqrt {p} \sqrt {q} \log \left (-\text {$\#$1} x+\sqrt {p x^4+q}-\sqrt {p} x^2-\sqrt {q}\right )-2 \text {$\#$1}^4 \sqrt {p} \sqrt {q} \log (x)+4 \text {$\#$1}^2 p q \log \left (-\text {$\#$1} x+\sqrt {p x^4+q}-\sqrt {p} x^2-\sqrt {q}\right )-4 \text {$\#$1}^2 p q \log (x)-8 p^{3/2} q^{3/2} \log \left (-\text {$\#$1} x+\sqrt {p x^4+q}-\sqrt {p} x^2-\sqrt {q}\right )+8 p^{3/2} q^{3/2} \log (x)}{\text {$\#$1}^7 (-a)+6 \text {$\#$1}^5 a \sqrt {p} \sqrt {q}-3 \text {$\#$1}^5 b-12 \text {$\#$1}^3 a p q+8 \text {$\#$1}^3 b \sqrt {p} \sqrt {q}-8 \text {$\#$1}^3 c+8 \text {$\#$1} a p^{3/2} q^{3/2}-4 \text {$\#$1} b p q}\& \right ] \]

________________________________________________________________________________________

Rubi [F]  time = 1.70, antiderivative size = 0, normalized size of antiderivative = 0.00, number of steps used = 0, number of rules used = 0, integrand size = 0, \(\frac {\text {number of rules}}{\text {integrand size}}\) = 0.000, Rules used = {} \begin {gather*} \int \frac {\left (-q+p x^4\right ) \sqrt {q+p x^4}}{c x^4+b x^2 \left (q+p x^4\right )+a \left (q+p x^4\right )^2} \, dx \end {gather*}

Verification is not applicable to the result.

[In]

Int[((-q + p*x^4)*Sqrt[q + p*x^4])/(c*x^4 + b*x^2*(q + p*x^4) + a*(q + p*x^4)^2),x]

[Out]

q*Defer[Int][Sqrt[q + p*x^4]/(-(c*x^4) - b*x^2*(q + p*x^4) - a*(q + p*x^4)^2), x] + p*Defer[Int][(x^4*Sqrt[q +
 p*x^4])/(c*x^4 + b*x^2*(q + p*x^4) + a*(q + p*x^4)^2), x]

Rubi steps

\begin {align*} \int \frac {\left (-q+p x^4\right ) \sqrt {q+p x^4}}{c x^4+b x^2 \left (q+p x^4\right )+a \left (q+p x^4\right )^2} \, dx &=\int \left (\frac {q \sqrt {q+p x^4}}{-a q^2-b q x^2-c \left (1+\frac {2 a p q}{c}\right ) x^4-b p x^6-a p^2 x^8}+\frac {p x^4 \sqrt {q+p x^4}}{a q^2+b q x^2+c \left (1+\frac {2 a p q}{c}\right ) x^4+b p x^6+a p^2 x^8}\right ) \, dx\\ &=p \int \frac {x^4 \sqrt {q+p x^4}}{a q^2+b q x^2+c \left (1+\frac {2 a p q}{c}\right ) x^4+b p x^6+a p^2 x^8} \, dx+q \int \frac {\sqrt {q+p x^4}}{-a q^2-b q x^2-c \left (1+\frac {2 a p q}{c}\right ) x^4-b p x^6-a p^2 x^8} \, dx\\ &=p \int \frac {x^4 \sqrt {q+p x^4}}{c x^4+b x^2 \left (q+p x^4\right )+a \left (q+p x^4\right )^2} \, dx+q \int \frac {\sqrt {q+p x^4}}{-c x^4-b x^2 \left (q+p x^4\right )-a \left (q+p x^4\right )^2} \, dx\\ \end {align*}

________________________________________________________________________________________

Mathematica [C]  time = 6.71, size = 3958, normalized size = 9.25 \begin {gather*} \text {Result too large to show} \end {gather*}

Warning: Unable to verify antiderivative.

[In]

Integrate[((-q + p*x^4)*Sqrt[q + p*x^4])/(c*x^4 + b*x^2*(q + p*x^4) + a*(q + p*x^4)^2),x]

[Out]

((-I)*Sqrt[1 - (I*Sqrt[p]*x^2)/Sqrt[q]]*Sqrt[1 + (I*Sqrt[p]*x^2)/Sqrt[q]]*EllipticF[I*ArcSinh[Sqrt[(I*Sqrt[p])
/Sqrt[q]]*x], -1])/(a*Sqrt[(I*Sqrt[p])/Sqrt[q]]*Sqrt[q + p*x^4]) - (I*q^2*Sqrt[1 - (I*Sqrt[p]*x^2)/Sqrt[q]]*Sq
rt[1 + (I*Sqrt[p]*x^2)/Sqrt[q]]*EllipticPi[(I*Sqrt[q])/(Sqrt[p]*Root[a*q^2 + b*q*#1 + (c + 2*a*p*q)*#1^2 + b*p
*#1^3 + a*p^2*#1^4 & , 1]), I*ArcSinh[Sqrt[(I*Sqrt[p])/Sqrt[q]]*x], -1])/(a*p^2*Sqrt[(I*Sqrt[p])/Sqrt[q]]*Sqrt
[q + p*x^4]*Root[a*q^2 + b*q*#1 + (c + 2*a*p*q)*#1^2 + b*p*#1^3 + a*p^2*#1^4 & , 1]*(Root[a*q^2 + b*q*#1 + (c
+ 2*a*p*q)*#1^2 + b*p*#1^3 + a*p^2*#1^4 & , 1] - Root[a*q^2 + b*q*#1 + (c + 2*a*p*q)*#1^2 + b*p*#1^3 + a*p^2*#
1^4 & , 2])*(Root[a*q^2 + b*q*#1 + (c + 2*a*p*q)*#1^2 + b*p*#1^3 + a*p^2*#1^4 & , 1] - Root[a*q^2 + b*q*#1 + (
c + 2*a*p*q)*#1^2 + b*p*#1^3 + a*p^2*#1^4 & , 3])*(Root[a*q^2 + b*q*#1 + (c + 2*a*p*q)*#1^2 + b*p*#1^3 + a*p^2
*#1^4 & , 1] - Root[a*q^2 + b*q*#1 + (c + 2*a*p*q)*#1^2 + b*p*#1^3 + a*p^2*#1^4 & , 4])) + (I*Sqrt[1 - (I*Sqrt
[p]*x^2)/Sqrt[q]]*Sqrt[1 + (I*Sqrt[p]*x^2)/Sqrt[q]]*EllipticPi[(I*Sqrt[q])/(Sqrt[p]*Root[a*q^2 + b*q*#1 + (c +
 2*a*p*q)*#1^2 + b*p*#1^3 + a*p^2*#1^4 & , 1]), I*ArcSinh[Sqrt[(I*Sqrt[p])/Sqrt[q]]*x], -1]*Root[a*q^2 + b*q*#
1 + (c + 2*a*p*q)*#1^2 + b*p*#1^3 + a*p^2*#1^4 & , 1]^3)/(a*Sqrt[(I*Sqrt[p])/Sqrt[q]]*Sqrt[q + p*x^4]*(Root[a*
q^2 + b*q*#1 + (c + 2*a*p*q)*#1^2 + b*p*#1^3 + a*p^2*#1^4 & , 1] - Root[a*q^2 + b*q*#1 + (c + 2*a*p*q)*#1^2 +
b*p*#1^3 + a*p^2*#1^4 & , 2])*(Root[a*q^2 + b*q*#1 + (c + 2*a*p*q)*#1^2 + b*p*#1^3 + a*p^2*#1^4 & , 1] - Root[
a*q^2 + b*q*#1 + (c + 2*a*p*q)*#1^2 + b*p*#1^3 + a*p^2*#1^4 & , 3])*(Root[a*q^2 + b*q*#1 + (c + 2*a*p*q)*#1^2
+ b*p*#1^3 + a*p^2*#1^4 & , 1] - Root[a*q^2 + b*q*#1 + (c + 2*a*p*q)*#1^2 + b*p*#1^3 + a*p^2*#1^4 & , 4])) - (
I*q^2*Sqrt[1 - (I*Sqrt[p]*x^2)/Sqrt[q]]*Sqrt[1 + (I*Sqrt[p]*x^2)/Sqrt[q]]*EllipticPi[(I*Sqrt[q])/(Sqrt[p]*Root
[a*q^2 + b*q*#1 + (c + 2*a*p*q)*#1^2 + b*p*#1^3 + a*p^2*#1^4 & , 2]), I*ArcSinh[Sqrt[(I*Sqrt[p])/Sqrt[q]]*x],
-1])/(a*p^2*Sqrt[(I*Sqrt[p])/Sqrt[q]]*Sqrt[q + p*x^4]*Root[a*q^2 + b*q*#1 + (c + 2*a*p*q)*#1^2 + b*p*#1^3 + a*
p^2*#1^4 & , 2]*(-Root[a*q^2 + b*q*#1 + (c + 2*a*p*q)*#1^2 + b*p*#1^3 + a*p^2*#1^4 & , 1] + Root[a*q^2 + b*q*#
1 + (c + 2*a*p*q)*#1^2 + b*p*#1^3 + a*p^2*#1^4 & , 2])*(Root[a*q^2 + b*q*#1 + (c + 2*a*p*q)*#1^2 + b*p*#1^3 +
a*p^2*#1^4 & , 2] - Root[a*q^2 + b*q*#1 + (c + 2*a*p*q)*#1^2 + b*p*#1^3 + a*p^2*#1^4 & , 3])*(Root[a*q^2 + b*q
*#1 + (c + 2*a*p*q)*#1^2 + b*p*#1^3 + a*p^2*#1^4 & , 2] - Root[a*q^2 + b*q*#1 + (c + 2*a*p*q)*#1^2 + b*p*#1^3
+ a*p^2*#1^4 & , 4])) + (I*Sqrt[1 - (I*Sqrt[p]*x^2)/Sqrt[q]]*Sqrt[1 + (I*Sqrt[p]*x^2)/Sqrt[q]]*EllipticPi[(I*S
qrt[q])/(Sqrt[p]*Root[a*q^2 + b*q*#1 + (c + 2*a*p*q)*#1^2 + b*p*#1^3 + a*p^2*#1^4 & , 2]), I*ArcSinh[Sqrt[(I*S
qrt[p])/Sqrt[q]]*x], -1]*Root[a*q^2 + b*q*#1 + (c + 2*a*p*q)*#1^2 + b*p*#1^3 + a*p^2*#1^4 & , 2]^3)/(a*Sqrt[(I
*Sqrt[p])/Sqrt[q]]*Sqrt[q + p*x^4]*(-Root[a*q^2 + b*q*#1 + (c + 2*a*p*q)*#1^2 + b*p*#1^3 + a*p^2*#1^4 & , 1] +
 Root[a*q^2 + b*q*#1 + (c + 2*a*p*q)*#1^2 + b*p*#1^3 + a*p^2*#1^4 & , 2])*(Root[a*q^2 + b*q*#1 + (c + 2*a*p*q)
*#1^2 + b*p*#1^3 + a*p^2*#1^4 & , 2] - Root[a*q^2 + b*q*#1 + (c + 2*a*p*q)*#1^2 + b*p*#1^3 + a*p^2*#1^4 & , 3]
)*(Root[a*q^2 + b*q*#1 + (c + 2*a*p*q)*#1^2 + b*p*#1^3 + a*p^2*#1^4 & , 2] - Root[a*q^2 + b*q*#1 + (c + 2*a*p*
q)*#1^2 + b*p*#1^3 + a*p^2*#1^4 & , 4])) - (I*q^2*Sqrt[1 - (I*Sqrt[p]*x^2)/Sqrt[q]]*Sqrt[1 + (I*Sqrt[p]*x^2)/S
qrt[q]]*EllipticPi[(I*Sqrt[q])/(Sqrt[p]*Root[a*q^2 + b*q*#1 + (c + 2*a*p*q)*#1^2 + b*p*#1^3 + a*p^2*#1^4 & , 3
]), I*ArcSinh[Sqrt[(I*Sqrt[p])/Sqrt[q]]*x], -1])/(a*p^2*Sqrt[(I*Sqrt[p])/Sqrt[q]]*Sqrt[q + p*x^4]*Root[a*q^2 +
 b*q*#1 + (c + 2*a*p*q)*#1^2 + b*p*#1^3 + a*p^2*#1^4 & , 3]*(-Root[a*q^2 + b*q*#1 + (c + 2*a*p*q)*#1^2 + b*p*#
1^3 + a*p^2*#1^4 & , 1] + Root[a*q^2 + b*q*#1 + (c + 2*a*p*q)*#1^2 + b*p*#1^3 + a*p^2*#1^4 & , 3])*(-Root[a*q^
2 + b*q*#1 + (c + 2*a*p*q)*#1^2 + b*p*#1^3 + a*p^2*#1^4 & , 2] + Root[a*q^2 + b*q*#1 + (c + 2*a*p*q)*#1^2 + b*
p*#1^3 + a*p^2*#1^4 & , 3])*(Root[a*q^2 + b*q*#1 + (c + 2*a*p*q)*#1^2 + b*p*#1^3 + a*p^2*#1^4 & , 3] - Root[a*
q^2 + b*q*#1 + (c + 2*a*p*q)*#1^2 + b*p*#1^3 + a*p^2*#1^4 & , 4])) + (I*Sqrt[1 - (I*Sqrt[p]*x^2)/Sqrt[q]]*Sqrt
[1 + (I*Sqrt[p]*x^2)/Sqrt[q]]*EllipticPi[(I*Sqrt[q])/(Sqrt[p]*Root[a*q^2 + b*q*#1 + (c + 2*a*p*q)*#1^2 + b*p*#
1^3 + a*p^2*#1^4 & , 3]), I*ArcSinh[Sqrt[(I*Sqrt[p])/Sqrt[q]]*x], -1]*Root[a*q^2 + b*q*#1 + (c + 2*a*p*q)*#1^2
 + b*p*#1^3 + a*p^2*#1^4 & , 3]^3)/(a*Sqrt[(I*Sqrt[p])/Sqrt[q]]*Sqrt[q + p*x^4]*(-Root[a*q^2 + b*q*#1 + (c + 2
*a*p*q)*#1^2 + b*p*#1^3 + a*p^2*#1^4 & , 1] + Root[a*q^2 + b*q*#1 + (c + 2*a*p*q)*#1^2 + b*p*#1^3 + a*p^2*#1^4
 & , 3])*(-Root[a*q^2 + b*q*#1 + (c + 2*a*p*q)*#1^2 + b*p*#1^3 + a*p^2*#1^4 & , 2] + Root[a*q^2 + b*q*#1 + (c
+ 2*a*p*q)*#1^2 + b*p*#1^3 + a*p^2*#1^4 & , 3])*(Root[a*q^2 + b*q*#1 + (c + 2*a*p*q)*#1^2 + b*p*#1^3 + a*p^2*#
1^4 & , 3] - Root[a*q^2 + b*q*#1 + (c + 2*a*p*q)*#1^2 + b*p*#1^3 + a*p^2*#1^4 & , 4])) - (I*q^2*Sqrt[1 - (I*Sq
rt[p]*x^2)/Sqrt[q]]*Sqrt[1 + (I*Sqrt[p]*x^2)/Sqrt[q]]*EllipticPi[(I*Sqrt[q])/(Sqrt[p]*Root[a*q^2 + b*q*#1 + (c
 + 2*a*p*q)*#1^2 + b*p*#1^3 + a*p^2*#1^4 & , 4]), I*ArcSinh[Sqrt[(I*Sqrt[p])/Sqrt[q]]*x], -1])/(a*p^2*Sqrt[(I*
Sqrt[p])/Sqrt[q]]*Sqrt[q + p*x^4]*Root[a*q^2 + b*q*#1 + (c + 2*a*p*q)*#1^2 + b*p*#1^3 + a*p^2*#1^4 & , 4]*(-Ro
ot[a*q^2 + b*q*#1 + (c + 2*a*p*q)*#1^2 + b*p*#1^3 + a*p^2*#1^4 & , 1] + Root[a*q^2 + b*q*#1 + (c + 2*a*p*q)*#1
^2 + b*p*#1^3 + a*p^2*#1^4 & , 4])*(-Root[a*q^2 + b*q*#1 + (c + 2*a*p*q)*#1^2 + b*p*#1^3 + a*p^2*#1^4 & , 2] +
 Root[a*q^2 + b*q*#1 + (c + 2*a*p*q)*#1^2 + b*p*#1^3 + a*p^2*#1^4 & , 4])*(-Root[a*q^2 + b*q*#1 + (c + 2*a*p*q
)*#1^2 + b*p*#1^3 + a*p^2*#1^4 & , 3] + Root[a*q^2 + b*q*#1 + (c + 2*a*p*q)*#1^2 + b*p*#1^3 + a*p^2*#1^4 & , 4
])) + (I*Sqrt[1 - (I*Sqrt[p]*x^2)/Sqrt[q]]*Sqrt[1 + (I*Sqrt[p]*x^2)/Sqrt[q]]*EllipticPi[(I*Sqrt[q])/(Sqrt[p]*R
oot[a*q^2 + b*q*#1 + (c + 2*a*p*q)*#1^2 + b*p*#1^3 + a*p^2*#1^4 & , 4]), I*ArcSinh[Sqrt[(I*Sqrt[p])/Sqrt[q]]*x
], -1]*Root[a*q^2 + b*q*#1 + (c + 2*a*p*q)*#1^2 + b*p*#1^3 + a*p^2*#1^4 & , 4]^3)/(a*Sqrt[(I*Sqrt[p])/Sqrt[q]]
*Sqrt[q + p*x^4]*(-Root[a*q^2 + b*q*#1 + (c + 2*a*p*q)*#1^2 + b*p*#1^3 + a*p^2*#1^4 & , 1] + Root[a*q^2 + b*q*
#1 + (c + 2*a*p*q)*#1^2 + b*p*#1^3 + a*p^2*#1^4 & , 4])*(-Root[a*q^2 + b*q*#1 + (c + 2*a*p*q)*#1^2 + b*p*#1^3
+ a*p^2*#1^4 & , 2] + Root[a*q^2 + b*q*#1 + (c + 2*a*p*q)*#1^2 + b*p*#1^3 + a*p^2*#1^4 & , 4])*(-Root[a*q^2 +
b*q*#1 + (c + 2*a*p*q)*#1^2 + b*p*#1^3 + a*p^2*#1^4 & , 3] + Root[a*q^2 + b*q*#1 + (c + 2*a*p*q)*#1^2 + b*p*#1
^3 + a*p^2*#1^4 & , 4]))

________________________________________________________________________________________

IntegrateAlgebraic [A]  time = 3.58, size = 189, normalized size = 0.44 \begin {gather*} -\frac {\left (-b+\sqrt {b^2-4 a c}\right ) \tan ^{-1}\left (\frac {\sqrt {b-\sqrt {b^2-4 a c}} x}{\sqrt {2} \sqrt {a} \sqrt {q+p x^4}}\right )}{\sqrt {2} \sqrt {a} \sqrt {b^2-4 a c} \sqrt {b-\sqrt {b^2-4 a c}}}-\frac {\sqrt {b+\sqrt {b^2-4 a c}} \tan ^{-1}\left (\frac {\sqrt {b+\sqrt {b^2-4 a c}} x}{\sqrt {2} \sqrt {a} \sqrt {q+p x^4}}\right )}{\sqrt {2} \sqrt {a} \sqrt {b^2-4 a c}} \end {gather*}

Antiderivative was successfully verified.

[In]

IntegrateAlgebraic[((-q + p*x^4)*Sqrt[q + p*x^4])/(c*x^4 + b*x^2*(q + p*x^4) + a*(q + p*x^4)^2),x]

[Out]

-(((-b + Sqrt[b^2 - 4*a*c])*ArcTan[(Sqrt[b - Sqrt[b^2 - 4*a*c]]*x)/(Sqrt[2]*Sqrt[a]*Sqrt[q + p*x^4])])/(Sqrt[2
]*Sqrt[a]*Sqrt[b^2 - 4*a*c]*Sqrt[b - Sqrt[b^2 - 4*a*c]])) - (Sqrt[b + Sqrt[b^2 - 4*a*c]]*ArcTan[(Sqrt[b + Sqrt
[b^2 - 4*a*c]]*x)/(Sqrt[2]*Sqrt[a]*Sqrt[q + p*x^4])])/(Sqrt[2]*Sqrt[a]*Sqrt[b^2 - 4*a*c])

________________________________________________________________________________________

fricas [B]  time = 178.32, size = 2083, normalized size = 4.87

result too large to display

Verification of antiderivative is not currently implemented for this CAS.

[In]

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

[Out]

-1/4*sqrt(1/2)*sqrt(-(b + (a*b^2 - 4*a^2*c)/sqrt(a^2*b^2 - 4*a^3*c))/(a*b^2 - 4*a^2*c))*log(-(sqrt(1/2)*((a*b^
2 - 4*a^2*c)*p^2*x^8 - (b^3 - 4*a*b*c)*p*x^6 - (b^2*c - 4*a*c^2 - 2*(a*b^2 - 4*a^2*c)*p*q)*x^4 - (b^3 - 4*a*b*
c)*q*x^2 + (a*b^2 - 4*a^2*c)*q^2 + ((a^2*b^3 - 4*a^3*b*c)*p^2*x^8 - (a*b^4 - 8*a^2*b^2*c + 16*a^3*c^2)*p*x^6 -
 (a*b^3*c - 4*a^2*b*c^2 - 2*(a^2*b^3 - 4*a^3*b*c)*p*q)*x^4 - (a*b^4 - 8*a^2*b^2*c + 16*a^3*c^2)*q*x^2 + (a^2*b
^3 - 4*a^3*b*c)*q^2)/sqrt(a^2*b^2 - 4*a^3*c))*sqrt(-(b + (a*b^2 - 4*a^2*c)/sqrt(a^2*b^2 - 4*a^3*c))/(a*b^2 - 4
*a^2*c)) + 2*((b^2 - 2*a*c)*p*x^5 + b*c*x^3 + (b^2 - 2*a*c)*q*x + ((a*b^3 - 4*a^2*b*c)*p*x^5 + (a*b^2*c - 4*a^
2*c^2)*x^3 + (a*b^3 - 4*a^2*b*c)*q*x)/sqrt(a^2*b^2 - 4*a^3*c))*sqrt(p*x^4 + q))/(a*p^2*x^8 + b*p*x^6 + (2*a*p*
q + c)*x^4 + b*q*x^2 + a*q^2)) + 1/4*sqrt(1/2)*sqrt(-(b + (a*b^2 - 4*a^2*c)/sqrt(a^2*b^2 - 4*a^3*c))/(a*b^2 -
4*a^2*c))*log((sqrt(1/2)*((a*b^2 - 4*a^2*c)*p^2*x^8 - (b^3 - 4*a*b*c)*p*x^6 - (b^2*c - 4*a*c^2 - 2*(a*b^2 - 4*
a^2*c)*p*q)*x^4 - (b^3 - 4*a*b*c)*q*x^2 + (a*b^2 - 4*a^2*c)*q^2 + ((a^2*b^3 - 4*a^3*b*c)*p^2*x^8 - (a*b^4 - 8*
a^2*b^2*c + 16*a^3*c^2)*p*x^6 - (a*b^3*c - 4*a^2*b*c^2 - 2*(a^2*b^3 - 4*a^3*b*c)*p*q)*x^4 - (a*b^4 - 8*a^2*b^2
*c + 16*a^3*c^2)*q*x^2 + (a^2*b^3 - 4*a^3*b*c)*q^2)/sqrt(a^2*b^2 - 4*a^3*c))*sqrt(-(b + (a*b^2 - 4*a^2*c)/sqrt
(a^2*b^2 - 4*a^3*c))/(a*b^2 - 4*a^2*c)) - 2*((b^2 - 2*a*c)*p*x^5 + b*c*x^3 + (b^2 - 2*a*c)*q*x + ((a*b^3 - 4*a
^2*b*c)*p*x^5 + (a*b^2*c - 4*a^2*c^2)*x^3 + (a*b^3 - 4*a^2*b*c)*q*x)/sqrt(a^2*b^2 - 4*a^3*c))*sqrt(p*x^4 + q))
/(a*p^2*x^8 + b*p*x^6 + (2*a*p*q + c)*x^4 + b*q*x^2 + a*q^2)) - 1/4*sqrt(1/2)*sqrt(-(b - (a*b^2 - 4*a^2*c)/sqr
t(a^2*b^2 - 4*a^3*c))/(a*b^2 - 4*a^2*c))*log(-(sqrt(1/2)*((a*b^2 - 4*a^2*c)*p^2*x^8 - (b^3 - 4*a*b*c)*p*x^6 -
(b^2*c - 4*a*c^2 - 2*(a*b^2 - 4*a^2*c)*p*q)*x^4 - (b^3 - 4*a*b*c)*q*x^2 + (a*b^2 - 4*a^2*c)*q^2 - ((a^2*b^3 -
4*a^3*b*c)*p^2*x^8 - (a*b^4 - 8*a^2*b^2*c + 16*a^3*c^2)*p*x^6 - (a*b^3*c - 4*a^2*b*c^2 - 2*(a^2*b^3 - 4*a^3*b*
c)*p*q)*x^4 - (a*b^4 - 8*a^2*b^2*c + 16*a^3*c^2)*q*x^2 + (a^2*b^3 - 4*a^3*b*c)*q^2)/sqrt(a^2*b^2 - 4*a^3*c))*s
qrt(-(b - (a*b^2 - 4*a^2*c)/sqrt(a^2*b^2 - 4*a^3*c))/(a*b^2 - 4*a^2*c)) + 2*((b^2 - 2*a*c)*p*x^5 + b*c*x^3 + (
b^2 - 2*a*c)*q*x - ((a*b^3 - 4*a^2*b*c)*p*x^5 + (a*b^2*c - 4*a^2*c^2)*x^3 + (a*b^3 - 4*a^2*b*c)*q*x)/sqrt(a^2*
b^2 - 4*a^3*c))*sqrt(p*x^4 + q))/(a*p^2*x^8 + b*p*x^6 + (2*a*p*q + c)*x^4 + b*q*x^2 + a*q^2)) + 1/4*sqrt(1/2)*
sqrt(-(b - (a*b^2 - 4*a^2*c)/sqrt(a^2*b^2 - 4*a^3*c))/(a*b^2 - 4*a^2*c))*log((sqrt(1/2)*((a*b^2 - 4*a^2*c)*p^2
*x^8 - (b^3 - 4*a*b*c)*p*x^6 - (b^2*c - 4*a*c^2 - 2*(a*b^2 - 4*a^2*c)*p*q)*x^4 - (b^3 - 4*a*b*c)*q*x^2 + (a*b^
2 - 4*a^2*c)*q^2 - ((a^2*b^3 - 4*a^3*b*c)*p^2*x^8 - (a*b^4 - 8*a^2*b^2*c + 16*a^3*c^2)*p*x^6 - (a*b^3*c - 4*a^
2*b*c^2 - 2*(a^2*b^3 - 4*a^3*b*c)*p*q)*x^4 - (a*b^4 - 8*a^2*b^2*c + 16*a^3*c^2)*q*x^2 + (a^2*b^3 - 4*a^3*b*c)*
q^2)/sqrt(a^2*b^2 - 4*a^3*c))*sqrt(-(b - (a*b^2 - 4*a^2*c)/sqrt(a^2*b^2 - 4*a^3*c))/(a*b^2 - 4*a^2*c)) - 2*((b
^2 - 2*a*c)*p*x^5 + b*c*x^3 + (b^2 - 2*a*c)*q*x - ((a*b^3 - 4*a^2*b*c)*p*x^5 + (a*b^2*c - 4*a^2*c^2)*x^3 + (a*
b^3 - 4*a^2*b*c)*q*x)/sqrt(a^2*b^2 - 4*a^3*c))*sqrt(p*x^4 + q))/(a*p^2*x^8 + b*p*x^6 + (2*a*p*q + c)*x^4 + b*q
*x^2 + a*q^2))

________________________________________________________________________________________

giac [F(-1)]  time = 0.00, size = 0, normalized size = 0.00 \begin {gather*} \text {Timed out} \end {gather*}

Verification of antiderivative is not currently implemented for this CAS.

[In]

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

[Out]

Timed out

________________________________________________________________________________________

maple [B]  time = 0.20, size = 242, normalized size = 0.57

method result size
default \(\frac {\left (\frac {\arctan \left (\frac {\sqrt {p \,x^{4}+q}\, \sqrt {2}\, a}{x \sqrt {\left (b +\sqrt {-4 a c +b^{2}}\right ) a}}\right )}{\sqrt {\left (b +\sqrt {-4 a c +b^{2}}\right ) a}}+\frac {\arctan \left (\frac {\sqrt {p \,x^{4}+q}\, \sqrt {2}\, a}{x \sqrt {\left (b +\sqrt {-4 a c +b^{2}}\right ) a}}\right ) b}{\sqrt {-4 a c +b^{2}}\, \sqrt {\left (b +\sqrt {-4 a c +b^{2}}\right ) a}}-\frac {\arctanh \left (\frac {\sqrt {p \,x^{4}+q}\, \sqrt {2}\, a}{x \sqrt {\left (-b +\sqrt {-4 a c +b^{2}}\right ) a}}\right )}{\sqrt {\left (-b +\sqrt {-4 a c +b^{2}}\right ) a}}+\frac {\arctanh \left (\frac {\sqrt {p \,x^{4}+q}\, \sqrt {2}\, a}{x \sqrt {\left (-b +\sqrt {-4 a c +b^{2}}\right ) a}}\right ) b}{\sqrt {-4 a c +b^{2}}\, \sqrt {\left (-b +\sqrt {-4 a c +b^{2}}\right ) a}}\right ) \sqrt {2}}{2}\) \(242\)
elliptic \(\frac {\left (\frac {\arctan \left (\frac {\sqrt {p \,x^{4}+q}\, \sqrt {2}\, a}{x \sqrt {\left (b +\sqrt {-4 a c +b^{2}}\right ) a}}\right )}{\sqrt {\left (b +\sqrt {-4 a c +b^{2}}\right ) a}}+\frac {\arctan \left (\frac {\sqrt {p \,x^{4}+q}\, \sqrt {2}\, a}{x \sqrt {\left (b +\sqrt {-4 a c +b^{2}}\right ) a}}\right ) b}{\sqrt {-4 a c +b^{2}}\, \sqrt {\left (b +\sqrt {-4 a c +b^{2}}\right ) a}}-\frac {\arctanh \left (\frac {\sqrt {p \,x^{4}+q}\, \sqrt {2}\, a}{x \sqrt {\left (-b +\sqrt {-4 a c +b^{2}}\right ) a}}\right )}{\sqrt {\left (-b +\sqrt {-4 a c +b^{2}}\right ) a}}+\frac {\arctanh \left (\frac {\sqrt {p \,x^{4}+q}\, \sqrt {2}\, a}{x \sqrt {\left (-b +\sqrt {-4 a c +b^{2}}\right ) a}}\right ) b}{\sqrt {-4 a c +b^{2}}\, \sqrt {\left (-b +\sqrt {-4 a c +b^{2}}\right ) a}}\right ) \sqrt {2}}{2}\) \(242\)

Verification of antiderivative is not currently implemented for this CAS.

[In]

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

[Out]

1/2*(1/((b+(-4*a*c+b^2)^(1/2))*a)^(1/2)*arctan((p*x^4+q)^(1/2)*2^(1/2)/x*a/((b+(-4*a*c+b^2)^(1/2))*a)^(1/2))+1
/(-4*a*c+b^2)^(1/2)/((b+(-4*a*c+b^2)^(1/2))*a)^(1/2)*arctan((p*x^4+q)^(1/2)*2^(1/2)/x*a/((b+(-4*a*c+b^2)^(1/2)
)*a)^(1/2))*b-1/((-b+(-4*a*c+b^2)^(1/2))*a)^(1/2)*arctanh((p*x^4+q)^(1/2)*2^(1/2)/x*a/((-b+(-4*a*c+b^2)^(1/2))
*a)^(1/2))+1/(-4*a*c+b^2)^(1/2)/((-b+(-4*a*c+b^2)^(1/2))*a)^(1/2)*arctanh((p*x^4+q)^(1/2)*2^(1/2)/x*a/((-b+(-4
*a*c+b^2)^(1/2))*a)^(1/2))*b)*2^(1/2)

________________________________________________________________________________________

maxima [F]  time = 0.00, size = 0, normalized size = 0.00 \begin {gather*} \int \frac {\sqrt {p x^{4} + q} {\left (p x^{4} - q\right )}}{c x^{4} + {\left (p x^{4} + q\right )} b x^{2} + {\left (p x^{4} + q\right )}^{2} a}\,{d x} \end {gather*}

Verification of antiderivative is not currently implemented for this CAS.

[In]

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

[Out]

integrate(sqrt(p*x^4 + q)*(p*x^4 - q)/(c*x^4 + (p*x^4 + q)*b*x^2 + (p*x^4 + q)^2*a), x)

________________________________________________________________________________________

mupad [F]  time = 0.00, size = -1, normalized size = -0.00 \begin {gather*} \int -\frac {\sqrt {p\,x^4+q}\,\left (q-p\,x^4\right )}{a\,{\left (p\,x^4+q\right )}^2+c\,x^4+b\,x^2\,\left (p\,x^4+q\right )} \,d x \end {gather*}

Verification of antiderivative is not currently implemented for this CAS.

[In]

int(-((q + p*x^4)^(1/2)*(q - p*x^4))/(a*(q + p*x^4)^2 + c*x^4 + b*x^2*(q + p*x^4)),x)

[Out]

int(-((q + p*x^4)^(1/2)*(q - p*x^4))/(a*(q + p*x^4)^2 + c*x^4 + b*x^2*(q + p*x^4)), x)

________________________________________________________________________________________

sympy [F(-1)]  time = 0.00, size = 0, normalized size = 0.00 \begin {gather*} \text {Timed out} \end {gather*}

Verification of antiderivative is not currently implemented for this CAS.

[In]

integrate((p*x**4-q)*(p*x**4+q)**(1/2)/(c*x**4+b*x**2*(p*x**4+q)+a*(p*x**4+q)**2),x)

[Out]

Timed out

________________________________________________________________________________________