\(\int \frac {x^3}{\sqrt {a+c x^2} (d+e x+f x^2)} \, dx\) [31]

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

Optimal result

Integrand size = 27, antiderivative size = 380 \[ \int \frac {x^3}{\sqrt {a+c x^2} \left (d+e x+f x^2\right )} \, dx=\frac {\sqrt {a+c x^2}}{c f}-\frac {e \text {arctanh}\left (\frac {\sqrt {c} x}{\sqrt {a+c x^2}}\right )}{\sqrt {c} f^2}-\frac {\left (2 d e f-\left (e^2-d f\right ) \left (e-\sqrt {e^2-4 d f}\right )\right ) \text {arctanh}\left (\frac {2 a f-c \left (e-\sqrt {e^2-4 d f}\right ) x}{\sqrt {2} \sqrt {2 a f^2+c \left (e^2-2 d f-e \sqrt {e^2-4 d f}\right )} \sqrt {a+c x^2}}\right )}{\sqrt {2} f^2 \sqrt {e^2-4 d f} \sqrt {2 a f^2+c \left (e^2-2 d f-e \sqrt {e^2-4 d f}\right )}}+\frac {\left (2 d e f-\left (e^2-d f\right ) \left (e+\sqrt {e^2-4 d f}\right )\right ) \text {arctanh}\left (\frac {2 a f-c \left (e+\sqrt {e^2-4 d f}\right ) x}{\sqrt {2} \sqrt {2 a f^2+c \left (e^2-2 d f+e \sqrt {e^2-4 d f}\right )} \sqrt {a+c x^2}}\right )}{\sqrt {2} f^2 \sqrt {e^2-4 d f} \sqrt {2 a f^2+c \left (e^2-2 d f+e \sqrt {e^2-4 d f}\right )}} \] Output:

(c*x^2+a)^(1/2)/c/f-e*arctanh(c^(1/2)*x/(c*x^2+a)^(1/2))/c^(1/2)/f^2-1/2*( 
2*d*e*f-(-d*f+e^2)*(e-(-4*d*f+e^2)^(1/2)))*arctanh(1/2*(2*a*f-c*(e-(-4*d*f 
+e^2)^(1/2))*x)*2^(1/2)/(2*a*f^2+c*(e^2-2*d*f-e*(-4*d*f+e^2)^(1/2)))^(1/2) 
/(c*x^2+a)^(1/2))*2^(1/2)/f^2/(-4*d*f+e^2)^(1/2)/(2*a*f^2+c*(e^2-2*d*f-e*( 
-4*d*f+e^2)^(1/2)))^(1/2)+1/2*(2*d*e*f-(-d*f+e^2)*(e+(-4*d*f+e^2)^(1/2)))* 
arctanh(1/2*(2*a*f-c*(e+(-4*d*f+e^2)^(1/2))*x)*2^(1/2)/(2*a*f^2+c*(e^2-2*d 
*f+e*(-4*d*f+e^2)^(1/2)))^(1/2)/(c*x^2+a)^(1/2))*2^(1/2)/f^2/(-4*d*f+e^2)^ 
(1/2)/(2*a*f^2+c*(e^2-2*d*f+e*(-4*d*f+e^2)^(1/2)))^(1/2)
 

Mathematica [C] (verified)

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

Time = 0.55 (sec) , antiderivative size = 312, normalized size of antiderivative = 0.82 \[ \int \frac {x^3}{\sqrt {a+c x^2} \left (d+e x+f x^2\right )} \, dx=\frac {f \sqrt {a+c x^2}+\sqrt {c} e \log \left (-\sqrt {c} x+\sqrt {a+c x^2}\right )-c \text {RootSum}\left [a^2 f+2 a \sqrt {c} e \text {$\#$1}+4 c d \text {$\#$1}^2-2 a f \text {$\#$1}^2-2 \sqrt {c} e \text {$\#$1}^3+f \text {$\#$1}^4\&,\frac {a e^2 \log \left (-\sqrt {c} x+\sqrt {a+c x^2}-\text {$\#$1}\right )-a d f \log \left (-\sqrt {c} x+\sqrt {a+c x^2}-\text {$\#$1}\right )+2 \sqrt {c} d e \log \left (-\sqrt {c} x+\sqrt {a+c x^2}-\text {$\#$1}\right ) \text {$\#$1}-e^2 \log \left (-\sqrt {c} x+\sqrt {a+c x^2}-\text {$\#$1}\right ) \text {$\#$1}^2+d f \log \left (-\sqrt {c} x+\sqrt {a+c x^2}-\text {$\#$1}\right ) \text {$\#$1}^2}{a \sqrt {c} e+4 c d \text {$\#$1}-2 a f \text {$\#$1}-3 \sqrt {c} e \text {$\#$1}^2+2 f \text {$\#$1}^3}\&\right ]}{c f^2} \] Input:

Integrate[x^3/(Sqrt[a + c*x^2]*(d + e*x + f*x^2)),x]
 

Output:

(f*Sqrt[a + c*x^2] + Sqrt[c]*e*Log[-(Sqrt[c]*x) + Sqrt[a + c*x^2]] - c*Roo 
tSum[a^2*f + 2*a*Sqrt[c]*e*#1 + 4*c*d*#1^2 - 2*a*f*#1^2 - 2*Sqrt[c]*e*#1^3 
 + f*#1^4 & , (a*e^2*Log[-(Sqrt[c]*x) + Sqrt[a + c*x^2] - #1] - a*d*f*Log[ 
-(Sqrt[c]*x) + Sqrt[a + c*x^2] - #1] + 2*Sqrt[c]*d*e*Log[-(Sqrt[c]*x) + Sq 
rt[a + c*x^2] - #1]*#1 - e^2*Log[-(Sqrt[c]*x) + Sqrt[a + c*x^2] - #1]*#1^2 
 + d*f*Log[-(Sqrt[c]*x) + Sqrt[a + c*x^2] - #1]*#1^2)/(a*Sqrt[c]*e + 4*c*d 
*#1 - 2*a*f*#1 - 3*Sqrt[c]*e*#1^2 + 2*f*#1^3) & ])/(c*f^2)
 

Rubi [A] (verified)

Time = 1.06 (sec) , antiderivative size = 380, normalized size of antiderivative = 1.00, number of steps used = 2, number of rules used = 2, \(\frac {\text {number of rules}}{\text {integrand size}}\) = 0.074, Rules used = {7279, 2009}

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^3}{\sqrt {a+c x^2} \left (d+e x+f x^2\right )} \, dx\)

\(\Big \downarrow \) 7279

\(\displaystyle \int \left (\frac {x \left (e^2-d f\right )+d e}{f^2 \sqrt {a+c x^2} \left (d+e x+f x^2\right )}-\frac {e}{f^2 \sqrt {a+c x^2}}+\frac {x}{f \sqrt {a+c x^2}}\right )dx\)

\(\Big \downarrow \) 2009

\(\displaystyle -\frac {\left (2 d e f-\left (e^2-d f\right ) \left (e-\sqrt {e^2-4 d f}\right )\right ) \text {arctanh}\left (\frac {2 a f-c x \left (e-\sqrt {e^2-4 d f}\right )}{\sqrt {2} \sqrt {a+c x^2} \sqrt {2 a f^2+c \left (-e \sqrt {e^2-4 d f}-2 d f+e^2\right )}}\right )}{\sqrt {2} f^2 \sqrt {e^2-4 d f} \sqrt {2 a f^2+c \left (-e \sqrt {e^2-4 d f}-2 d f+e^2\right )}}+\frac {\left (2 d e f-\left (e^2-d f\right ) \left (\sqrt {e^2-4 d f}+e\right )\right ) \text {arctanh}\left (\frac {2 a f-c x \left (\sqrt {e^2-4 d f}+e\right )}{\sqrt {2} \sqrt {a+c x^2} \sqrt {2 a f^2+c \left (e \sqrt {e^2-4 d f}-2 d f+e^2\right )}}\right )}{\sqrt {2} f^2 \sqrt {e^2-4 d f} \sqrt {2 a f^2+c \left (e \sqrt {e^2-4 d f}-2 d f+e^2\right )}}-\frac {e \text {arctanh}\left (\frac {\sqrt {c} x}{\sqrt {a+c x^2}}\right )}{\sqrt {c} f^2}+\frac {\sqrt {a+c x^2}}{c f}\)

Input:

Int[x^3/(Sqrt[a + c*x^2]*(d + e*x + f*x^2)),x]
 

Output:

Sqrt[a + c*x^2]/(c*f) - (e*ArcTanh[(Sqrt[c]*x)/Sqrt[a + c*x^2]])/(Sqrt[c]* 
f^2) - ((2*d*e*f - (e^2 - d*f)*(e - Sqrt[e^2 - 4*d*f]))*ArcTanh[(2*a*f - c 
*(e - Sqrt[e^2 - 4*d*f])*x)/(Sqrt[2]*Sqrt[2*a*f^2 + c*(e^2 - 2*d*f - e*Sqr 
t[e^2 - 4*d*f])]*Sqrt[a + c*x^2])])/(Sqrt[2]*f^2*Sqrt[e^2 - 4*d*f]*Sqrt[2* 
a*f^2 + c*(e^2 - 2*d*f - e*Sqrt[e^2 - 4*d*f])]) + ((2*d*e*f - (e^2 - d*f)* 
(e + Sqrt[e^2 - 4*d*f]))*ArcTanh[(2*a*f - c*(e + Sqrt[e^2 - 4*d*f])*x)/(Sq 
rt[2]*Sqrt[2*a*f^2 + c*(e^2 - 2*d*f + e*Sqrt[e^2 - 4*d*f])]*Sqrt[a + c*x^2 
])])/(Sqrt[2]*f^2*Sqrt[e^2 - 4*d*f]*Sqrt[2*a*f^2 + c*(e^2 - 2*d*f + e*Sqrt 
[e^2 - 4*d*f])])
 

Defintions of rubi rules used

rule 2009
Int[u_, x_Symbol] :> Simp[IntSum[u, x], x] /; SumQ[u]
 

rule 7279
Int[(u_)/((a_.) + (b_.)*(x_)^(n_.) + (c_.)*(x_)^(n2_.)), x_Symbol] :> With[ 
{v = RationalFunctionExpand[u/(a + b*x^n + c*x^(2*n)), x]}, Int[v, x] /; Su 
mQ[v]] /; FreeQ[{a, b, c}, x] && EqQ[n2, 2*n] && IGtQ[n, 0]
 
Maple [B] (verified)

Leaf count of result is larger than twice the leaf count of optimal. \(712\) vs. \(2(335)=670\).

Time = 2.24 (sec) , antiderivative size = 713, normalized size of antiderivative = 1.88

method result size
default \(\frac {\sqrt {c \,x^{2}+a}}{c f}+\frac {\left (d f \sqrt {-4 d f +e^{2}}-e^{2} \sqrt {-4 d f +e^{2}}+3 d e f -e^{3}\right ) \sqrt {2}\, \ln \left (\frac {\frac {\sqrt {-4 d f +e^{2}}\, c e +2 a \,f^{2}-2 d f c +c \,e^{2}}{f^{2}}-\frac {c \left (e +\sqrt {-4 d f +e^{2}}\right ) \left (x +\frac {e +\sqrt {-4 d f +e^{2}}}{2 f}\right )}{f}+\frac {\sqrt {2}\, \sqrt {\frac {\sqrt {-4 d f +e^{2}}\, c e +2 a \,f^{2}-2 d f c +c \,e^{2}}{f^{2}}}\, \sqrt {4 c {\left (x +\frac {e +\sqrt {-4 d f +e^{2}}}{2 f}\right )}^{2}-\frac {4 c \left (e +\sqrt {-4 d f +e^{2}}\right ) \left (x +\frac {e +\sqrt {-4 d f +e^{2}}}{2 f}\right )}{f}+\frac {2 \sqrt {-4 d f +e^{2}}\, c e +4 a \,f^{2}-4 d f c +2 c \,e^{2}}{f^{2}}}}{2}}{x +\frac {e +\sqrt {-4 d f +e^{2}}}{2 f}}\right )}{2 f^{3} \sqrt {-4 d f +e^{2}}\, \sqrt {\frac {\sqrt {-4 d f +e^{2}}\, c e +2 a \,f^{2}-2 d f c +c \,e^{2}}{f^{2}}}}+\frac {\left (d f \sqrt {-4 d f +e^{2}}-e^{2} \sqrt {-4 d f +e^{2}}-3 d e f +e^{3}\right ) \sqrt {2}\, \ln \left (\frac {\frac {-\sqrt {-4 d f +e^{2}}\, c e +2 a \,f^{2}-2 d f c +c \,e^{2}}{f^{2}}-\frac {c \left (e -\sqrt {-4 d f +e^{2}}\right ) \left (x -\frac {-e +\sqrt {-4 d f +e^{2}}}{2 f}\right )}{f}+\frac {\sqrt {2}\, \sqrt {\frac {-\sqrt {-4 d f +e^{2}}\, c e +2 a \,f^{2}-2 d f c +c \,e^{2}}{f^{2}}}\, \sqrt {4 c {\left (x -\frac {-e +\sqrt {-4 d f +e^{2}}}{2 f}\right )}^{2}-\frac {4 c \left (e -\sqrt {-4 d f +e^{2}}\right ) \left (x -\frac {-e +\sqrt {-4 d f +e^{2}}}{2 f}\right )}{f}+\frac {-2 \sqrt {-4 d f +e^{2}}\, c e +4 a \,f^{2}-4 d f c +2 c \,e^{2}}{f^{2}}}}{2}}{x -\frac {-e +\sqrt {-4 d f +e^{2}}}{2 f}}\right )}{2 f^{3} \sqrt {-4 d f +e^{2}}\, \sqrt {\frac {-\sqrt {-4 d f +e^{2}}\, c e +2 a \,f^{2}-2 d f c +c \,e^{2}}{f^{2}}}}-\frac {e \ln \left (\sqrt {c}\, x +\sqrt {c \,x^{2}+a}\right )}{f^{2} \sqrt {c}}\) \(713\)
risch \(\frac {\sqrt {c \,x^{2}+a}}{c f}-\frac {\frac {e \ln \left (\sqrt {c}\, x +\sqrt {c \,x^{2}+a}\right )}{f \sqrt {c}}-\frac {\left (d f \sqrt {-4 d f +e^{2}}-e^{2} \sqrt {-4 d f +e^{2}}-3 d e f +e^{3}\right ) \sqrt {2}\, \ln \left (\frac {\frac {-\sqrt {-4 d f +e^{2}}\, c e +2 a \,f^{2}-2 d f c +c \,e^{2}}{f^{2}}-\frac {c \left (e -\sqrt {-4 d f +e^{2}}\right ) \left (x -\frac {-e +\sqrt {-4 d f +e^{2}}}{2 f}\right )}{f}+\frac {\sqrt {2}\, \sqrt {\frac {-\sqrt {-4 d f +e^{2}}\, c e +2 a \,f^{2}-2 d f c +c \,e^{2}}{f^{2}}}\, \sqrt {4 c {\left (x -\frac {-e +\sqrt {-4 d f +e^{2}}}{2 f}\right )}^{2}-\frac {4 c \left (e -\sqrt {-4 d f +e^{2}}\right ) \left (x -\frac {-e +\sqrt {-4 d f +e^{2}}}{2 f}\right )}{f}+\frac {-2 \sqrt {-4 d f +e^{2}}\, c e +4 a \,f^{2}-4 d f c +2 c \,e^{2}}{f^{2}}}}{2}}{x -\frac {-e +\sqrt {-4 d f +e^{2}}}{2 f}}\right )}{2 f^{2} \sqrt {-4 d f +e^{2}}\, \sqrt {\frac {-\sqrt {-4 d f +e^{2}}\, c e +2 a \,f^{2}-2 d f c +c \,e^{2}}{f^{2}}}}-\frac {\left (d f \sqrt {-4 d f +e^{2}}-e^{2} \sqrt {-4 d f +e^{2}}+3 d e f -e^{3}\right ) \sqrt {2}\, \ln \left (\frac {\frac {\sqrt {-4 d f +e^{2}}\, c e +2 a \,f^{2}-2 d f c +c \,e^{2}}{f^{2}}-\frac {c \left (e +\sqrt {-4 d f +e^{2}}\right ) \left (x +\frac {e +\sqrt {-4 d f +e^{2}}}{2 f}\right )}{f}+\frac {\sqrt {2}\, \sqrt {\frac {\sqrt {-4 d f +e^{2}}\, c e +2 a \,f^{2}-2 d f c +c \,e^{2}}{f^{2}}}\, \sqrt {4 c {\left (x +\frac {e +\sqrt {-4 d f +e^{2}}}{2 f}\right )}^{2}-\frac {4 c \left (e +\sqrt {-4 d f +e^{2}}\right ) \left (x +\frac {e +\sqrt {-4 d f +e^{2}}}{2 f}\right )}{f}+\frac {2 \sqrt {-4 d f +e^{2}}\, c e +4 a \,f^{2}-4 d f c +2 c \,e^{2}}{f^{2}}}}{2}}{x +\frac {e +\sqrt {-4 d f +e^{2}}}{2 f}}\right )}{2 f^{2} \sqrt {-4 d f +e^{2}}\, \sqrt {\frac {\sqrt {-4 d f +e^{2}}\, c e +2 a \,f^{2}-2 d f c +c \,e^{2}}{f^{2}}}}}{f}\) \(718\)

Input:

int(x^3/(c*x^2+a)^(1/2)/(f*x^2+e*x+d),x,method=_RETURNVERBOSE)
 

Output:

(c*x^2+a)^(1/2)/c/f+1/2/f^3*(d*f*(-4*d*f+e^2)^(1/2)-e^2*(-4*d*f+e^2)^(1/2) 
+3*d*e*f-e^3)/(-4*d*f+e^2)^(1/2)*2^(1/2)/(((-4*d*f+e^2)^(1/2)*c*e+2*a*f^2- 
2*d*f*c+c*e^2)/f^2)^(1/2)*ln((((-4*d*f+e^2)^(1/2)*c*e+2*a*f^2-2*d*f*c+c*e^ 
2)/f^2-c*(e+(-4*d*f+e^2)^(1/2))/f*(x+1/2*(e+(-4*d*f+e^2)^(1/2))/f)+1/2*2^( 
1/2)*(((-4*d*f+e^2)^(1/2)*c*e+2*a*f^2-2*d*f*c+c*e^2)/f^2)^(1/2)*(4*c*(x+1/ 
2*(e+(-4*d*f+e^2)^(1/2))/f)^2-4*c*(e+(-4*d*f+e^2)^(1/2))/f*(x+1/2*(e+(-4*d 
*f+e^2)^(1/2))/f)+2*((-4*d*f+e^2)^(1/2)*c*e+2*a*f^2-2*d*f*c+c*e^2)/f^2)^(1 
/2))/(x+1/2*(e+(-4*d*f+e^2)^(1/2))/f))+1/2/f^3*(d*f*(-4*d*f+e^2)^(1/2)-e^2 
*(-4*d*f+e^2)^(1/2)-3*d*e*f+e^3)/(-4*d*f+e^2)^(1/2)*2^(1/2)/((-(-4*d*f+e^2 
)^(1/2)*c*e+2*a*f^2-2*d*f*c+c*e^2)/f^2)^(1/2)*ln(((-(-4*d*f+e^2)^(1/2)*c*e 
+2*a*f^2-2*d*f*c+c*e^2)/f^2-c*(e-(-4*d*f+e^2)^(1/2))/f*(x-1/2/f*(-e+(-4*d* 
f+e^2)^(1/2)))+1/2*2^(1/2)*((-(-4*d*f+e^2)^(1/2)*c*e+2*a*f^2-2*d*f*c+c*e^2 
)/f^2)^(1/2)*(4*c*(x-1/2/f*(-e+(-4*d*f+e^2)^(1/2)))^2-4*c*(e-(-4*d*f+e^2)^ 
(1/2))/f*(x-1/2/f*(-e+(-4*d*f+e^2)^(1/2)))+2*(-(-4*d*f+e^2)^(1/2)*c*e+2*a* 
f^2-2*d*f*c+c*e^2)/f^2)^(1/2))/(x-1/2/f*(-e+(-4*d*f+e^2)^(1/2))))-e/f^2*ln 
(c^(1/2)*x+(c*x^2+a)^(1/2))/c^(1/2)
 

Fricas [F(-1)]

Timed out. \[ \int \frac {x^3}{\sqrt {a+c x^2} \left (d+e x+f x^2\right )} \, dx=\text {Timed out} \] Input:

integrate(x^3/(c*x^2+a)^(1/2)/(f*x^2+e*x+d),x, algorithm="fricas")
 

Output:

Timed out
 

Sympy [F]

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

integrate(x**3/(c*x**2+a)**(1/2)/(f*x**2+e*x+d),x)
 

Output:

Integral(x**3/(sqrt(a + c*x**2)*(d + e*x + f*x**2)), x)
 

Maxima [F(-2)]

Exception generated. \[ \int \frac {x^3}{\sqrt {a+c x^2} \left (d+e x+f x^2\right )} \, dx=\text {Exception raised: ValueError} \] Input:

integrate(x^3/(c*x^2+a)^(1/2)/(f*x^2+e*x+d),x, algorithm="maxima")
 

Output:

Exception raised: ValueError >> Computation failed since Maxima requested 
additional constraints; using the 'assume' command before evaluation *may* 
 help (example of legal syntax is 'assume(4*d*f-e^2>0)', see `assume?` for 
 more deta
 

Giac [F(-2)]

Exception generated. \[ \int \frac {x^3}{\sqrt {a+c x^2} \left (d+e x+f x^2\right )} \, dx=\text {Exception raised: TypeError} \] Input:

integrate(x^3/(c*x^2+a)^(1/2)/(f*x^2+e*x+d),x, algorithm="giac")
 

Output:

Exception raised: TypeError >> an error occurred running a Giac command:IN 
PUT:sage2:=int(sage0,sageVARx):;OUTPUT:index.cc index_m i_lex_is_greater E 
rror: Bad Argument Value
                                                                                    
                                                                                    
 

Mupad [F(-1)]

Timed out. \[ \int \frac {x^3}{\sqrt {a+c x^2} \left (d+e x+f x^2\right )} \, dx=\int \frac {x^3}{\sqrt {c\,x^2+a}\,\left (f\,x^2+e\,x+d\right )} \,d x \] Input:

int(x^3/((a + c*x^2)^(1/2)*(d + e*x + f*x^2)),x)
 

Output:

int(x^3/((a + c*x^2)^(1/2)*(d + e*x + f*x^2)), x)
 

Reduce [F]

\[ \int \frac {x^3}{\sqrt {a+c x^2} \left (d+e x+f x^2\right )} \, dx=\frac {4 \sqrt {c \,x^{2}+a}\, d \,e^{2} f -2 \sqrt {c \,x^{2}+a}\, e^{4}+\sqrt {c}\, \mathrm {log}\left (\sqrt {c \,x^{2}+a}-\sqrt {c}\, x \right ) d^{2} e f -\sqrt {c}\, \mathrm {log}\left (\sqrt {c \,x^{2}+a}+\sqrt {c}\, x \right ) d^{2} e f +2 \left (\int \frac {\sqrt {c \,x^{2}+a}}{c f \,x^{4}+c e \,x^{3}+a f \,x^{2}+c d \,x^{2}+a e x +a d}d x \right ) c \,d^{3} e f +2 \left (\int \frac {\sqrt {c \,x^{2}+a}\, x^{3}}{c f \,x^{4}+c e \,x^{3}+a f \,x^{2}+c d \,x^{2}+a e x +a d}d x \right ) c \,d^{2} f^{3}-4 \left (\int \frac {\sqrt {c \,x^{2}+a}\, x^{3}}{c f \,x^{4}+c e \,x^{3}+a f \,x^{2}+c d \,x^{2}+a e x +a d}d x \right ) c d \,e^{2} f^{2}+2 \left (\int \frac {\sqrt {c \,x^{2}+a}\, x^{3}}{c f \,x^{4}+c e \,x^{3}+a f \,x^{2}+c d \,x^{2}+a e x +a d}d x \right ) c \,e^{4} f +2 \left (\int \frac {\sqrt {c \,x^{2}+a}\, x^{2}}{c f \,x^{4}+c e \,x^{3}+a f \,x^{2}+c d \,x^{2}+a e x +a d}d x \right ) c \,d^{2} e \,f^{2}-4 \left (\int \frac {\sqrt {c \,x^{2}+a}\, x^{2}}{c f \,x^{4}+c e \,x^{3}+a f \,x^{2}+c d \,x^{2}+a e x +a d}d x \right ) c d \,e^{3} f +2 \left (\int \frac {\sqrt {c \,x^{2}+a}\, x^{2}}{c f \,x^{4}+c e \,x^{3}+a f \,x^{2}+c d \,x^{2}+a e x +a d}d x \right ) c \,e^{5}-2 \left (\int \frac {\sqrt {c \,x^{2}+a}\, x}{c f \,x^{4}+c e \,x^{3}+a f \,x^{2}+c d \,x^{2}+a e x +a d}d x \right ) c \,d^{2} e^{2} f +2 \left (\int \frac {\sqrt {c \,x^{2}+a}\, x}{c f \,x^{4}+c e \,x^{3}+a f \,x^{2}+c d \,x^{2}+a e x +a d}d x \right ) c d \,e^{4}}{2 c \,d^{2} f^{3}} \] Input:

int(x^3/(c*x^2+a)^(1/2)/(f*x^2+e*x+d),x)
 

Output:

(4*sqrt(a + c*x**2)*d*e**2*f - 2*sqrt(a + c*x**2)*e**4 + sqrt(c)*log(sqrt( 
a + c*x**2) - sqrt(c)*x)*d**2*e*f - sqrt(c)*log(sqrt(a + c*x**2) + sqrt(c) 
*x)*d**2*e*f + 2*int(sqrt(a + c*x**2)/(a*d + a*e*x + a*f*x**2 + c*d*x**2 + 
 c*e*x**3 + c*f*x**4),x)*c*d**3*e*f + 2*int((sqrt(a + c*x**2)*x**3)/(a*d + 
 a*e*x + a*f*x**2 + c*d*x**2 + c*e*x**3 + c*f*x**4),x)*c*d**2*f**3 - 4*int 
((sqrt(a + c*x**2)*x**3)/(a*d + a*e*x + a*f*x**2 + c*d*x**2 + c*e*x**3 + c 
*f*x**4),x)*c*d*e**2*f**2 + 2*int((sqrt(a + c*x**2)*x**3)/(a*d + a*e*x + a 
*f*x**2 + c*d*x**2 + c*e*x**3 + c*f*x**4),x)*c*e**4*f + 2*int((sqrt(a + c* 
x**2)*x**2)/(a*d + a*e*x + a*f*x**2 + c*d*x**2 + c*e*x**3 + c*f*x**4),x)*c 
*d**2*e*f**2 - 4*int((sqrt(a + c*x**2)*x**2)/(a*d + a*e*x + a*f*x**2 + c*d 
*x**2 + c*e*x**3 + c*f*x**4),x)*c*d*e**3*f + 2*int((sqrt(a + c*x**2)*x**2) 
/(a*d + a*e*x + a*f*x**2 + c*d*x**2 + c*e*x**3 + c*f*x**4),x)*c*e**5 - 2*i 
nt((sqrt(a + c*x**2)*x)/(a*d + a*e*x + a*f*x**2 + c*d*x**2 + c*e*x**3 + c* 
f*x**4),x)*c*d**2*e**2*f + 2*int((sqrt(a + c*x**2)*x)/(a*d + a*e*x + a*f*x 
**2 + c*d*x**2 + c*e*x**3 + c*f*x**4),x)*c*d*e**4)/(2*c*d**2*f**3)