\(\int \frac {\sqrt {a-b x^4}}{x^2} \, dx\) [181]

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

Optimal result

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

-(-b*x^4+a)^(1/2)/x-2*a^(3/4)*b^(1/4)*(1-b*x^4/a)^(1/2)*EllipticE(b^(1/4)* 
x/a^(1/4),I)/(-b*x^4+a)^(1/2)+2*a^(3/4)*b^(1/4)*(1-b*x^4/a)^(1/2)*Elliptic 
F(b^(1/4)*x/a^(1/4),I)/(-b*x^4+a)^(1/2)
 

Mathematica [C] (verified)

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

Time = 8.47 (sec) , antiderivative size = 50, normalized size of antiderivative = 0.40 \[ \int \frac {\sqrt {a-b x^4}}{x^2} \, dx=-\frac {\sqrt {a-b x^4} \operatorname {Hypergeometric2F1}\left (-\frac {1}{2},-\frac {1}{4},\frac {3}{4},\frac {b x^4}{a}\right )}{x \sqrt {1-\frac {b x^4}{a}}} \] Input:

Integrate[Sqrt[a - b*x^4]/x^2,x]
 

Output:

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

Rubi [A] (verified)

Time = 0.50 (sec) , antiderivative size = 129, normalized size of antiderivative = 1.02, number of steps used = 8, number of rules used = 8, \(\frac {\text {number of rules}}{\text {integrand size}}\) = 0.500, Rules used = {809, 836, 27, 765, 762, 1390, 1389, 327}

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

\(\Big \downarrow \) 809

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

\(\Big \downarrow \) 836

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

\(\Big \downarrow \) 27

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

\(\Big \downarrow \) 765

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

\(\Big \downarrow \) 762

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

\(\Big \downarrow \) 1390

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

\(\Big \downarrow \) 1389

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

\(\Big \downarrow \) 327

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

Input:

Int[Sqrt[a - b*x^4]/x^2,x]
 

Output:

-(Sqrt[a - b*x^4]/x) - 2*b*((a^(3/4)*Sqrt[1 - (b*x^4)/a]*EllipticE[ArcSin[ 
(b^(1/4)*x)/a^(1/4)], -1])/(b^(3/4)*Sqrt[a - b*x^4]) - (a^(3/4)*Sqrt[1 - ( 
b*x^4)/a]*EllipticF[ArcSin[(b^(1/4)*x)/a^(1/4)], -1])/(b^(3/4)*Sqrt[a - b* 
x^4]))
 

Defintions of rubi rules used

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 327
Int[Sqrt[(a_) + (b_.)*(x_)^2]/Sqrt[(c_) + (d_.)*(x_)^2], x_Symbol] :> Simp[ 
(Sqrt[a]/(Sqrt[c]*Rt[-d/c, 2]))*EllipticE[ArcSin[Rt[-d/c, 2]*x], b*(c/(a*d) 
)], x] /; FreeQ[{a, b, c, d}, x] && NegQ[d/c] && GtQ[c, 0] && GtQ[a, 0]
 

rule 762
Int[1/Sqrt[(a_) + (b_.)*(x_)^4], x_Symbol] :> Simp[(1/(Sqrt[a]*Rt[-b/a, 4]) 
)*EllipticF[ArcSin[Rt[-b/a, 4]*x], -1], x] /; FreeQ[{a, b}, x] && NegQ[b/a] 
 && GtQ[a, 0]
 

rule 765
Int[1/Sqrt[(a_) + (b_.)*(x_)^4], x_Symbol] :> Simp[Sqrt[1 + b*(x^4/a)]/Sqrt 
[a + b*x^4]   Int[1/Sqrt[1 + b*(x^4/a)], x], x] /; FreeQ[{a, b}, x] && NegQ 
[b/a] &&  !GtQ[a, 0]
 

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

rule 836
Int[(x_)^2/Sqrt[(a_) + (b_.)*(x_)^4], x_Symbol] :> With[{q = Rt[-b/a, 2]}, 
Simp[-q^(-1)   Int[1/Sqrt[a + b*x^4], x], x] + Simp[1/q   Int[(1 + q*x^2)/S 
qrt[a + b*x^4], x], x]] /; FreeQ[{a, b}, x] && NegQ[b/a]
 

rule 1389
Int[((d_) + (e_.)*(x_)^2)/Sqrt[(a_) + (c_.)*(x_)^4], x_Symbol] :> Simp[d/Sq 
rt[a]   Int[Sqrt[1 + e*(x^2/d)]/Sqrt[1 - e*(x^2/d)], x], x] /; FreeQ[{a, c, 
 d, e}, x] && EqQ[c*d^2 + a*e^2, 0] && NegQ[c/a] && GtQ[a, 0]
 

rule 1390
Int[((d_) + (e_.)*(x_)^2)/Sqrt[(a_) + (c_.)*(x_)^4], x_Symbol] :> Simp[Sqrt 
[1 + c*(x^4/a)]/Sqrt[a + c*x^4]   Int[(d + e*x^2)/Sqrt[1 + c*(x^4/a)], x], 
x] /; FreeQ[{a, c, d, e}, x] && EqQ[c*d^2 + a*e^2, 0] && NegQ[c/a] &&  !GtQ 
[a, 0] &&  !(LtQ[a, 0] && GtQ[c, 0])
 
Maple [A] (verified)

Time = 0.61 (sec) , antiderivative size = 104, normalized size of antiderivative = 0.83

method result size
default \(-\frac {\sqrt {-b \,x^{4}+a}}{x}+\frac {2 \sqrt {b}\, \sqrt {a}\, \sqrt {1-\frac {\sqrt {b}\, x^{2}}{\sqrt {a}}}\, \sqrt {1+\frac {\sqrt {b}\, x^{2}}{\sqrt {a}}}\, \left (\operatorname {EllipticF}\left (x \sqrt {\frac {\sqrt {b}}{\sqrt {a}}}, i\right )-\operatorname {EllipticE}\left (x \sqrt {\frac {\sqrt {b}}{\sqrt {a}}}, i\right )\right )}{\sqrt {\frac {\sqrt {b}}{\sqrt {a}}}\, \sqrt {-b \,x^{4}+a}}\) \(104\)
risch \(-\frac {\sqrt {-b \,x^{4}+a}}{x}+\frac {2 \sqrt {b}\, \sqrt {a}\, \sqrt {1-\frac {\sqrt {b}\, x^{2}}{\sqrt {a}}}\, \sqrt {1+\frac {\sqrt {b}\, x^{2}}{\sqrt {a}}}\, \left (\operatorname {EllipticF}\left (x \sqrt {\frac {\sqrt {b}}{\sqrt {a}}}, i\right )-\operatorname {EllipticE}\left (x \sqrt {\frac {\sqrt {b}}{\sqrt {a}}}, i\right )\right )}{\sqrt {\frac {\sqrt {b}}{\sqrt {a}}}\, \sqrt {-b \,x^{4}+a}}\) \(104\)
elliptic \(-\frac {\sqrt {-b \,x^{4}+a}}{x}+\frac {2 \sqrt {b}\, \sqrt {a}\, \sqrt {1-\frac {\sqrt {b}\, x^{2}}{\sqrt {a}}}\, \sqrt {1+\frac {\sqrt {b}\, x^{2}}{\sqrt {a}}}\, \left (\operatorname {EllipticF}\left (x \sqrt {\frac {\sqrt {b}}{\sqrt {a}}}, i\right )-\operatorname {EllipticE}\left (x \sqrt {\frac {\sqrt {b}}{\sqrt {a}}}, i\right )\right )}{\sqrt {\frac {\sqrt {b}}{\sqrt {a}}}\, \sqrt {-b \,x^{4}+a}}\) \(104\)

Input:

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

Output:

-(-b*x^4+a)^(1/2)/x+2*b^(1/2)*a^(1/2)/(1/a^(1/2)*b^(1/2))^(1/2)*(1-b^(1/2) 
*x^2/a^(1/2))^(1/2)*(1+b^(1/2)*x^2/a^(1/2))^(1/2)/(-b*x^4+a)^(1/2)*(Ellipt 
icF(x*(1/a^(1/2)*b^(1/2))^(1/2),I)-EllipticE(x*(1/a^(1/2)*b^(1/2))^(1/2),I 
))
 

Fricas [F]

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

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

Output:

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

Sympy [A] (verification not implemented)

Time = 0.59 (sec) , antiderivative size = 36, normalized size of antiderivative = 0.29 \[ \int \frac {\sqrt {a-b x^4}}{x^2} \, dx=\frac {i \sqrt {b} x \Gamma \left (\frac {1}{4}\right ) {{}_{2}F_{1}\left (\begin {matrix} - \frac {1}{2}, - \frac {1}{4} \\ \frac {3}{4} \end {matrix}\middle | {\frac {a}{b x^{4}}} \right )}}{4 \Gamma \left (\frac {5}{4}\right )} \] Input:

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

Output:

I*sqrt(b)*x*gamma(1/4)*hyper((-1/2, -1/4), (3/4,), a/(b*x**4))/(4*gamma(5/ 
4))
 

Maxima [F]

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

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

Output:

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

Giac [F]

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

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

Output:

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

Mupad [B] (verification not implemented)

Time = 0.52 (sec) , antiderivative size = 40, normalized size of antiderivative = 0.32 \[ \int \frac {\sqrt {a-b x^4}}{x^2} \, dx=\frac {\sqrt {a-b\,x^4}\,{{}}_2{\mathrm {F}}_1\left (-\frac {1}{2},-\frac {1}{4};\ \frac {3}{4};\ \frac {a}{b\,x^4}\right )}{x\,\sqrt {1-\frac {a}{b\,x^4}}} \] Input:

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

Output:

((a - b*x^4)^(1/2)*hypergeom([-1/2, -1/4], 3/4, a/(b*x^4)))/(x*(1 - a/(b*x 
^4))^(1/2))
 

Reduce [F]

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

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

Output:

(sqrt(a - b*x**4) + 2*int(sqrt(a - b*x**4)/(a*x**2 - b*x**6),x)*a*x)/x