\(\int \frac {x^4}{\sqrt {a-b x^{3/2}} \sqrt {a+b x^{3/2}}} \, dx\) [231]

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

Optimal result

Integrand size = 31, antiderivative size = 579 \[ \int \frac {x^4}{\sqrt {a-b x^{3/2}} \sqrt {a+b x^{3/2}}} \, dx=-\frac {2 x^2 \sqrt {a-b x^{3/2}} \sqrt {a+b x^{3/2}}}{7 b^2}+\frac {8 a^2 \left (a^2-b^2 x^3\right )}{7 b^{10/3} \left (\left (1+\sqrt {3}\right ) a^{2/3}-b^{2/3} x\right ) \sqrt {a-b x^{3/2}} \sqrt {a+b x^{3/2}}}-\frac {4 \sqrt [4]{3} \sqrt {2-\sqrt {3}} a^{8/3} \left (a^{2/3}-b^{2/3} x\right ) \sqrt {\frac {a^{4/3}+a^{2/3} b^{2/3} x+b^{4/3} x^2}{\left (\left (1+\sqrt {3}\right ) a^{2/3}-b^{2/3} x\right )^2}} E\left (\arcsin \left (\frac {\left (1-\sqrt {3}\right ) a^{2/3}-b^{2/3} x}{\left (1+\sqrt {3}\right ) a^{2/3}-b^{2/3} x}\right )|-7-4 \sqrt {3}\right )}{7 b^{10/3} \sqrt {\frac {a^{2/3} \left (a^{2/3}-b^{2/3} x\right )}{\left (\left (1+\sqrt {3}\right ) a^{2/3}-b^{2/3} x\right )^2}} \sqrt {a-b x^{3/2}} \sqrt {a+b x^{3/2}}}+\frac {8 \sqrt {2} a^{8/3} \left (a^{2/3}-b^{2/3} x\right ) \sqrt {\frac {a^{4/3}+a^{2/3} b^{2/3} x+b^{4/3} x^2}{\left (\left (1+\sqrt {3}\right ) a^{2/3}-b^{2/3} x\right )^2}} \operatorname {EllipticF}\left (\arcsin \left (\frac {\left (1-\sqrt {3}\right ) a^{2/3}-b^{2/3} x}{\left (1+\sqrt {3}\right ) a^{2/3}-b^{2/3} x}\right ),-7-4 \sqrt {3}\right )}{7 \sqrt [4]{3} b^{10/3} \sqrt {\frac {a^{2/3} \left (a^{2/3}-b^{2/3} x\right )}{\left (\left (1+\sqrt {3}\right ) a^{2/3}-b^{2/3} x\right )^2}} \sqrt {a-b x^{3/2}} \sqrt {a+b x^{3/2}}} \] Output:

-2/7*x^2*(a-b*x^(3/2))^(1/2)*(a+b*x^(3/2))^(1/2)/b^2+8/7*a^2*(-b^2*x^3+a^2 
)/b^(10/3)/((1+3^(1/2))*a^(2/3)-b^(2/3)*x)/(a-b*x^(3/2))^(1/2)/(a+b*x^(3/2 
))^(1/2)-4/7*3^(1/4)*(1/2*6^(1/2)-1/2*2^(1/2))*a^(8/3)*(a^(2/3)-b^(2/3)*x) 
*((a^(4/3)+a^(2/3)*b^(2/3)*x+b^(4/3)*x^2)/((1+3^(1/2))*a^(2/3)-b^(2/3)*x)^ 
2)^(1/2)*EllipticE(((1-3^(1/2))*a^(2/3)-b^(2/3)*x)/((1+3^(1/2))*a^(2/3)-b^ 
(2/3)*x),I*3^(1/2)+2*I)/b^(10/3)/(a^(2/3)*(a^(2/3)-b^(2/3)*x)/((1+3^(1/2)) 
*a^(2/3)-b^(2/3)*x)^2)^(1/2)/(a-b*x^(3/2))^(1/2)/(a+b*x^(3/2))^(1/2)+8/21* 
2^(1/2)*a^(8/3)*(a^(2/3)-b^(2/3)*x)*((a^(4/3)+a^(2/3)*b^(2/3)*x+b^(4/3)*x^ 
2)/((1+3^(1/2))*a^(2/3)-b^(2/3)*x)^2)^(1/2)*EllipticF(((1-3^(1/2))*a^(2/3) 
-b^(2/3)*x)/((1+3^(1/2))*a^(2/3)-b^(2/3)*x),I*3^(1/2)+2*I)*3^(3/4)/b^(10/3 
)/(a^(2/3)*(a^(2/3)-b^(2/3)*x)/((1+3^(1/2))*a^(2/3)-b^(2/3)*x)^2)^(1/2)/(a 
-b*x^(3/2))^(1/2)/(a+b*x^(3/2))^(1/2)
 

Mathematica [C] (verified)

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

Time = 6.54 (sec) , antiderivative size = 73, normalized size of antiderivative = 0.13 \[ \int \frac {x^4}{\sqrt {a-b x^{3/2}} \sqrt {a+b x^{3/2}}} \, dx=\frac {x^5 \sqrt {1-\frac {b^2 x^3}{a^2}} \, _2F_1\left (\frac {1}{2},\frac {5}{3};\frac {8}{3};\frac {b^2 x^3}{a^2}\right )}{5 \sqrt {a-b x^{3/2}} \sqrt {a+b x^{3/2}}} \] Input:

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

Output:

(x^5*Sqrt[1 - (b^2*x^3)/a^2]*HypergeometricPFQ[{1/2, 5/3}, {8/3}, (b^2*x^3 
)/a^2])/(5*Sqrt[a - b*x^(3/2)]*Sqrt[a + b*x^(3/2)])
 

Rubi [A] (warning: unable to verify)

Time = 1.03 (sec) , antiderivative size = 601, normalized size of antiderivative = 1.04, number of steps used = 5, number of rules used = 5, \(\frac {\text {number of rules}}{\text {integrand size}}\) = 0.161, Rules used = {845, 890, 832, 759, 2416}

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^4}{\sqrt {a-b x^{3/2}} \sqrt {a+b x^{3/2}}} \, dx\)

\(\Big \downarrow \) 845

\(\displaystyle \frac {4 a^2 \int \frac {x}{\sqrt {a-b x^{3/2}} \sqrt {b x^{3/2}+a}}dx}{7 b^2}-\frac {2 x^2 \sqrt {a-b x^{3/2}} \sqrt {a+b x^{3/2}}}{7 b^2}\)

\(\Big \downarrow \) 890

\(\displaystyle \frac {4 a^2 \sqrt {a^2-b^2 x^3} \int \frac {x}{\sqrt {a^2-b^2 x^3}}dx}{7 b^2 \sqrt {a-b x^{3/2}} \sqrt {a+b x^{3/2}}}-\frac {2 x^2 \sqrt {a-b x^{3/2}} \sqrt {a+b x^{3/2}}}{7 b^2}\)

\(\Big \downarrow \) 832

\(\displaystyle \frac {4 a^2 \sqrt {a^2-b^2 x^3} \left (\frac {\left (1-\sqrt {3}\right ) a^{2/3} \int \frac {1}{\sqrt {a^2-b^2 x^3}}dx}{b^{2/3}}-\frac {\int \frac {\left (1-\sqrt {3}\right ) a^{2/3}-b^{2/3} x}{\sqrt {a^2-b^2 x^3}}dx}{b^{2/3}}\right )}{7 b^2 \sqrt {a-b x^{3/2}} \sqrt {a+b x^{3/2}}}-\frac {2 x^2 \sqrt {a-b x^{3/2}} \sqrt {a+b x^{3/2}}}{7 b^2}\)

\(\Big \downarrow \) 759

\(\displaystyle \frac {4 a^2 \sqrt {a^2-b^2 x^3} \left (-\frac {\int \frac {\left (1-\sqrt {3}\right ) a^{2/3}-b^{2/3} x}{\sqrt {a^2-b^2 x^3}}dx}{b^{2/3}}-\frac {2 \left (1-\sqrt {3}\right ) \sqrt {2+\sqrt {3}} a^{2/3} \left (a^{2/3}-b^{2/3} x\right ) \sqrt {\frac {a^{2/3} b^{2/3} x+a^{4/3}+b^{4/3} x^2}{\left (\left (1+\sqrt {3}\right ) a^{2/3}-b^{2/3} x\right )^2}} \operatorname {EllipticF}\left (\arcsin \left (\frac {\left (1-\sqrt {3}\right ) a^{2/3}-b^{2/3} x}{\left (1+\sqrt {3}\right ) a^{2/3}-b^{2/3} x}\right ),-7-4 \sqrt {3}\right )}{\sqrt [4]{3} b^{4/3} \sqrt {\frac {a^{2/3} \left (a^{2/3}-b^{2/3} x\right )}{\left (\left (1+\sqrt {3}\right ) a^{2/3}-b^{2/3} x\right )^2}} \sqrt {a^2-b^2 x^3}}\right )}{7 b^2 \sqrt {a-b x^{3/2}} \sqrt {a+b x^{3/2}}}-\frac {2 x^2 \sqrt {a-b x^{3/2}} \sqrt {a+b x^{3/2}}}{7 b^2}\)

\(\Big \downarrow \) 2416

\(\displaystyle \frac {4 a^2 \sqrt {a^2-b^2 x^3} \left (-\frac {2 \left (1-\sqrt {3}\right ) \sqrt {2+\sqrt {3}} a^{2/3} \left (a^{2/3}-b^{2/3} x\right ) \sqrt {\frac {a^{2/3} b^{2/3} x+a^{4/3}+b^{4/3} x^2}{\left (\left (1+\sqrt {3}\right ) a^{2/3}-b^{2/3} x\right )^2}} \operatorname {EllipticF}\left (\arcsin \left (\frac {\left (1-\sqrt {3}\right ) a^{2/3}-b^{2/3} x}{\left (1+\sqrt {3}\right ) a^{2/3}-b^{2/3} x}\right ),-7-4 \sqrt {3}\right )}{\sqrt [4]{3} b^{4/3} \sqrt {\frac {a^{2/3} \left (a^{2/3}-b^{2/3} x\right )}{\left (\left (1+\sqrt {3}\right ) a^{2/3}-b^{2/3} x\right )^2}} \sqrt {a^2-b^2 x^3}}-\frac {\frac {\sqrt [4]{3} \sqrt {2-\sqrt {3}} a^{2/3} \left (a^{2/3}-b^{2/3} x\right ) \sqrt {\frac {a^{2/3} b^{2/3} x+a^{4/3}+b^{4/3} x^2}{\left (\left (1+\sqrt {3}\right ) a^{2/3}-b^{2/3} x\right )^2}} E\left (\arcsin \left (\frac {\left (1-\sqrt {3}\right ) a^{2/3}-b^{2/3} x}{\left (1+\sqrt {3}\right ) a^{2/3}-b^{2/3} x}\right )|-7-4 \sqrt {3}\right )}{b^{2/3} \sqrt {\frac {a^{2/3} \left (a^{2/3}-b^{2/3} x\right )}{\left (\left (1+\sqrt {3}\right ) a^{2/3}-b^{2/3} x\right )^2}} \sqrt {a^2-b^2 x^3}}-\frac {2 \sqrt {a^2-b^2 x^3}}{b^{2/3} \left (\left (1+\sqrt {3}\right ) a^{2/3}-b^{2/3} x\right )}}{b^{2/3}}\right )}{7 b^2 \sqrt {a-b x^{3/2}} \sqrt {a+b x^{3/2}}}-\frac {2 x^2 \sqrt {a-b x^{3/2}} \sqrt {a+b x^{3/2}}}{7 b^2}\)

Input:

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

Output:

(-2*x^2*Sqrt[a - b*x^(3/2)]*Sqrt[a + b*x^(3/2)])/(7*b^2) + (4*a^2*Sqrt[a^2 
 - b^2*x^3]*(-(((-2*Sqrt[a^2 - b^2*x^3])/(b^(2/3)*((1 + Sqrt[3])*a^(2/3) - 
 b^(2/3)*x)) + (3^(1/4)*Sqrt[2 - Sqrt[3]]*a^(2/3)*(a^(2/3) - b^(2/3)*x)*Sq 
rt[(a^(4/3) + a^(2/3)*b^(2/3)*x + b^(4/3)*x^2)/((1 + Sqrt[3])*a^(2/3) - b^ 
(2/3)*x)^2]*EllipticE[ArcSin[((1 - Sqrt[3])*a^(2/3) - b^(2/3)*x)/((1 + Sqr 
t[3])*a^(2/3) - b^(2/3)*x)], -7 - 4*Sqrt[3]])/(b^(2/3)*Sqrt[(a^(2/3)*(a^(2 
/3) - b^(2/3)*x))/((1 + Sqrt[3])*a^(2/3) - b^(2/3)*x)^2]*Sqrt[a^2 - b^2*x^ 
3]))/b^(2/3)) - (2*(1 - Sqrt[3])*Sqrt[2 + Sqrt[3]]*a^(2/3)*(a^(2/3) - b^(2 
/3)*x)*Sqrt[(a^(4/3) + a^(2/3)*b^(2/3)*x + b^(4/3)*x^2)/((1 + Sqrt[3])*a^( 
2/3) - b^(2/3)*x)^2]*EllipticF[ArcSin[((1 - Sqrt[3])*a^(2/3) - b^(2/3)*x)/ 
((1 + Sqrt[3])*a^(2/3) - b^(2/3)*x)], -7 - 4*Sqrt[3]])/(3^(1/4)*b^(4/3)*Sq 
rt[(a^(2/3)*(a^(2/3) - b^(2/3)*x))/((1 + Sqrt[3])*a^(2/3) - b^(2/3)*x)^2]* 
Sqrt[a^2 - b^2*x^3])))/(7*b^2*Sqrt[a - b*x^(3/2)]*Sqrt[a + b*x^(3/2)])
 

Defintions of rubi rules used

rule 759
Int[1/Sqrt[(a_) + (b_.)*(x_)^3], x_Symbol] :> With[{r = Numer[Rt[b/a, 3]], 
s = Denom[Rt[b/a, 3]]}, Simp[2*Sqrt[2 + Sqrt[3]]*(s + r*x)*(Sqrt[(s^2 - r*s 
*x + r^2*x^2)/((1 + Sqrt[3])*s + r*x)^2]/(3^(1/4)*r*Sqrt[a + b*x^3]*Sqrt[s* 
((s + r*x)/((1 + Sqrt[3])*s + r*x)^2)]))*EllipticF[ArcSin[((1 - Sqrt[3])*s 
+ r*x)/((1 + Sqrt[3])*s + r*x)], -7 - 4*Sqrt[3]], x]] /; FreeQ[{a, b}, x] & 
& PosQ[a]
 

rule 832
Int[(x_)/Sqrt[(a_) + (b_.)*(x_)^3], x_Symbol] :> With[{r = Numer[Rt[b/a, 3] 
], s = Denom[Rt[b/a, 3]]}, Simp[(-(1 - Sqrt[3]))*(s/r)   Int[1/Sqrt[a + b*x 
^3], x], x] + Simp[1/r   Int[((1 - Sqrt[3])*s + r*x)/Sqrt[a + b*x^3], x], x 
]] /; FreeQ[{a, b}, x] && PosQ[a]
 

rule 845
Int[((c_.)*(x_))^(m_)*((a1_) + (b1_.)*(x_)^(n_))^(p_)*((a2_) + (b2_.)*(x_)^ 
(n_))^(p_), x_Symbol] :> Simp[c^(2*n - 1)*(c*x)^(m - 2*n + 1)*(a1 + b1*x^n) 
^(p + 1)*((a2 + b2*x^n)^(p + 1)/(b1*b2*(m + 2*n*p + 1))), x] - Simp[a1*a2*c 
^(2*n)*((m - 2*n + 1)/(b1*b2*(m + 2*n*p + 1)))   Int[(c*x)^(m - 2*n)*(a1 + 
b1*x^n)^p*(a2 + b2*x^n)^p, x], x] /; FreeQ[{a1, b1, a2, b2, c, p}, x] && Eq 
Q[a2*b1 + a1*b2, 0] && IGtQ[2*n, 0] && GtQ[m, 2*n - 1] && NeQ[m + 2*n*p + 1 
, 0] && IntBinomialQ[a1*a2, b1*b2, c, 2*n, m, p, x]
 

rule 890
Int[((c_.)*(x_))^(m_.)*((a1_) + (b1_.)*(x_)^(n_))^(p_)*((a2_) + (b2_.)*(x_) 
^(n_))^(p_), x_Symbol] :> Simp[(a1 + b1*x^n)^FracPart[p]*((a2 + b2*x^n)^Fra 
cPart[p]/(a1*a2 + b1*b2*x^(2*n))^FracPart[p])   Int[(c*x)^m*(a1*a2 + b1*b2* 
x^(2*n))^p, x], x] /; FreeQ[{a1, b1, a2, b2, c, m, n, p}, x] && EqQ[a2*b1 + 
 a1*b2, 0] &&  !IntegerQ[p]
 

rule 2416
Int[((c_) + (d_.)*(x_))/Sqrt[(a_) + (b_.)*(x_)^3], x_Symbol] :> With[{r = N 
umer[Simplify[(1 - Sqrt[3])*(d/c)]], s = Denom[Simplify[(1 - Sqrt[3])*(d/c) 
]]}, Simp[2*d*s^3*(Sqrt[a + b*x^3]/(a*r^2*((1 + Sqrt[3])*s + r*x))), x] - S 
imp[3^(1/4)*Sqrt[2 - Sqrt[3]]*d*s*(s + r*x)*(Sqrt[(s^2 - r*s*x + r^2*x^2)/( 
(1 + Sqrt[3])*s + r*x)^2]/(r^2*Sqrt[a + b*x^3]*Sqrt[s*((s + r*x)/((1 + Sqrt 
[3])*s + r*x)^2)]))*EllipticE[ArcSin[((1 - Sqrt[3])*s + r*x)/((1 + Sqrt[3]) 
*s + r*x)], -7 - 4*Sqrt[3]], x]] /; FreeQ[{a, b, c, d}, x] && PosQ[a] && Eq 
Q[b*c^3 - 2*(5 - 3*Sqrt[3])*a*d^3, 0]
 
Maple [F]

\[\int \frac {x^{4}}{\sqrt {a -b \,x^{\frac {3}{2}}}\, \sqrt {a +b \,x^{\frac {3}{2}}}}d x\]

Input:

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

Output:

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

Fricas [A] (verification not implemented)

Time = 0.09 (sec) , antiderivative size = 65, normalized size of antiderivative = 0.11 \[ \int \frac {x^4}{\sqrt {a-b x^{3/2}} \sqrt {a+b x^{3/2}}} \, dx=-\frac {2 \, {\left (\sqrt {b x^{\frac {3}{2}} + a} \sqrt {-b x^{\frac {3}{2}} + a} b^{2} x^{2} - 4 \, \sqrt {-b^{2}} a^{2} {\rm weierstrassZeta}\left (0, \frac {4 \, a^{2}}{b^{2}}, {\rm weierstrassPInverse}\left (0, \frac {4 \, a^{2}}{b^{2}}, x\right )\right )\right )}}{7 \, b^{4}} \] Input:

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

Output:

-2/7*(sqrt(b*x^(3/2) + a)*sqrt(-b*x^(3/2) + a)*b^2*x^2 - 4*sqrt(-b^2)*a^2* 
weierstrassZeta(0, 4*a^2/b^2, weierstrassPInverse(0, 4*a^2/b^2, x)))/b^4
 

Sympy [A] (verification not implemented)

Time = 15.35 (sec) , antiderivative size = 141, normalized size of antiderivative = 0.24 \[ \int \frac {x^4}{\sqrt {a-b x^{3/2}} \sqrt {a+b x^{3/2}}} \, dx=\frac {i a^{\frac {7}{3}} {G_{6, 6}^{5, 3}\left (\begin {matrix} - \frac {11}{12}, - \frac {5}{12}, 1 & - \frac {2}{3}, - \frac {2}{3}, - \frac {1}{6} \\- \frac {7}{6}, - \frac {11}{12}, - \frac {2}{3}, - \frac {5}{12}, - \frac {1}{6} & 0 \end {matrix} \middle | {\frac {a^{2}}{b^{2} x^{3}}} \right )}}{6 \pi ^{\frac {3}{2}} b^{\frac {10}{3}}} + \frac {a^{\frac {7}{3}} {G_{6, 6}^{2, 6}\left (\begin {matrix} - \frac {5}{3}, - \frac {17}{12}, - \frac {7}{6}, - \frac {11}{12}, - \frac {2}{3}, 1 & \\- \frac {17}{12}, - \frac {11}{12} & - \frac {5}{3}, - \frac {7}{6}, - \frac {7}{6}, 0 \end {matrix} \middle | {\frac {a^{2} e^{- 2 i \pi }}{b^{2} x^{3}}} \right )} e^{- \frac {i \pi }{3}}}{6 \pi ^{\frac {3}{2}} b^{\frac {10}{3}}} \] Input:

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

Output:

I*a**(7/3)*meijerg(((-11/12, -5/12, 1), (-2/3, -2/3, -1/6)), ((-7/6, -11/1 
2, -2/3, -5/12, -1/6), (0,)), a**2/(b**2*x**3))/(6*pi**(3/2)*b**(10/3)) + 
a**(7/3)*meijerg(((-5/3, -17/12, -7/6, -11/12, -2/3, 1), ()), ((-17/12, -1 
1/12), (-5/3, -7/6, -7/6, 0)), a**2*exp_polar(-2*I*pi)/(b**2*x**3))*exp(-I 
*pi/3)/(6*pi**(3/2)*b**(10/3))
 

Maxima [F]

\[ \int \frac {x^4}{\sqrt {a-b x^{3/2}} \sqrt {a+b x^{3/2}}} \, dx=\int { \frac {x^{4}}{\sqrt {b x^{\frac {3}{2}} + a} \sqrt {-b x^{\frac {3}{2}} + a}} \,d x } \] Input:

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

Output:

integrate(x^4/(sqrt(b*x^(3/2) + a)*sqrt(-b*x^(3/2) + a)), x)
 

Giac [F]

\[ \int \frac {x^4}{\sqrt {a-b x^{3/2}} \sqrt {a+b x^{3/2}}} \, dx=\int { \frac {x^{4}}{\sqrt {b x^{\frac {3}{2}} + a} \sqrt {-b x^{\frac {3}{2}} + a}} \,d x } \] Input:

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

Output:

integrate(x^4/(sqrt(b*x^(3/2) + a)*sqrt(-b*x^(3/2) + a)), x)
 

Mupad [F(-1)]

Timed out. \[ \int \frac {x^4}{\sqrt {a-b x^{3/2}} \sqrt {a+b x^{3/2}}} \, dx=\int \frac {x^4}{\sqrt {a+b\,x^{3/2}}\,\sqrt {a-b\,x^{3/2}}} \,d x \] Input:

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

Output:

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

Reduce [F]

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

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

Output:

(2*( - sqrt(sqrt(x)*b*x + a)*sqrt( - sqrt(x)*b*x + a)*x**2 + 2*int((sqrt(s 
qrt(x)*b*x + a)*sqrt( - sqrt(x)*b*x + a)*x)/(a**2 - b**2*x**3),x)*a**2))/( 
7*b**2)