\(\int \frac {x^{9/2}}{(a+c x^4)^3} \, dx\) [153]

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

Optimal result

Integrand size = 15, antiderivative size = 273 \[ \int \frac {x^{9/2}}{\left (a+c x^4\right )^3} \, dx=-\frac {x^{3/2}}{8 c \left (a+c x^4\right )^2}+\frac {3 x^{3/2}}{64 a c \left (a+c x^4\right )}-\frac {15 \arctan \left (1-\frac {\sqrt {2} \sqrt [8]{c} \sqrt {x}}{\sqrt [8]{-a}}\right )}{256 \sqrt {2} (-a)^{13/8} c^{11/8}}+\frac {15 \arctan \left (1+\frac {\sqrt {2} \sqrt [8]{c} \sqrt {x}}{\sqrt [8]{-a}}\right )}{256 \sqrt {2} (-a)^{13/8} c^{11/8}}-\frac {15 \arctan \left (\frac {\sqrt [8]{c} \sqrt {x}}{\sqrt [8]{-a}}\right )}{256 (-a)^{13/8} c^{11/8}}+\frac {15 \text {arctanh}\left (\frac {\sqrt [8]{c} \sqrt {x}}{\sqrt [8]{-a}}\right )}{256 (-a)^{13/8} c^{11/8}}-\frac {15 \text {arctanh}\left (\frac {\sqrt {2} \sqrt [8]{-a} \sqrt [8]{c} \sqrt {x}}{\sqrt [4]{-a}+\sqrt [4]{c} x}\right )}{256 \sqrt {2} (-a)^{13/8} c^{11/8}} \] Output:

-1/8*x^(3/2)/c/(c*x^4+a)^2+3/64*x^(3/2)/a/c/(c*x^4+a)+15/512*arctan(-1+2^( 
1/2)*c^(1/8)*x^(1/2)/(-a)^(1/8))*2^(1/2)/(-a)^(13/8)/c^(11/8)+15/512*arcta 
n(1+2^(1/2)*c^(1/8)*x^(1/2)/(-a)^(1/8))*2^(1/2)/(-a)^(13/8)/c^(11/8)-15/25 
6*arctan(c^(1/8)*x^(1/2)/(-a)^(1/8))/(-a)^(13/8)/c^(11/8)+15/256*arctanh(c 
^(1/8)*x^(1/2)/(-a)^(1/8))/(-a)^(13/8)/c^(11/8)-15/512*arctanh(2^(1/2)*(-a 
)^(1/8)*c^(1/8)*x^(1/2)/((-a)^(1/4)+c^(1/4)*x))*2^(1/2)/(-a)^(13/8)/c^(11/ 
8)
 

Mathematica [A] (verified)

Time = 1.54 (sec) , antiderivative size = 277, normalized size of antiderivative = 1.01 \[ \int \frac {x^{9/2}}{\left (a+c x^4\right )^3} \, dx=\frac {\frac {8 a^{5/8} c^{3/8} x^{3/2} \left (-5 a+3 c x^4\right )}{\left (a+c x^4\right )^2}+15 \sqrt {2-\sqrt {2}} \arctan \left (\frac {\sqrt {1-\frac {1}{\sqrt {2}}} \left (\sqrt [4]{a}-\sqrt [4]{c} x\right )}{\sqrt [8]{a} \sqrt [8]{c} \sqrt {x}}\right )-15 \sqrt {2+\sqrt {2}} \arctan \left (\frac {\sqrt {1+\frac {1}{\sqrt {2}}} \left (\sqrt [4]{a}-\sqrt [4]{c} x\right )}{\sqrt [8]{a} \sqrt [8]{c} \sqrt {x}}\right )+15 \sqrt {2-\sqrt {2}} \text {arctanh}\left (\frac {\sqrt {2+\sqrt {2}} \sqrt [8]{a} \sqrt [8]{c} \sqrt {x}}{\sqrt [4]{a}+\sqrt [4]{c} x}\right )-15 \sqrt {2+\sqrt {2}} \text {arctanh}\left (\frac {\sqrt [8]{a} \sqrt [8]{c} \sqrt {-\left (\left (-2+\sqrt {2}\right ) x\right )}}{\sqrt [4]{a}+\sqrt [4]{c} x}\right )}{512 a^{13/8} c^{11/8}} \] Input:

Integrate[x^(9/2)/(a + c*x^4)^3,x]
 

Output:

((8*a^(5/8)*c^(3/8)*x^(3/2)*(-5*a + 3*c*x^4))/(a + c*x^4)^2 + 15*Sqrt[2 - 
Sqrt[2]]*ArcTan[(Sqrt[1 - 1/Sqrt[2]]*(a^(1/4) - c^(1/4)*x))/(a^(1/8)*c^(1/ 
8)*Sqrt[x])] - 15*Sqrt[2 + Sqrt[2]]*ArcTan[(Sqrt[1 + 1/Sqrt[2]]*(a^(1/4) - 
 c^(1/4)*x))/(a^(1/8)*c^(1/8)*Sqrt[x])] + 15*Sqrt[2 - Sqrt[2]]*ArcTanh[(Sq 
rt[2 + Sqrt[2]]*a^(1/8)*c^(1/8)*Sqrt[x])/(a^(1/4) + c^(1/4)*x)] - 15*Sqrt[ 
2 + Sqrt[2]]*ArcTanh[(a^(1/8)*c^(1/8)*Sqrt[-((-2 + Sqrt[2])*x)])/(a^(1/4) 
+ c^(1/4)*x)])/(512*a^(13/8)*c^(11/8))
 

Rubi [A] (verified)

Time = 1.06 (sec) , antiderivative size = 384, normalized size of antiderivative = 1.41, number of steps used = 16, number of rules used = 15, \(\frac {\text {number of rules}}{\text {integrand size}}\) = 1.000, Rules used = {817, 819, 851, 829, 826, 827, 218, 221, 1476, 1082, 217, 1479, 25, 27, 1103}

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^{9/2}}{\left (a+c x^4\right )^3} \, dx\)

\(\Big \downarrow \) 817

\(\displaystyle \frac {3 \int \frac {\sqrt {x}}{\left (c x^4+a\right )^2}dx}{16 c}-\frac {x^{3/2}}{8 c \left (a+c x^4\right )^2}\)

\(\Big \downarrow \) 819

\(\displaystyle \frac {3 \left (\frac {5 \int \frac {\sqrt {x}}{c x^4+a}dx}{8 a}+\frac {x^{3/2}}{4 a \left (a+c x^4\right )}\right )}{16 c}-\frac {x^{3/2}}{8 c \left (a+c x^4\right )^2}\)

\(\Big \downarrow \) 851

\(\displaystyle \frac {3 \left (\frac {5 \int \frac {x}{c x^4+a}d\sqrt {x}}{4 a}+\frac {x^{3/2}}{4 a \left (a+c x^4\right )}\right )}{16 c}-\frac {x^{3/2}}{8 c \left (a+c x^4\right )^2}\)

\(\Big \downarrow \) 829

\(\displaystyle \frac {3 \left (\frac {5 \left (-\frac {\int \frac {x}{\sqrt {-a}-\sqrt {c} x^2}d\sqrt {x}}{2 \sqrt {-a}}-\frac {\int \frac {x}{\sqrt {c} x^2+\sqrt {-a}}d\sqrt {x}}{2 \sqrt {-a}}\right )}{4 a}+\frac {x^{3/2}}{4 a \left (a+c x^4\right )}\right )}{16 c}-\frac {x^{3/2}}{8 c \left (a+c x^4\right )^2}\)

\(\Big \downarrow \) 826

\(\displaystyle \frac {3 \left (\frac {5 \left (-\frac {\int \frac {x}{\sqrt {-a}-\sqrt {c} x^2}d\sqrt {x}}{2 \sqrt {-a}}-\frac {\frac {\int \frac {\sqrt [4]{c} x+\sqrt [4]{-a}}{\sqrt {c} x^2+\sqrt {-a}}d\sqrt {x}}{2 \sqrt [4]{c}}-\frac {\int \frac {\sqrt [4]{-a}-\sqrt [4]{c} x}{\sqrt {c} x^2+\sqrt {-a}}d\sqrt {x}}{2 \sqrt [4]{c}}}{2 \sqrt {-a}}\right )}{4 a}+\frac {x^{3/2}}{4 a \left (a+c x^4\right )}\right )}{16 c}-\frac {x^{3/2}}{8 c \left (a+c x^4\right )^2}\)

\(\Big \downarrow \) 827

\(\displaystyle \frac {3 \left (\frac {5 \left (-\frac {\frac {\int \frac {\sqrt [4]{c} x+\sqrt [4]{-a}}{\sqrt {c} x^2+\sqrt {-a}}d\sqrt {x}}{2 \sqrt [4]{c}}-\frac {\int \frac {\sqrt [4]{-a}-\sqrt [4]{c} x}{\sqrt {c} x^2+\sqrt {-a}}d\sqrt {x}}{2 \sqrt [4]{c}}}{2 \sqrt {-a}}-\frac {\frac {\int \frac {1}{\sqrt [4]{-a}-\sqrt [4]{c} x}d\sqrt {x}}{2 \sqrt [4]{c}}-\frac {\int \frac {1}{\sqrt [4]{c} x+\sqrt [4]{-a}}d\sqrt {x}}{2 \sqrt [4]{c}}}{2 \sqrt {-a}}\right )}{4 a}+\frac {x^{3/2}}{4 a \left (a+c x^4\right )}\right )}{16 c}-\frac {x^{3/2}}{8 c \left (a+c x^4\right )^2}\)

\(\Big \downarrow \) 218

\(\displaystyle \frac {3 \left (\frac {5 \left (-\frac {\frac {\int \frac {1}{\sqrt [4]{-a}-\sqrt [4]{c} x}d\sqrt {x}}{2 \sqrt [4]{c}}-\frac {\arctan \left (\frac {\sqrt [8]{c} \sqrt {x}}{\sqrt [8]{-a}}\right )}{2 \sqrt [8]{-a} c^{3/8}}}{2 \sqrt {-a}}-\frac {\frac {\int \frac {\sqrt [4]{c} x+\sqrt [4]{-a}}{\sqrt {c} x^2+\sqrt {-a}}d\sqrt {x}}{2 \sqrt [4]{c}}-\frac {\int \frac {\sqrt [4]{-a}-\sqrt [4]{c} x}{\sqrt {c} x^2+\sqrt {-a}}d\sqrt {x}}{2 \sqrt [4]{c}}}{2 \sqrt {-a}}\right )}{4 a}+\frac {x^{3/2}}{4 a \left (a+c x^4\right )}\right )}{16 c}-\frac {x^{3/2}}{8 c \left (a+c x^4\right )^2}\)

\(\Big \downarrow \) 221

\(\displaystyle \frac {3 \left (\frac {5 \left (-\frac {\frac {\int \frac {\sqrt [4]{c} x+\sqrt [4]{-a}}{\sqrt {c} x^2+\sqrt {-a}}d\sqrt {x}}{2 \sqrt [4]{c}}-\frac {\int \frac {\sqrt [4]{-a}-\sqrt [4]{c} x}{\sqrt {c} x^2+\sqrt {-a}}d\sqrt {x}}{2 \sqrt [4]{c}}}{2 \sqrt {-a}}-\frac {\frac {\text {arctanh}\left (\frac {\sqrt [8]{c} \sqrt {x}}{\sqrt [8]{-a}}\right )}{2 \sqrt [8]{-a} c^{3/8}}-\frac {\arctan \left (\frac {\sqrt [8]{c} \sqrt {x}}{\sqrt [8]{-a}}\right )}{2 \sqrt [8]{-a} c^{3/8}}}{2 \sqrt {-a}}\right )}{4 a}+\frac {x^{3/2}}{4 a \left (a+c x^4\right )}\right )}{16 c}-\frac {x^{3/2}}{8 c \left (a+c x^4\right )^2}\)

\(\Big \downarrow \) 1476

\(\displaystyle \frac {3 \left (\frac {5 \left (-\frac {\frac {\frac {\int \frac {1}{x-\frac {\sqrt {2} \sqrt [8]{-a} \sqrt {x}}{\sqrt [8]{c}}+\frac {\sqrt [4]{-a}}{\sqrt [4]{c}}}d\sqrt {x}}{2 \sqrt [4]{c}}+\frac {\int \frac {1}{x+\frac {\sqrt {2} \sqrt [8]{-a} \sqrt {x}}{\sqrt [8]{c}}+\frac {\sqrt [4]{-a}}{\sqrt [4]{c}}}d\sqrt {x}}{2 \sqrt [4]{c}}}{2 \sqrt [4]{c}}-\frac {\int \frac {\sqrt [4]{-a}-\sqrt [4]{c} x}{\sqrt {c} x^2+\sqrt {-a}}d\sqrt {x}}{2 \sqrt [4]{c}}}{2 \sqrt {-a}}-\frac {\frac {\text {arctanh}\left (\frac {\sqrt [8]{c} \sqrt {x}}{\sqrt [8]{-a}}\right )}{2 \sqrt [8]{-a} c^{3/8}}-\frac {\arctan \left (\frac {\sqrt [8]{c} \sqrt {x}}{\sqrt [8]{-a}}\right )}{2 \sqrt [8]{-a} c^{3/8}}}{2 \sqrt {-a}}\right )}{4 a}+\frac {x^{3/2}}{4 a \left (a+c x^4\right )}\right )}{16 c}-\frac {x^{3/2}}{8 c \left (a+c x^4\right )^2}\)

\(\Big \downarrow \) 1082

\(\displaystyle \frac {3 \left (\frac {5 \left (-\frac {\frac {\frac {\int \frac {1}{-x-1}d\left (1-\frac {\sqrt {2} \sqrt [8]{c} \sqrt {x}}{\sqrt [8]{-a}}\right )}{\sqrt {2} \sqrt [8]{-a} \sqrt [8]{c}}-\frac {\int \frac {1}{-x-1}d\left (\frac {\sqrt {2} \sqrt [8]{c} \sqrt {x}}{\sqrt [8]{-a}}+1\right )}{\sqrt {2} \sqrt [8]{-a} \sqrt [8]{c}}}{2 \sqrt [4]{c}}-\frac {\int \frac {\sqrt [4]{-a}-\sqrt [4]{c} x}{\sqrt {c} x^2+\sqrt {-a}}d\sqrt {x}}{2 \sqrt [4]{c}}}{2 \sqrt {-a}}-\frac {\frac {\text {arctanh}\left (\frac {\sqrt [8]{c} \sqrt {x}}{\sqrt [8]{-a}}\right )}{2 \sqrt [8]{-a} c^{3/8}}-\frac {\arctan \left (\frac {\sqrt [8]{c} \sqrt {x}}{\sqrt [8]{-a}}\right )}{2 \sqrt [8]{-a} c^{3/8}}}{2 \sqrt {-a}}\right )}{4 a}+\frac {x^{3/2}}{4 a \left (a+c x^4\right )}\right )}{16 c}-\frac {x^{3/2}}{8 c \left (a+c x^4\right )^2}\)

\(\Big \downarrow \) 217

\(\displaystyle \frac {3 \left (\frac {5 \left (-\frac {\frac {\frac {\arctan \left (\frac {\sqrt {2} \sqrt [8]{c} \sqrt {x}}{\sqrt [8]{-a}}+1\right )}{\sqrt {2} \sqrt [8]{-a} \sqrt [8]{c}}-\frac {\arctan \left (1-\frac {\sqrt {2} \sqrt [8]{c} \sqrt {x}}{\sqrt [8]{-a}}\right )}{\sqrt {2} \sqrt [8]{-a} \sqrt [8]{c}}}{2 \sqrt [4]{c}}-\frac {\int \frac {\sqrt [4]{-a}-\sqrt [4]{c} x}{\sqrt {c} x^2+\sqrt {-a}}d\sqrt {x}}{2 \sqrt [4]{c}}}{2 \sqrt {-a}}-\frac {\frac {\text {arctanh}\left (\frac {\sqrt [8]{c} \sqrt {x}}{\sqrt [8]{-a}}\right )}{2 \sqrt [8]{-a} c^{3/8}}-\frac {\arctan \left (\frac {\sqrt [8]{c} \sqrt {x}}{\sqrt [8]{-a}}\right )}{2 \sqrt [8]{-a} c^{3/8}}}{2 \sqrt {-a}}\right )}{4 a}+\frac {x^{3/2}}{4 a \left (a+c x^4\right )}\right )}{16 c}-\frac {x^{3/2}}{8 c \left (a+c x^4\right )^2}\)

\(\Big \downarrow \) 1479

\(\displaystyle \frac {3 \left (\frac {5 \left (-\frac {\frac {\frac {\arctan \left (\frac {\sqrt {2} \sqrt [8]{c} \sqrt {x}}{\sqrt [8]{-a}}+1\right )}{\sqrt {2} \sqrt [8]{-a} \sqrt [8]{c}}-\frac {\arctan \left (1-\frac {\sqrt {2} \sqrt [8]{c} \sqrt {x}}{\sqrt [8]{-a}}\right )}{\sqrt {2} \sqrt [8]{-a} \sqrt [8]{c}}}{2 \sqrt [4]{c}}-\frac {-\frac {\int -\frac {\sqrt {2} \sqrt [8]{-a}-2 \sqrt [8]{c} \sqrt {x}}{\sqrt [8]{c} \left (x-\frac {\sqrt {2} \sqrt [8]{-a} \sqrt {x}}{\sqrt [8]{c}}+\frac {\sqrt [4]{-a}}{\sqrt [4]{c}}\right )}d\sqrt {x}}{2 \sqrt {2} \sqrt [8]{-a} \sqrt [8]{c}}-\frac {\int -\frac {\sqrt {2} \left (\sqrt {2} \sqrt [8]{c} \sqrt {x}+\sqrt [8]{-a}\right )}{\sqrt [8]{c} \left (x+\frac {\sqrt {2} \sqrt [8]{-a} \sqrt {x}}{\sqrt [8]{c}}+\frac {\sqrt [4]{-a}}{\sqrt [4]{c}}\right )}d\sqrt {x}}{2 \sqrt {2} \sqrt [8]{-a} \sqrt [8]{c}}}{2 \sqrt [4]{c}}}{2 \sqrt {-a}}-\frac {\frac {\text {arctanh}\left (\frac {\sqrt [8]{c} \sqrt {x}}{\sqrt [8]{-a}}\right )}{2 \sqrt [8]{-a} c^{3/8}}-\frac {\arctan \left (\frac {\sqrt [8]{c} \sqrt {x}}{\sqrt [8]{-a}}\right )}{2 \sqrt [8]{-a} c^{3/8}}}{2 \sqrt {-a}}\right )}{4 a}+\frac {x^{3/2}}{4 a \left (a+c x^4\right )}\right )}{16 c}-\frac {x^{3/2}}{8 c \left (a+c x^4\right )^2}\)

\(\Big \downarrow \) 25

\(\displaystyle \frac {3 \left (\frac {5 \left (-\frac {\frac {\frac {\arctan \left (\frac {\sqrt {2} \sqrt [8]{c} \sqrt {x}}{\sqrt [8]{-a}}+1\right )}{\sqrt {2} \sqrt [8]{-a} \sqrt [8]{c}}-\frac {\arctan \left (1-\frac {\sqrt {2} \sqrt [8]{c} \sqrt {x}}{\sqrt [8]{-a}}\right )}{\sqrt {2} \sqrt [8]{-a} \sqrt [8]{c}}}{2 \sqrt [4]{c}}-\frac {\frac {\int \frac {\sqrt {2} \sqrt [8]{-a}-2 \sqrt [8]{c} \sqrt {x}}{\sqrt [8]{c} \left (x-\frac {\sqrt {2} \sqrt [8]{-a} \sqrt {x}}{\sqrt [8]{c}}+\frac {\sqrt [4]{-a}}{\sqrt [4]{c}}\right )}d\sqrt {x}}{2 \sqrt {2} \sqrt [8]{-a} \sqrt [8]{c}}+\frac {\int \frac {\sqrt {2} \left (\sqrt {2} \sqrt [8]{c} \sqrt {x}+\sqrt [8]{-a}\right )}{\sqrt [8]{c} \left (x+\frac {\sqrt {2} \sqrt [8]{-a} \sqrt {x}}{\sqrt [8]{c}}+\frac {\sqrt [4]{-a}}{\sqrt [4]{c}}\right )}d\sqrt {x}}{2 \sqrt {2} \sqrt [8]{-a} \sqrt [8]{c}}}{2 \sqrt [4]{c}}}{2 \sqrt {-a}}-\frac {\frac {\text {arctanh}\left (\frac {\sqrt [8]{c} \sqrt {x}}{\sqrt [8]{-a}}\right )}{2 \sqrt [8]{-a} c^{3/8}}-\frac {\arctan \left (\frac {\sqrt [8]{c} \sqrt {x}}{\sqrt [8]{-a}}\right )}{2 \sqrt [8]{-a} c^{3/8}}}{2 \sqrt {-a}}\right )}{4 a}+\frac {x^{3/2}}{4 a \left (a+c x^4\right )}\right )}{16 c}-\frac {x^{3/2}}{8 c \left (a+c x^4\right )^2}\)

\(\Big \downarrow \) 27

\(\displaystyle \frac {3 \left (\frac {5 \left (-\frac {\frac {\frac {\arctan \left (\frac {\sqrt {2} \sqrt [8]{c} \sqrt {x}}{\sqrt [8]{-a}}+1\right )}{\sqrt {2} \sqrt [8]{-a} \sqrt [8]{c}}-\frac {\arctan \left (1-\frac {\sqrt {2} \sqrt [8]{c} \sqrt {x}}{\sqrt [8]{-a}}\right )}{\sqrt {2} \sqrt [8]{-a} \sqrt [8]{c}}}{2 \sqrt [4]{c}}-\frac {\frac {\int \frac {\sqrt {2} \sqrt [8]{-a}-2 \sqrt [8]{c} \sqrt {x}}{x-\frac {\sqrt {2} \sqrt [8]{-a} \sqrt {x}}{\sqrt [8]{c}}+\frac {\sqrt [4]{-a}}{\sqrt [4]{c}}}d\sqrt {x}}{2 \sqrt {2} \sqrt [8]{-a} \sqrt [4]{c}}+\frac {\int \frac {\sqrt {2} \sqrt [8]{c} \sqrt {x}+\sqrt [8]{-a}}{x+\frac {\sqrt {2} \sqrt [8]{-a} \sqrt {x}}{\sqrt [8]{c}}+\frac {\sqrt [4]{-a}}{\sqrt [4]{c}}}d\sqrt {x}}{2 \sqrt [8]{-a} \sqrt [4]{c}}}{2 \sqrt [4]{c}}}{2 \sqrt {-a}}-\frac {\frac {\text {arctanh}\left (\frac {\sqrt [8]{c} \sqrt {x}}{\sqrt [8]{-a}}\right )}{2 \sqrt [8]{-a} c^{3/8}}-\frac {\arctan \left (\frac {\sqrt [8]{c} \sqrt {x}}{\sqrt [8]{-a}}\right )}{2 \sqrt [8]{-a} c^{3/8}}}{2 \sqrt {-a}}\right )}{4 a}+\frac {x^{3/2}}{4 a \left (a+c x^4\right )}\right )}{16 c}-\frac {x^{3/2}}{8 c \left (a+c x^4\right )^2}\)

\(\Big \downarrow \) 1103

\(\displaystyle \frac {3 \left (\frac {5 \left (-\frac {\frac {\text {arctanh}\left (\frac {\sqrt [8]{c} \sqrt {x}}{\sqrt [8]{-a}}\right )}{2 \sqrt [8]{-a} c^{3/8}}-\frac {\arctan \left (\frac {\sqrt [8]{c} \sqrt {x}}{\sqrt [8]{-a}}\right )}{2 \sqrt [8]{-a} c^{3/8}}}{2 \sqrt {-a}}-\frac {\frac {\frac {\arctan \left (\frac {\sqrt {2} \sqrt [8]{c} \sqrt {x}}{\sqrt [8]{-a}}+1\right )}{\sqrt {2} \sqrt [8]{-a} \sqrt [8]{c}}-\frac {\arctan \left (1-\frac {\sqrt {2} \sqrt [8]{c} \sqrt {x}}{\sqrt [8]{-a}}\right )}{\sqrt {2} \sqrt [8]{-a} \sqrt [8]{c}}}{2 \sqrt [4]{c}}-\frac {\frac {\log \left (\sqrt {2} \sqrt [8]{-a} \sqrt [8]{c} \sqrt {x}+\sqrt [4]{-a}+\sqrt [4]{c} x\right )}{2 \sqrt {2} \sqrt [8]{-a} \sqrt [8]{c}}-\frac {\log \left (-\sqrt {2} \sqrt [8]{-a} \sqrt [8]{c} \sqrt {x}+\sqrt [4]{-a}+\sqrt [4]{c} x\right )}{2 \sqrt {2} \sqrt [8]{-a} \sqrt [8]{c}}}{2 \sqrt [4]{c}}}{2 \sqrt {-a}}\right )}{4 a}+\frac {x^{3/2}}{4 a \left (a+c x^4\right )}\right )}{16 c}-\frac {x^{3/2}}{8 c \left (a+c x^4\right )^2}\)

Input:

Int[x^(9/2)/(a + c*x^4)^3,x]
 

Output:

-1/8*x^(3/2)/(c*(a + c*x^4)^2) + (3*(x^(3/2)/(4*a*(a + c*x^4)) + (5*(-1/2* 
(-1/2*ArcTan[(c^(1/8)*Sqrt[x])/(-a)^(1/8)]/((-a)^(1/8)*c^(3/8)) + ArcTanh[ 
(c^(1/8)*Sqrt[x])/(-a)^(1/8)]/(2*(-a)^(1/8)*c^(3/8)))/Sqrt[-a] - ((-(ArcTa 
n[1 - (Sqrt[2]*c^(1/8)*Sqrt[x])/(-a)^(1/8)]/(Sqrt[2]*(-a)^(1/8)*c^(1/8))) 
+ ArcTan[1 + (Sqrt[2]*c^(1/8)*Sqrt[x])/(-a)^(1/8)]/(Sqrt[2]*(-a)^(1/8)*c^( 
1/8)))/(2*c^(1/4)) - (-1/2*Log[(-a)^(1/4) - Sqrt[2]*(-a)^(1/8)*c^(1/8)*Sqr 
t[x] + c^(1/4)*x]/(Sqrt[2]*(-a)^(1/8)*c^(1/8)) + Log[(-a)^(1/4) + Sqrt[2]* 
(-a)^(1/8)*c^(1/8)*Sqrt[x] + c^(1/4)*x]/(2*Sqrt[2]*(-a)^(1/8)*c^(1/8)))/(2 
*c^(1/4)))/(2*Sqrt[-a])))/(4*a)))/(16*c)
 

Defintions of rubi rules used

rule 25
Int[-(Fx_), x_Symbol] :> Simp[Identity[-1]   Int[Fx, x], x]
 

rule 27
Int[(a_)*(Fx_), x_Symbol] :> Simp[a   Int[Fx, x], x] /; FreeQ[a, x] &&  !Ma 
tchQ[Fx, (b_)*(Gx_) /; FreeQ[b, x]]
 

rule 217
Int[((a_) + (b_.)*(x_)^2)^(-1), x_Symbol] :> Simp[(-(Rt[-a, 2]*Rt[-b, 2])^( 
-1))*ArcTan[Rt[-b, 2]*(x/Rt[-a, 2])], x] /; FreeQ[{a, b}, x] && PosQ[a/b] & 
& (LtQ[a, 0] || LtQ[b, 0])
 

rule 218
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]
 

rule 221
Int[((a_) + (b_.)*(x_)^2)^(-1), x_Symbol] :> Simp[(Rt[-a/b, 2]/a)*ArcTanh[x 
/Rt[-a/b, 2]], x] /; FreeQ[{a, b}, x] && NegQ[a/b]
 

rule 817
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*n*(p + 1))), x] - Simp[c^n 
*((m - n + 1)/(b*n*(p + 1)))   Int[(c*x)^(m - n)*(a + b*x^n)^(p + 1), x], x 
] /; FreeQ[{a, b, c}, x] && IGtQ[n, 0] && LtQ[p, -1] && GtQ[m + 1, n] &&  ! 
ILtQ[(m + n*(p + 1) + 1)/n, 0] && IntBinomialQ[a, b, c, n, m, p, x]
 

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

rule 826
Int[(x_)^2/((a_) + (b_.)*(x_)^4), x_Symbol] :> With[{r = Numerator[Rt[a/b, 
2]], s = Denominator[Rt[a/b, 2]]}, Simp[1/(2*s)   Int[(r + s*x^2)/(a + b*x^ 
4), x], x] - Simp[1/(2*s)   Int[(r - s*x^2)/(a + b*x^4), x], x]] /; FreeQ[{ 
a, b}, x] && (GtQ[a/b, 0] || (PosQ[a/b] && AtomQ[SplitProduct[SumBaseQ, a]] 
 && AtomQ[SplitProduct[SumBaseQ, b]]))
 

rule 827
Int[(x_)^2/((a_) + (b_.)*(x_)^4), x_Symbol] :> With[{r = Numerator[Rt[-a/b, 
 2]], s = Denominator[Rt[-a/b, 2]]}, Simp[s/(2*b)   Int[1/(r + s*x^2), x], 
x] - Simp[s/(2*b)   Int[1/(r - s*x^2), x], x]] /; FreeQ[{a, b}, x] &&  !GtQ 
[a/b, 0]
 

rule 829
Int[(x_)^(m_)/((a_) + (b_.)*(x_)^(n_)), x_Symbol] :> With[{r = Numerator[Rt 
[-a/b, 2]], s = Denominator[Rt[-a/b, 2]]}, Simp[r/(2*a)   Int[x^m/(r + s*x^ 
(n/2)), x], x] + Simp[r/(2*a)   Int[x^m/(r - s*x^(n/2)), x], x]] /; FreeQ[{ 
a, b}, x] && IGtQ[n/4, 0] && IGtQ[m, 0] && LtQ[m, n/2] &&  !GtQ[a/b, 0]
 

rule 851
Int[((c_.)*(x_))^(m_)*((a_) + (b_.)*(x_)^(n_))^(p_), x_Symbol] :> With[{k = 
 Denominator[m]}, Simp[k/c   Subst[Int[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 1082
Int[((a_) + (b_.)*(x_) + (c_.)*(x_)^2)^(-1), x_Symbol] :> With[{q = 1 - 4*S 
implify[a*(c/b^2)]}, Simp[-2/b   Subst[Int[1/(q - x^2), x], x, 1 + 2*c*(x/b 
)], x] /; RationalQ[q] && (EqQ[q^2, 1] ||  !RationalQ[b^2 - 4*a*c])] /; Fre 
eQ[{a, b, c}, x]
 

rule 1103
Int[((d_) + (e_.)*(x_))/((a_.) + (b_.)*(x_) + (c_.)*(x_)^2), x_Symbol] :> S 
imp[d*(Log[RemoveContent[a + b*x + c*x^2, x]]/b), x] /; FreeQ[{a, b, c, d, 
e}, x] && EqQ[2*c*d - b*e, 0]
 

rule 1476
Int[((d_) + (e_.)*(x_)^2)/((a_) + (c_.)*(x_)^4), x_Symbol] :> With[{q = Rt[ 
2*(d/e), 2]}, Simp[e/(2*c)   Int[1/Simp[d/e + q*x + x^2, x], x], x] + Simp[ 
e/(2*c)   Int[1/Simp[d/e - q*x + x^2, x], x], x]] /; FreeQ[{a, c, d, e}, x] 
 && EqQ[c*d^2 - a*e^2, 0] && PosQ[d*e]
 

rule 1479
Int[((d_) + (e_.)*(x_)^2)/((a_) + (c_.)*(x_)^4), x_Symbol] :> With[{q = Rt[ 
-2*(d/e), 2]}, Simp[e/(2*c*q)   Int[(q - 2*x)/Simp[d/e + q*x - x^2, x], x], 
 x] + Simp[e/(2*c*q)   Int[(q + 2*x)/Simp[d/e - q*x - x^2, x], x], x]] /; F 
reeQ[{a, c, d, e}, x] && EqQ[c*d^2 - a*e^2, 0] && NegQ[d*e]
 
Maple [C] (verified)

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

Time = 0.76 (sec) , antiderivative size = 61, normalized size of antiderivative = 0.22

method result size
derivativedivides \(\frac {-\frac {5 x^{\frac {3}{2}}}{64 c}+\frac {3 x^{\frac {11}{2}}}{64 a}}{\left (c \,x^{4}+a \right )^{2}}+\frac {15 \left (\munderset {\textit {\_R} =\operatorname {RootOf}\left (c \,\textit {\_Z}^{8}+a \right )}{\sum }\frac {\ln \left (\sqrt {x}-\textit {\_R} \right )}{\textit {\_R}^{5}}\right )}{512 a \,c^{2}}\) \(61\)
default \(\frac {-\frac {5 x^{\frac {3}{2}}}{64 c}+\frac {3 x^{\frac {11}{2}}}{64 a}}{\left (c \,x^{4}+a \right )^{2}}+\frac {15 \left (\munderset {\textit {\_R} =\operatorname {RootOf}\left (c \,\textit {\_Z}^{8}+a \right )}{\sum }\frac {\ln \left (\sqrt {x}-\textit {\_R} \right )}{\textit {\_R}^{5}}\right )}{512 a \,c^{2}}\) \(61\)

Input:

int(x^(9/2)/(c*x^4+a)^3,x,method=_RETURNVERBOSE)
 

Output:

2*(-5/128*x^(3/2)/c+3/128/a*x^(11/2))/(c*x^4+a)^2+15/512/a/c^2*sum(1/_R^5* 
ln(x^(1/2)-_R),_R=RootOf(_Z^8*c+a))
 

Fricas [C] (verification not implemented)

Result contains complex when optimal does not.

Time = 0.11 (sec) , antiderivative size = 561, normalized size of antiderivative = 2.05 \[ \int \frac {x^{9/2}}{\left (a+c x^4\right )^3} \, dx =\text {Too large to display} \] Input:

integrate(x^(9/2)/(c*x^4+a)^3,x, algorithm="fricas")
 

Output:

-1/1024*(15*sqrt(2)*((I - 1)*a*c^3*x^8 + (2*I - 2)*a^2*c^2*x^4 + (I - 1)*a 
^3*c)*(-1/(a^13*c^11))^(1/8)*log((1/2*I + 1/2)*sqrt(2)*a^5*c^4*(-1/(a^13*c 
^11))^(3/8) + sqrt(x)) + 15*sqrt(2)*(-(I + 1)*a*c^3*x^8 - (2*I + 2)*a^2*c^ 
2*x^4 - (I + 1)*a^3*c)*(-1/(a^13*c^11))^(1/8)*log(-(1/2*I - 1/2)*sqrt(2)*a 
^5*c^4*(-1/(a^13*c^11))^(3/8) + sqrt(x)) + 15*sqrt(2)*((I + 1)*a*c^3*x^8 + 
 (2*I + 2)*a^2*c^2*x^4 + (I + 1)*a^3*c)*(-1/(a^13*c^11))^(1/8)*log((1/2*I 
- 1/2)*sqrt(2)*a^5*c^4*(-1/(a^13*c^11))^(3/8) + sqrt(x)) + 15*sqrt(2)*(-(I 
 - 1)*a*c^3*x^8 - (2*I - 2)*a^2*c^2*x^4 - (I - 1)*a^3*c)*(-1/(a^13*c^11))^ 
(1/8)*log(-(1/2*I + 1/2)*sqrt(2)*a^5*c^4*(-1/(a^13*c^11))^(3/8) + sqrt(x)) 
 + 30*(a*c^3*x^8 + 2*a^2*c^2*x^4 + a^3*c)*(-1/(a^13*c^11))^(1/8)*log(a^5*c 
^4*(-1/(a^13*c^11))^(3/8) + sqrt(x)) + 30*(-I*a*c^3*x^8 - 2*I*a^2*c^2*x^4 
- I*a^3*c)*(-1/(a^13*c^11))^(1/8)*log(I*a^5*c^4*(-1/(a^13*c^11))^(3/8) + s 
qrt(x)) + 30*(I*a*c^3*x^8 + 2*I*a^2*c^2*x^4 + I*a^3*c)*(-1/(a^13*c^11))^(1 
/8)*log(-I*a^5*c^4*(-1/(a^13*c^11))^(3/8) + sqrt(x)) - 30*(a*c^3*x^8 + 2*a 
^2*c^2*x^4 + a^3*c)*(-1/(a^13*c^11))^(1/8)*log(-a^5*c^4*(-1/(a^13*c^11))^( 
3/8) + sqrt(x)) - 16*(3*c*x^5 - 5*a*x)*sqrt(x))/(a*c^3*x^8 + 2*a^2*c^2*x^4 
 + a^3*c)
 

Sympy [F(-1)]

Timed out. \[ \int \frac {x^{9/2}}{\left (a+c x^4\right )^3} \, dx=\text {Timed out} \] Input:

integrate(x**(9/2)/(c*x**4+a)**3,x)
 

Output:

Timed out
 

Maxima [F]

\[ \int \frac {x^{9/2}}{\left (a+c x^4\right )^3} \, dx=\int { \frac {x^{\frac {9}{2}}}{{\left (c x^{4} + a\right )}^{3}} \,d x } \] Input:

integrate(x^(9/2)/(c*x^4+a)^3,x, algorithm="maxima")
                                                                                    
                                                                                    
 

Output:

1/64*(3*c*x^(11/2) - 5*a*x^(3/2))/(a*c^3*x^8 + 2*a^2*c^2*x^4 + a^3*c) + 15 
*integrate(1/128*sqrt(x)/(a*c^2*x^4 + a^2*c), x)
 

Giac [B] (verification not implemented)

Leaf count of result is larger than twice the leaf count of optimal. 499 vs. \(2 (188) = 376\).

Time = 0.27 (sec) , antiderivative size = 499, normalized size of antiderivative = 1.83 \[ \int \frac {x^{9/2}}{\left (a+c x^4\right )^3} \, dx =\text {Too large to display} \] Input:

integrate(x^(9/2)/(c*x^4+a)^3,x, algorithm="giac")
 

Output:

-15/256*(a/c)^(3/8)*arctan((sqrt(-sqrt(2) + 2)*(a/c)^(1/8) + 2*sqrt(x))/(s 
qrt(sqrt(2) + 2)*(a/c)^(1/8)))/(a^2*c*sqrt(2*sqrt(2) + 4)) - 15/256*(a/c)^ 
(3/8)*arctan(-(sqrt(-sqrt(2) + 2)*(a/c)^(1/8) - 2*sqrt(x))/(sqrt(sqrt(2) + 
 2)*(a/c)^(1/8)))/(a^2*c*sqrt(2*sqrt(2) + 4)) + 15/256*(a/c)^(3/8)*arctan( 
(sqrt(sqrt(2) + 2)*(a/c)^(1/8) + 2*sqrt(x))/(sqrt(-sqrt(2) + 2)*(a/c)^(1/8 
)))/(a^2*c*sqrt(-2*sqrt(2) + 4)) + 15/256*(a/c)^(3/8)*arctan(-(sqrt(sqrt(2 
) + 2)*(a/c)^(1/8) - 2*sqrt(x))/(sqrt(-sqrt(2) + 2)*(a/c)^(1/8)))/(a^2*c*s 
qrt(-2*sqrt(2) + 4)) + 15/512*(a/c)^(3/8)*log(sqrt(x)*sqrt(sqrt(2) + 2)*(a 
/c)^(1/8) + x + (a/c)^(1/4))/(a^2*c*sqrt(2*sqrt(2) + 4)) - 15/512*(a/c)^(3 
/8)*log(-sqrt(x)*sqrt(sqrt(2) + 2)*(a/c)^(1/8) + x + (a/c)^(1/4))/(a^2*c*s 
qrt(2*sqrt(2) + 4)) - 15/512*(a/c)^(3/8)*log(sqrt(x)*sqrt(-sqrt(2) + 2)*(a 
/c)^(1/8) + x + (a/c)^(1/4))/(a^2*c*sqrt(-2*sqrt(2) + 4)) + 15/512*(a/c)^( 
3/8)*log(-sqrt(x)*sqrt(-sqrt(2) + 2)*(a/c)^(1/8) + x + (a/c)^(1/4))/(a^2*c 
*sqrt(-2*sqrt(2) + 4)) + 1/64*(3*c*x^(11/2) - 5*a*x^(3/2))/((c*x^4 + a)^2* 
a*c)
 

Mupad [B] (verification not implemented)

Time = 0.25 (sec) , antiderivative size = 156, normalized size of antiderivative = 0.57 \[ \int \frac {x^{9/2}}{\left (a+c x^4\right )^3} \, dx=\frac {\frac {3\,x^{11/2}}{64\,a}-\frac {5\,x^{3/2}}{64\,c}}{a^2+2\,a\,c\,x^4+c^2\,x^8}-\frac {15\,\mathrm {atan}\left (\frac {c^{1/8}\,\sqrt {x}}{{\left (-a\right )}^{1/8}}\right )}{256\,{\left (-a\right )}^{13/8}\,c^{11/8}}-\frac {\mathrm {atan}\left (\frac {c^{1/8}\,\sqrt {x}\,1{}\mathrm {i}}{{\left (-a\right )}^{1/8}}\right )\,15{}\mathrm {i}}{256\,{\left (-a\right )}^{13/8}\,c^{11/8}}+\frac {\sqrt {2}\,\mathrm {atan}\left (\frac {\sqrt {2}\,c^{1/8}\,\sqrt {x}\,\left (\frac {1}{2}-\frac {1}{2}{}\mathrm {i}\right )}{{\left (-a\right )}^{1/8}}\right )\,\left (\frac {15}{512}-\frac {15}{512}{}\mathrm {i}\right )}{{\left (-a\right )}^{13/8}\,c^{11/8}}+\frac {\sqrt {2}\,\mathrm {atan}\left (\frac {\sqrt {2}\,c^{1/8}\,\sqrt {x}\,\left (\frac {1}{2}+\frac {1}{2}{}\mathrm {i}\right )}{{\left (-a\right )}^{1/8}}\right )\,\left (\frac {15}{512}+\frac {15}{512}{}\mathrm {i}\right )}{{\left (-a\right )}^{13/8}\,c^{11/8}} \] Input:

int(x^(9/2)/(a + c*x^4)^3,x)
 

Output:

((3*x^(11/2))/(64*a) - (5*x^(3/2))/(64*c))/(a^2 + c^2*x^8 + 2*a*c*x^4) - ( 
15*atan((c^(1/8)*x^(1/2))/(-a)^(1/8)))/(256*(-a)^(13/8)*c^(11/8)) - (atan( 
(c^(1/8)*x^(1/2)*1i)/(-a)^(1/8))*15i)/(256*(-a)^(13/8)*c^(11/8)) + (2^(1/2 
)*atan((2^(1/2)*c^(1/8)*x^(1/2)*(1/2 - 1i/2))/(-a)^(1/8))*(15/512 - 15i/51 
2))/((-a)^(13/8)*c^(11/8)) + (2^(1/2)*atan((2^(1/2)*c^(1/8)*x^(1/2)*(1/2 + 
 1i/2))/(-a)^(1/8))*(15/512 + 15i/512))/((-a)^(13/8)*c^(11/8))
 

Reduce [B] (verification not implemented)

Time = 0.32 (sec) , antiderivative size = 2365, normalized size of antiderivative = 8.66 \[ \int \frac {x^{9/2}}{\left (a+c x^4\right )^3} \, dx =\text {Too large to display} \] Input:

int(x^(9/2)/(c*x^4+a)^3,x)
 

Output:

(30*c**(5/8)*a**(3/8)*sqrt(sqrt(2) + 2)*sqrt(2)*atan((c**(1/8)*a**(1/8)*sq 
rt( - sqrt(2) + 2) - 2*sqrt(x)*c**(1/4))/(c**(1/8)*a**(1/8)*sqrt(sqrt(2) + 
 2)))*a**2 + 60*c**(5/8)*a**(3/8)*sqrt(sqrt(2) + 2)*sqrt(2)*atan((c**(1/8) 
*a**(1/8)*sqrt( - sqrt(2) + 2) - 2*sqrt(x)*c**(1/4))/(c**(1/8)*a**(1/8)*sq 
rt(sqrt(2) + 2)))*a*c*x**4 + 30*c**(5/8)*a**(3/8)*sqrt(sqrt(2) + 2)*sqrt(2 
)*atan((c**(1/8)*a**(1/8)*sqrt( - sqrt(2) + 2) - 2*sqrt(x)*c**(1/4))/(c**( 
1/8)*a**(1/8)*sqrt(sqrt(2) + 2)))*c**2*x**8 - 30*c**(5/8)*a**(3/8)*sqrt(sq 
rt(2) + 2)*atan((c**(1/8)*a**(1/8)*sqrt( - sqrt(2) + 2) - 2*sqrt(x)*c**(1/ 
4))/(c**(1/8)*a**(1/8)*sqrt(sqrt(2) + 2)))*a**2 - 60*c**(5/8)*a**(3/8)*sqr 
t(sqrt(2) + 2)*atan((c**(1/8)*a**(1/8)*sqrt( - sqrt(2) + 2) - 2*sqrt(x)*c* 
*(1/4))/(c**(1/8)*a**(1/8)*sqrt(sqrt(2) + 2)))*a*c*x**4 - 30*c**(5/8)*a**( 
3/8)*sqrt(sqrt(2) + 2)*atan((c**(1/8)*a**(1/8)*sqrt( - sqrt(2) + 2) - 2*sq 
rt(x)*c**(1/4))/(c**(1/8)*a**(1/8)*sqrt(sqrt(2) + 2)))*c**2*x**8 - 30*c**( 
5/8)*a**(3/8)*sqrt(sqrt(2) + 2)*sqrt(2)*atan((c**(1/8)*a**(1/8)*sqrt( - sq 
rt(2) + 2) + 2*sqrt(x)*c**(1/4))/(c**(1/8)*a**(1/8)*sqrt(sqrt(2) + 2)))*a* 
*2 - 60*c**(5/8)*a**(3/8)*sqrt(sqrt(2) + 2)*sqrt(2)*atan((c**(1/8)*a**(1/8 
)*sqrt( - sqrt(2) + 2) + 2*sqrt(x)*c**(1/4))/(c**(1/8)*a**(1/8)*sqrt(sqrt( 
2) + 2)))*a*c*x**4 - 30*c**(5/8)*a**(3/8)*sqrt(sqrt(2) + 2)*sqrt(2)*atan(( 
c**(1/8)*a**(1/8)*sqrt( - sqrt(2) + 2) + 2*sqrt(x)*c**(1/4))/(c**(1/8)*a** 
(1/8)*sqrt(sqrt(2) + 2)))*c**2*x**8 + 30*c**(5/8)*a**(3/8)*sqrt(sqrt(2)...