Integrand size = 20, antiderivative size = 196 \[ \int \frac {x^2}{a+(b+d) x^2+c x^4} \, dx=-\frac {\sqrt {b+d-\sqrt {b^2-4 a c+2 b d+d^2}} \arctan \left (\frac {\sqrt {2} \sqrt {c} x}{\sqrt {b+d-\sqrt {b^2-4 a c+2 b d+d^2}}}\right )}{\sqrt {2} \sqrt {c} \sqrt {b^2-4 a c+2 b d+d^2}}+\frac {\sqrt {b+d+\sqrt {b^2-4 a c+2 b d+d^2}} \arctan \left (\frac {\sqrt {2} \sqrt {c} x}{\sqrt {b+d+\sqrt {b^2-4 a c+2 b d+d^2}}}\right )}{\sqrt {2} \sqrt {c} \sqrt {b^2-4 a c+2 b d+d^2}} \] Output:
-1/2*(b+d-(-4*a*c+b^2+2*b*d+d^2)^(1/2))^(1/2)*arctan(2^(1/2)*c^(1/2)*x/(b+ d-(-4*a*c+b^2+2*b*d+d^2)^(1/2))^(1/2))*2^(1/2)/c^(1/2)/(-4*a*c+b^2+2*b*d+d ^2)^(1/2)+1/2*(b+d+(-4*a*c+b^2+2*b*d+d^2)^(1/2))^(1/2)*arctan(2^(1/2)*c^(1 /2)*x/(b+d+(-4*a*c+b^2+2*b*d+d^2)^(1/2))^(1/2))*2^(1/2)/c^(1/2)/(-4*a*c+b^ 2+2*b*d+d^2)^(1/2)
Time = 0.12 (sec) , antiderivative size = 221, normalized size of antiderivative = 1.13 \[ \int \frac {x^2}{a+(b+d) x^2+c x^4} \, dx=\frac {\left (-b-d+\sqrt {b^2-4 a c+2 b d+d^2}\right ) \arctan \left (\frac {\sqrt {2} \sqrt {c} x}{\sqrt {b+d-\sqrt {b^2-4 a c+2 b d+d^2}}}\right )}{\sqrt {2} \sqrt {c} \sqrt {b^2-4 a c+2 b d+d^2} \sqrt {b+d-\sqrt {b^2-4 a c+2 b d+d^2}}}+\frac {\sqrt {b+d+\sqrt {b^2-4 a c+2 b d+d^2}} \arctan \left (\frac {\sqrt {2} \sqrt {c} x}{\sqrt {b+d+\sqrt {b^2-4 a c+2 b d+d^2}}}\right )}{\sqrt {2} \sqrt {c} \sqrt {b^2-4 a c+2 b d+d^2}} \] Input:
Integrate[x^2/(a + (b + d)*x^2 + c*x^4),x]
Output:
((-b - d + Sqrt[b^2 - 4*a*c + 2*b*d + d^2])*ArcTan[(Sqrt[2]*Sqrt[c]*x)/Sqr t[b + d - Sqrt[b^2 - 4*a*c + 2*b*d + d^2]]])/(Sqrt[2]*Sqrt[c]*Sqrt[b^2 - 4 *a*c + 2*b*d + d^2]*Sqrt[b + d - Sqrt[b^2 - 4*a*c + 2*b*d + d^2]]) + (Sqrt [b + d + Sqrt[b^2 - 4*a*c + 2*b*d + d^2]]*ArcTan[(Sqrt[2]*Sqrt[c]*x)/Sqrt[ b + d + Sqrt[b^2 - 4*a*c + 2*b*d + d^2]]])/(Sqrt[2]*Sqrt[c]*Sqrt[b^2 - 4*a *c + 2*b*d + d^2])
Time = 0.73 (sec) , antiderivative size = 208, normalized size of antiderivative = 1.06, number of steps used = 2, number of rules used = 2, \(\frac {\text {number of rules}}{\text {integrand size}}\) = 0.100, Rules used = {1450, 218}
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 \frac {x^2}{a+x^2 (b+d)+c x^4} \, dx\) |
\(\Big \downarrow \) 1450 |
\(\displaystyle \frac {1}{2} \left (1-\frac {b+d}{\sqrt {-4 a c+b^2+2 b d+d^2}}\right ) \int \frac {1}{c x^2+\frac {1}{2} \left (b+d-\sqrt {b^2+2 d b+d^2-4 a c}\right )}dx+\frac {1}{2} \left (\frac {b+d}{\sqrt {-4 a c+b^2+2 b d+d^2}}+1\right ) \int \frac {1}{c x^2+\frac {1}{2} \left (b+d+\sqrt {b^2+2 d b+d^2-4 a c}\right )}dx\) |
\(\Big \downarrow \) 218 |
\(\displaystyle \frac {\left (1-\frac {b+d}{\sqrt {-4 a c+b^2+2 b d+d^2}}\right ) \arctan \left (\frac {\sqrt {2} \sqrt {c} x}{\sqrt {-\sqrt {-4 a c+b^2+2 b d+d^2}+b+d}}\right )}{\sqrt {2} \sqrt {c} \sqrt {-\sqrt {-4 a c+b^2+2 b d+d^2}+b+d}}+\frac {\left (\frac {b+d}{\sqrt {-4 a c+b^2+2 b d+d^2}}+1\right ) \arctan \left (\frac {\sqrt {2} \sqrt {c} x}{\sqrt {\sqrt {-4 a c+b^2+2 b d+d^2}+b+d}}\right )}{\sqrt {2} \sqrt {c} \sqrt {\sqrt {-4 a c+b^2+2 b d+d^2}+b+d}}\) |
Input:
Int[x^2/(a + (b + d)*x^2 + c*x^4),x]
Output:
((1 - (b + d)/Sqrt[b^2 - 4*a*c + 2*b*d + d^2])*ArcTan[(Sqrt[2]*Sqrt[c]*x)/ Sqrt[b + d - Sqrt[b^2 - 4*a*c + 2*b*d + d^2]]])/(Sqrt[2]*Sqrt[c]*Sqrt[b + d - Sqrt[b^2 - 4*a*c + 2*b*d + d^2]]) + ((1 + (b + d)/Sqrt[b^2 - 4*a*c + 2 *b*d + d^2])*ArcTan[(Sqrt[2]*Sqrt[c]*x)/Sqrt[b + d + Sqrt[b^2 - 4*a*c + 2* b*d + d^2]]])/(Sqrt[2]*Sqrt[c]*Sqrt[b + d + Sqrt[b^2 - 4*a*c + 2*b*d + d^2 ]])
Int[((a_) + (b_.)*(x_)^2)^(-1), x_Symbol] :> Simp[(Rt[a/b, 2]/a)*ArcTan[x/R t[a/b, 2]], x] /; FreeQ[{a, b}, x] && PosQ[a/b]
Int[((d_.)*(x_))^(m_)/((a_) + (b_.)*(x_)^2 + (c_.)*(x_)^4), x_Symbol] :> Wi th[{q = Rt[b^2 - 4*a*c, 2]}, Simp[(d^2/2)*(b/q + 1) Int[(d*x)^(m - 2)/(b/ 2 + q/2 + c*x^2), x], x] - Simp[(d^2/2)*(b/q - 1) Int[(d*x)^(m - 2)/(b/2 - q/2 + c*x^2), x], x]] /; FreeQ[{a, b, c, d}, x] && NeQ[b^2 - 4*a*c, 0] && GeQ[m, 2]
Result contains higher order function than in optimal. Order 9 vs. order 3.
Time = 0.16 (sec) , antiderivative size = 46, normalized size of antiderivative = 0.23
method | result | size |
risch | \(\frac {\left (\munderset {\textit {\_R} =\operatorname {RootOf}\left (c \,\textit {\_Z}^{4}+\left (b +d \right ) \textit {\_Z}^{2}+a \right )}{\sum }\frac {\textit {\_R}^{2} \ln \left (x -\textit {\_R} \right )}{2 \textit {\_R}^{3} c +\textit {\_R} b +\textit {\_R} d}\right )}{2}\) | \(46\) |
default | \(4 c \left (-\frac {\left (\sqrt {-4 a c +b^{2}+2 b d +d^{2}}-b -d \right ) \sqrt {2}\, \operatorname {arctanh}\left (\frac {c x \sqrt {2}}{\sqrt {\left (\sqrt {-4 a c +b^{2}+2 b d +d^{2}}-b -d \right ) c}}\right )}{8 c \sqrt {-4 a c +b^{2}+2 b d +d^{2}}\, \sqrt {\left (\sqrt {-4 a c +b^{2}+2 b d +d^{2}}-b -d \right ) c}}+\frac {\left (b +d +\sqrt {-4 a c +b^{2}+2 b d +d^{2}}\right ) \sqrt {2}\, \arctan \left (\frac {c x \sqrt {2}}{\sqrt {\left (b +d +\sqrt {-4 a c +b^{2}+2 b d +d^{2}}\right ) c}}\right )}{8 \sqrt {-4 a c +b^{2}+2 b d +d^{2}}\, c \sqrt {\left (b +d +\sqrt {-4 a c +b^{2}+2 b d +d^{2}}\right ) c}}\right )\) | \(217\) |
Input:
int(x^2/(a+(b+d)*x^2+c*x^4),x,method=_RETURNVERBOSE)
Output:
1/2*sum(_R^2/(2*_R^3*c+_R*b+_R*d)*ln(x-_R),_R=RootOf(c*_Z^4+(b+d)*_Z^2+a))
Leaf count of result is larger than twice the leaf count of optimal. 935 vs. \(2 (161) = 322\).
Time = 0.08 (sec) , antiderivative size = 935, normalized size of antiderivative = 4.77 \[ \int \frac {x^2}{a+(b+d) x^2+c x^4} \, dx =\text {Too large to display} \] Input:
integrate(x^2/(a+(b+d)*x^2+c*x^4),x, algorithm="fricas")
Output:
1/2*sqrt(1/2)*sqrt(-(b + d + (b^2*c - 4*a*c^2 + 2*b*c*d + c*d^2)/sqrt(b^2* c^2 - 4*a*c^3 + 2*b*c^2*d + c^2*d^2))/(b^2*c - 4*a*c^2 + 2*b*c*d + c*d^2)) *log(sqrt(1/2)*(b^2*c - 4*a*c^2 + 2*b*c*d + c*d^2)*sqrt(-(b + d + (b^2*c - 4*a*c^2 + 2*b*c*d + c*d^2)/sqrt(b^2*c^2 - 4*a*c^3 + 2*b*c^2*d + c^2*d^2)) /(b^2*c - 4*a*c^2 + 2*b*c*d + c*d^2))/sqrt(b^2*c^2 - 4*a*c^3 + 2*b*c^2*d + c^2*d^2) + x) - 1/2*sqrt(1/2)*sqrt(-(b + d + (b^2*c - 4*a*c^2 + 2*b*c*d + c*d^2)/sqrt(b^2*c^2 - 4*a*c^3 + 2*b*c^2*d + c^2*d^2))/(b^2*c - 4*a*c^2 + 2*b*c*d + c*d^2))*log(-sqrt(1/2)*(b^2*c - 4*a*c^2 + 2*b*c*d + c*d^2)*sqrt( -(b + d + (b^2*c - 4*a*c^2 + 2*b*c*d + c*d^2)/sqrt(b^2*c^2 - 4*a*c^3 + 2*b *c^2*d + c^2*d^2))/(b^2*c - 4*a*c^2 + 2*b*c*d + c*d^2))/sqrt(b^2*c^2 - 4*a *c^3 + 2*b*c^2*d + c^2*d^2) + x) - 1/2*sqrt(1/2)*sqrt(-(b + d - (b^2*c - 4 *a*c^2 + 2*b*c*d + c*d^2)/sqrt(b^2*c^2 - 4*a*c^3 + 2*b*c^2*d + c^2*d^2))/( b^2*c - 4*a*c^2 + 2*b*c*d + c*d^2))*log(sqrt(1/2)*(b^2*c - 4*a*c^2 + 2*b*c *d + c*d^2)*sqrt(-(b + d - (b^2*c - 4*a*c^2 + 2*b*c*d + c*d^2)/sqrt(b^2*c^ 2 - 4*a*c^3 + 2*b*c^2*d + c^2*d^2))/(b^2*c - 4*a*c^2 + 2*b*c*d + c*d^2))/s qrt(b^2*c^2 - 4*a*c^3 + 2*b*c^2*d + c^2*d^2) + x) + 1/2*sqrt(1/2)*sqrt(-(b + d - (b^2*c - 4*a*c^2 + 2*b*c*d + c*d^2)/sqrt(b^2*c^2 - 4*a*c^3 + 2*b*c^ 2*d + c^2*d^2))/(b^2*c - 4*a*c^2 + 2*b*c*d + c*d^2))*log(-sqrt(1/2)*(b^2*c - 4*a*c^2 + 2*b*c*d + c*d^2)*sqrt(-(b + d - (b^2*c - 4*a*c^2 + 2*b*c*d + c*d^2)/sqrt(b^2*c^2 - 4*a*c^3 + 2*b*c^2*d + c^2*d^2))/(b^2*c - 4*a*c^2 ...
Time = 1.53 (sec) , antiderivative size = 180, normalized size of antiderivative = 0.92 \[ \int \frac {x^2}{a+(b+d) x^2+c x^4} \, dx=\operatorname {RootSum} {\left (t^{4} \cdot \left (256 a^{2} c^{3} - 128 a b^{2} c^{2} - 256 a b c^{2} d - 128 a c^{2} d^{2} + 16 b^{4} c + 64 b^{3} c d + 96 b^{2} c d^{2} + 64 b c d^{3} + 16 c d^{4}\right ) + t^{2} \left (- 16 a b c - 16 a c d + 4 b^{3} + 12 b^{2} d + 12 b d^{2} + 4 d^{3}\right ) + a, \left ( t \mapsto t \log {\left (64 t^{3} a c^{2} - 16 t^{3} b^{2} c - 32 t^{3} b c d - 16 t^{3} c d^{2} - 2 t b - 2 t d + x \right )} \right )\right )} \] Input:
integrate(x**2/(a+(b+d)*x**2+c*x**4),x)
Output:
RootSum(_t**4*(256*a**2*c**3 - 128*a*b**2*c**2 - 256*a*b*c**2*d - 128*a*c* *2*d**2 + 16*b**4*c + 64*b**3*c*d + 96*b**2*c*d**2 + 64*b*c*d**3 + 16*c*d* *4) + _t**2*(-16*a*b*c - 16*a*c*d + 4*b**3 + 12*b**2*d + 12*b*d**2 + 4*d** 3) + a, Lambda(_t, _t*log(64*_t**3*a*c**2 - 16*_t**3*b**2*c - 32*_t**3*b*c *d - 16*_t**3*c*d**2 - 2*_t*b - 2*_t*d + x)))
\[ \int \frac {x^2}{a+(b+d) x^2+c x^4} \, dx=\int { \frac {x^{2}}{c x^{4} + {\left (b + d\right )} x^{2} + a} \,d x } \] Input:
integrate(x^2/(a+(b+d)*x^2+c*x^4),x, algorithm="maxima")
Output:
integrate(x^2/(c*x^4 + (b + d)*x^2 + a), x)
Leaf count of result is larger than twice the leaf count of optimal. 1160 vs. \(2 (161) = 322\).
Time = 2.42 (sec) , antiderivative size = 1160, normalized size of antiderivative = 5.92 \[ \int \frac {x^2}{a+(b+d) x^2+c x^4} \, dx =\text {Too large to display} \] Input:
integrate(x^2/(a+(b+d)*x^2+c*x^4),x, algorithm="giac")
Output:
1/2*(2*b^2*c^2 - 8*a*c^3 + 4*b*c^2*d + 2*c^2*d^2 - sqrt(2)*sqrt(b^2 - 4*a* c + 2*b*d + d^2)*sqrt(b*c + c*d + sqrt(b^2 - 4*a*c + 2*b*d + d^2)*c)*b^2 + 4*sqrt(2)*sqrt(b^2 - 4*a*c + 2*b*d + d^2)*sqrt(b*c + c*d + sqrt(b^2 - 4*a *c + 2*b*d + d^2)*c)*a*c + 2*sqrt(2)*sqrt(b^2 - 4*a*c + 2*b*d + d^2)*sqrt( b*c + c*d + sqrt(b^2 - 4*a*c + 2*b*d + d^2)*c)*b*c - sqrt(2)*sqrt(b^2 - 4* a*c + 2*b*d + d^2)*sqrt(b*c + c*d + sqrt(b^2 - 4*a*c + 2*b*d + d^2)*c)*c^2 - 2*sqrt(2)*sqrt(b^2 - 4*a*c + 2*b*d + d^2)*sqrt(b*c + c*d + sqrt(b^2 - 4 *a*c + 2*b*d + d^2)*c)*b*d + 2*sqrt(2)*sqrt(b^2 - 4*a*c + 2*b*d + d^2)*sqr t(b*c + c*d + sqrt(b^2 - 4*a*c + 2*b*d + d^2)*c)*c*d - sqrt(2)*sqrt(b^2 - 4*a*c + 2*b*d + d^2)*sqrt(b*c + c*d + sqrt(b^2 - 4*a*c + 2*b*d + d^2)*c)*d ^2 - 2*(b^2 - 4*a*c + 2*b*d + d^2)*c^2)*arctan(2*sqrt(1/2)*x/sqrt((b + d + sqrt((b + d)^2 - 4*a*c))/c))/((b^4 - 8*a*b^2*c - 2*b^3*c + 16*a^2*c^2 + 8 *a*b*c^2 + b^2*c^2 - 4*a*c^3 + 4*b^3*d - 16*a*b*c*d - 6*b^2*c*d + 8*a*c^2* d + 2*b*c^2*d + 6*b^2*d^2 - 8*a*c*d^2 - 6*b*c*d^2 + c^2*d^2 + 4*b*d^3 - 2* c*d^3 + d^4)*abs(c)) + 1/2*(2*b^2*c^2 - 8*a*c^3 + 4*b*c^2*d + 2*c^2*d^2 - sqrt(2)*sqrt(b^2 - 4*a*c + 2*b*d + d^2)*sqrt(b*c + c*d - sqrt(b^2 - 4*a*c + 2*b*d + d^2)*c)*b^2 + 4*sqrt(2)*sqrt(b^2 - 4*a*c + 2*b*d + d^2)*sqrt(b*c + c*d - sqrt(b^2 - 4*a*c + 2*b*d + d^2)*c)*a*c + 2*sqrt(2)*sqrt(b^2 - 4*a *c + 2*b*d + d^2)*sqrt(b*c + c*d - sqrt(b^2 - 4*a*c + 2*b*d + d^2)*c)*b*c - sqrt(2)*sqrt(b^2 - 4*a*c + 2*b*d + d^2)*sqrt(b*c + c*d - sqrt(b^2 - 4...
Time = 18.92 (sec) , antiderivative size = 921, normalized size of antiderivative = 4.70 \[ \int \frac {x^2}{a+(b+d) x^2+c x^4} \, dx=2\,\mathrm {atanh}\left (\frac {\left (x\,\left (2\,b^2\,c+4\,b\,c\,d-4\,a\,c^2+2\,c\,d^2\right )-\frac {x\,\left (8\,b^3\,c^2+24\,b^2\,c^2\,d-32\,a\,b\,c^3+24\,b\,c^2\,d^2-32\,a\,c^3\,d+8\,c^2\,d^3\right )\,\left (3\,b\,d^2+3\,b^2\,d+b^3+d^3+\sqrt {{\left (b^2+2\,b\,d+d^2-4\,a\,c\right )}^3}-4\,a\,b\,c-4\,a\,c\,d\right )}{8\,\left (16\,a^2\,c^3-8\,a\,b^2\,c^2-16\,a\,b\,c^2\,d-8\,a\,c^2\,d^2+b^4\,c+4\,b^3\,c\,d+6\,b^2\,c\,d^2+4\,b\,c\,d^3+c\,d^4\right )}\right )\,\sqrt {-\frac {3\,b\,d^2+3\,b^2\,d+b^3+d^3+\sqrt {{\left (b^2+2\,b\,d+d^2-4\,a\,c\right )}^3}-4\,a\,b\,c-4\,a\,c\,d}{8\,\left (16\,a^2\,c^3-8\,a\,b^2\,c^2-16\,a\,b\,c^2\,d-8\,a\,c^2\,d^2+b^4\,c+4\,b^3\,c\,d+6\,b^2\,c\,d^2+4\,b\,c\,d^3+c\,d^4\right )}}}{a\,c}\right )\,\sqrt {-\frac {3\,b\,d^2+3\,b^2\,d+b^3+d^3+\sqrt {{\left (b^2+2\,b\,d+d^2-4\,a\,c\right )}^3}-4\,a\,b\,c-4\,a\,c\,d}{8\,\left (16\,a^2\,c^3-8\,a\,b^2\,c^2-16\,a\,b\,c^2\,d-8\,a\,c^2\,d^2+b^4\,c+4\,b^3\,c\,d+6\,b^2\,c\,d^2+4\,b\,c\,d^3+c\,d^4\right )}}+2\,\mathrm {atanh}\left (\frac {\left (x\,\left (2\,b^2\,c+4\,b\,c\,d-4\,a\,c^2+2\,c\,d^2\right )-\frac {x\,\left (8\,b^3\,c^2+24\,b^2\,c^2\,d-32\,a\,b\,c^3+24\,b\,c^2\,d^2-32\,a\,c^3\,d+8\,c^2\,d^3\right )\,\left (3\,b\,d^2+3\,b^2\,d+b^3+d^3-\sqrt {{\left (b^2+2\,b\,d+d^2-4\,a\,c\right )}^3}-4\,a\,b\,c-4\,a\,c\,d\right )}{8\,\left (16\,a^2\,c^3-8\,a\,b^2\,c^2-16\,a\,b\,c^2\,d-8\,a\,c^2\,d^2+b^4\,c+4\,b^3\,c\,d+6\,b^2\,c\,d^2+4\,b\,c\,d^3+c\,d^4\right )}\right )\,\sqrt {-\frac {3\,b\,d^2+3\,b^2\,d+b^3+d^3-\sqrt {{\left (b^2+2\,b\,d+d^2-4\,a\,c\right )}^3}-4\,a\,b\,c-4\,a\,c\,d}{8\,\left (16\,a^2\,c^3-8\,a\,b^2\,c^2-16\,a\,b\,c^2\,d-8\,a\,c^2\,d^2+b^4\,c+4\,b^3\,c\,d+6\,b^2\,c\,d^2+4\,b\,c\,d^3+c\,d^4\right )}}}{a\,c}\right )\,\sqrt {-\frac {3\,b\,d^2+3\,b^2\,d+b^3+d^3-\sqrt {{\left (b^2+2\,b\,d+d^2-4\,a\,c\right )}^3}-4\,a\,b\,c-4\,a\,c\,d}{8\,\left (16\,a^2\,c^3-8\,a\,b^2\,c^2-16\,a\,b\,c^2\,d-8\,a\,c^2\,d^2+b^4\,c+4\,b^3\,c\,d+6\,b^2\,c\,d^2+4\,b\,c\,d^3+c\,d^4\right )}} \] Input:
int(x^2/(a + x^2*(b + d) + c*x^4),x)
Output:
2*atanh(((x*(2*b^2*c - 4*a*c^2 + 2*c*d^2 + 4*b*c*d) - (x*(8*b^3*c^2 + 8*c^ 2*d^3 + 24*b*c^2*d^2 + 24*b^2*c^2*d - 32*a*b*c^3 - 32*a*c^3*d)*(3*b*d^2 + 3*b^2*d + b^3 + d^3 + ((2*b*d - 4*a*c + b^2 + d^2)^3)^(1/2) - 4*a*b*c - 4* a*c*d))/(8*(b^4*c + c*d^4 + 16*a^2*c^3 - 8*a*b^2*c^2 - 8*a*c^2*d^2 + 6*b^2 *c*d^2 + 4*b*c*d^3 + 4*b^3*c*d - 16*a*b*c^2*d)))*(-(3*b*d^2 + 3*b^2*d + b^ 3 + d^3 + ((2*b*d - 4*a*c + b^2 + d^2)^3)^(1/2) - 4*a*b*c - 4*a*c*d)/(8*(b ^4*c + c*d^4 + 16*a^2*c^3 - 8*a*b^2*c^2 - 8*a*c^2*d^2 + 6*b^2*c*d^2 + 4*b* c*d^3 + 4*b^3*c*d - 16*a*b*c^2*d)))^(1/2))/(a*c))*(-(3*b*d^2 + 3*b^2*d + b ^3 + d^3 + ((2*b*d - 4*a*c + b^2 + d^2)^3)^(1/2) - 4*a*b*c - 4*a*c*d)/(8*( b^4*c + c*d^4 + 16*a^2*c^3 - 8*a*b^2*c^2 - 8*a*c^2*d^2 + 6*b^2*c*d^2 + 4*b *c*d^3 + 4*b^3*c*d - 16*a*b*c^2*d)))^(1/2) + 2*atanh(((x*(2*b^2*c - 4*a*c^ 2 + 2*c*d^2 + 4*b*c*d) - (x*(8*b^3*c^2 + 8*c^2*d^3 + 24*b*c^2*d^2 + 24*b^2 *c^2*d - 32*a*b*c^3 - 32*a*c^3*d)*(3*b*d^2 + 3*b^2*d + b^3 + d^3 - ((2*b*d - 4*a*c + b^2 + d^2)^3)^(1/2) - 4*a*b*c - 4*a*c*d))/(8*(b^4*c + c*d^4 + 1 6*a^2*c^3 - 8*a*b^2*c^2 - 8*a*c^2*d^2 + 6*b^2*c*d^2 + 4*b*c*d^3 + 4*b^3*c* d - 16*a*b*c^2*d)))*(-(3*b*d^2 + 3*b^2*d + b^3 + d^3 - ((2*b*d - 4*a*c + b ^2 + d^2)^3)^(1/2) - 4*a*b*c - 4*a*c*d)/(8*(b^4*c + c*d^4 + 16*a^2*c^3 - 8 *a*b^2*c^2 - 8*a*c^2*d^2 + 6*b^2*c*d^2 + 4*b*c*d^3 + 4*b^3*c*d - 16*a*b*c^ 2*d)))^(1/2))/(a*c))*(-(3*b*d^2 + 3*b^2*d + b^3 + d^3 - ((2*b*d - 4*a*c + b^2 + d^2)^3)^(1/2) - 4*a*b*c - 4*a*c*d)/(8*(b^4*c + c*d^4 + 16*a^2*c^3...
Time = 0.21 (sec) , antiderivative size = 592, normalized size of antiderivative = 3.02 \[ \int \frac {x^2}{a+(b+d) x^2+c x^4} \, dx=\frac {-4 \sqrt {a}\, \sqrt {2 \sqrt {c}\, \sqrt {a}+b +d}\, \mathit {atan} \left (\frac {\sqrt {2 \sqrt {c}\, \sqrt {a}-b -d}-2 \sqrt {c}\, x}{\sqrt {2 \sqrt {c}\, \sqrt {a}+b +d}}\right ) c +2 \sqrt {c}\, \sqrt {2 \sqrt {c}\, \sqrt {a}+b +d}\, \mathit {atan} \left (\frac {\sqrt {2 \sqrt {c}\, \sqrt {a}-b -d}-2 \sqrt {c}\, x}{\sqrt {2 \sqrt {c}\, \sqrt {a}+b +d}}\right ) b +2 \sqrt {c}\, \sqrt {2 \sqrt {c}\, \sqrt {a}+b +d}\, \mathit {atan} \left (\frac {\sqrt {2 \sqrt {c}\, \sqrt {a}-b -d}-2 \sqrt {c}\, x}{\sqrt {2 \sqrt {c}\, \sqrt {a}+b +d}}\right ) d +4 \sqrt {a}\, \sqrt {2 \sqrt {c}\, \sqrt {a}+b +d}\, \mathit {atan} \left (\frac {\sqrt {2 \sqrt {c}\, \sqrt {a}-b -d}+2 \sqrt {c}\, x}{\sqrt {2 \sqrt {c}\, \sqrt {a}+b +d}}\right ) c -2 \sqrt {c}\, \sqrt {2 \sqrt {c}\, \sqrt {a}+b +d}\, \mathit {atan} \left (\frac {\sqrt {2 \sqrt {c}\, \sqrt {a}-b -d}+2 \sqrt {c}\, x}{\sqrt {2 \sqrt {c}\, \sqrt {a}+b +d}}\right ) b -2 \sqrt {c}\, \sqrt {2 \sqrt {c}\, \sqrt {a}+b +d}\, \mathit {atan} \left (\frac {\sqrt {2 \sqrt {c}\, \sqrt {a}-b -d}+2 \sqrt {c}\, x}{\sqrt {2 \sqrt {c}\, \sqrt {a}+b +d}}\right ) d +2 \sqrt {a}\, \sqrt {2 \sqrt {c}\, \sqrt {a}-b -d}\, \mathrm {log}\left (-\sqrt {2 \sqrt {c}\, \sqrt {a}-b -d}\, x +\sqrt {a}+\sqrt {c}\, x^{2}\right ) c -2 \sqrt {a}\, \sqrt {2 \sqrt {c}\, \sqrt {a}-b -d}\, \mathrm {log}\left (\sqrt {2 \sqrt {c}\, \sqrt {a}-b -d}\, x +\sqrt {a}+\sqrt {c}\, x^{2}\right ) c +\sqrt {c}\, \sqrt {2 \sqrt {c}\, \sqrt {a}-b -d}\, \mathrm {log}\left (-\sqrt {2 \sqrt {c}\, \sqrt {a}-b -d}\, x +\sqrt {a}+\sqrt {c}\, x^{2}\right ) b +\sqrt {c}\, \sqrt {2 \sqrt {c}\, \sqrt {a}-b -d}\, \mathrm {log}\left (-\sqrt {2 \sqrt {c}\, \sqrt {a}-b -d}\, x +\sqrt {a}+\sqrt {c}\, x^{2}\right ) d -\sqrt {c}\, \sqrt {2 \sqrt {c}\, \sqrt {a}-b -d}\, \mathrm {log}\left (\sqrt {2 \sqrt {c}\, \sqrt {a}-b -d}\, x +\sqrt {a}+\sqrt {c}\, x^{2}\right ) b -\sqrt {c}\, \sqrt {2 \sqrt {c}\, \sqrt {a}-b -d}\, \mathrm {log}\left (\sqrt {2 \sqrt {c}\, \sqrt {a}-b -d}\, x +\sqrt {a}+\sqrt {c}\, x^{2}\right ) d}{4 c \left (4 a c -b^{2}-2 b d -d^{2}\right )} \] Input:
int(x^2/(a+(b+d)*x^2+c*x^4),x)
Output:
( - 4*sqrt(a)*sqrt(2*sqrt(c)*sqrt(a) + b + d)*atan((sqrt(2*sqrt(c)*sqrt(a) - b - d) - 2*sqrt(c)*x)/sqrt(2*sqrt(c)*sqrt(a) + b + d))*c + 2*sqrt(c)*sq rt(2*sqrt(c)*sqrt(a) + b + d)*atan((sqrt(2*sqrt(c)*sqrt(a) - b - d) - 2*sq rt(c)*x)/sqrt(2*sqrt(c)*sqrt(a) + b + d))*b + 2*sqrt(c)*sqrt(2*sqrt(c)*sqr t(a) + b + d)*atan((sqrt(2*sqrt(c)*sqrt(a) - b - d) - 2*sqrt(c)*x)/sqrt(2* sqrt(c)*sqrt(a) + b + d))*d + 4*sqrt(a)*sqrt(2*sqrt(c)*sqrt(a) + b + d)*at an((sqrt(2*sqrt(c)*sqrt(a) - b - d) + 2*sqrt(c)*x)/sqrt(2*sqrt(c)*sqrt(a) + b + d))*c - 2*sqrt(c)*sqrt(2*sqrt(c)*sqrt(a) + b + d)*atan((sqrt(2*sqrt( c)*sqrt(a) - b - d) + 2*sqrt(c)*x)/sqrt(2*sqrt(c)*sqrt(a) + b + d))*b - 2* sqrt(c)*sqrt(2*sqrt(c)*sqrt(a) + b + d)*atan((sqrt(2*sqrt(c)*sqrt(a) - b - d) + 2*sqrt(c)*x)/sqrt(2*sqrt(c)*sqrt(a) + b + d))*d + 2*sqrt(a)*sqrt(2*s qrt(c)*sqrt(a) - b - d)*log( - sqrt(2*sqrt(c)*sqrt(a) - b - d)*x + sqrt(a) + sqrt(c)*x**2)*c - 2*sqrt(a)*sqrt(2*sqrt(c)*sqrt(a) - b - d)*log(sqrt(2* sqrt(c)*sqrt(a) - b - d)*x + sqrt(a) + sqrt(c)*x**2)*c + sqrt(c)*sqrt(2*sq rt(c)*sqrt(a) - b - d)*log( - sqrt(2*sqrt(c)*sqrt(a) - b - d)*x + sqrt(a) + sqrt(c)*x**2)*b + sqrt(c)*sqrt(2*sqrt(c)*sqrt(a) - b - d)*log( - sqrt(2* sqrt(c)*sqrt(a) - b - d)*x + sqrt(a) + sqrt(c)*x**2)*d - sqrt(c)*sqrt(2*sq rt(c)*sqrt(a) - b - d)*log(sqrt(2*sqrt(c)*sqrt(a) - b - d)*x + sqrt(a) + s qrt(c)*x**2)*b - sqrt(c)*sqrt(2*sqrt(c)*sqrt(a) - b - d)*log(sqrt(2*sqrt(c )*sqrt(a) - b - d)*x + sqrt(a) + sqrt(c)*x**2)*d)/(4*c*(4*a*c - b**2 - ...