\(\int \frac {(a+b x^2)^{7/6}}{x^8} \, dx\) [1097]

Optimal result

Integrand size = 15, antiderivative size = 337 \[ \int \frac {\left (a+b x^2\right )^{7/6}}{x^8} \, dx=-\frac {b \sqrt [6]{a+b x^2}}{15 x^5}-\frac {b^2 \sqrt [6]{a+b x^2}}{135 a x^3}+\frac {8 b^3 \sqrt [6]{a+b x^2}}{405 a^2 x}-\frac {\left (a+b x^2\right )^{7/6}}{7 x^7}+\frac {8 b^3 \sqrt [6]{a+b x^2} \left (\sqrt [3]{a}-\sqrt [3]{a+b x^2}\right ) \sqrt {\frac {a^{2/3}+\sqrt [3]{a} \sqrt [3]{a+b x^2}+\left (a+b x^2\right )^{2/3}}{\left (\sqrt [3]{a}-\left (1+\sqrt {3}\right ) \sqrt [3]{a+b x^2}\right )^2}} \operatorname {EllipticF}\left (\arccos \left (\frac {\sqrt [3]{a}-\left (1-\sqrt {3}\right ) \sqrt [3]{a+b x^2}}{\sqrt [3]{a}-\left (1+\sqrt {3}\right ) \sqrt [3]{a+b x^2}}\right ),\frac {1}{4} \left (2+\sqrt {3}\right )\right )}{405 \sqrt [4]{3} a^{7/3} x \sqrt {-\frac {\sqrt [3]{a+b x^2} \left (\sqrt [3]{a}-\sqrt [3]{a+b x^2}\right )}{\left (\sqrt [3]{a}-\left (1+\sqrt {3}\right ) \sqrt [3]{a+b x^2}\right )^2}}} \] Output:

-1/15*b*(b*x^2+a)^(1/6)/x^5-1/135*b^2*(b*x^2+a)^(1/6)/a/x^3+8/405*b^3*(b*x 
^2+a)^(1/6)/a^2/x-1/7*(b*x^2+a)^(7/6)/x^7+8/1215*b^3*(b*x^2+a)^(1/6)*(a^(1 
/3)-(b*x^2+a)^(1/3))*((a^(2/3)+a^(1/3)*(b*x^2+a)^(1/3)+(b*x^2+a)^(2/3))/(a 
^(1/3)-(1+3^(1/2))*(b*x^2+a)^(1/3))^2)^(1/2)*InverseJacobiAM(arccos((a^(1/ 
3)-(1-3^(1/2))*(b*x^2+a)^(1/3))/(a^(1/3)-(1+3^(1/2))*(b*x^2+a)^(1/3))),1/4 
*6^(1/2)+1/4*2^(1/2))*3^(3/4)/a^(7/3)/x/(-(b*x^2+a)^(1/3)*(a^(1/3)-(b*x^2+ 
a)^(1/3))/(a^(1/3)-(1+3^(1/2))*(b*x^2+a)^(1/3))^2)^(1/2)
 

Mathematica [C] (verified)

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

Time = 10.01 (sec) , antiderivative size = 52, normalized size of antiderivative = 0.15 \[ \int \frac {\left (a+b x^2\right )^{7/6}}{x^8} \, dx=-\frac {a \sqrt [6]{a+b x^2} \operatorname {Hypergeometric2F1}\left (-\frac {7}{2},-\frac {7}{6},-\frac {5}{2},-\frac {b x^2}{a}\right )}{7 x^7 \sqrt [6]{1+\frac {b x^2}{a}}} \] Input:

Integrate[(a + b*x^2)^(7/6)/x^8,x]
 

Output:

-1/7*(a*(a + b*x^2)^(1/6)*Hypergeometric2F1[-7/2, -7/6, -5/2, -((b*x^2)/a) 
])/(x^7*(1 + (b*x^2)/a)^(1/6))
 

Rubi [A] (warning: unable to verify)

Time = 0.32 (sec) , antiderivative size = 431, normalized size of antiderivative = 1.28, number of steps used = 8, number of rules used = 7, \(\frac {\text {number of rules}}{\text {integrand size}}\) = 0.467, Rules used = {247, 247, 264, 264, 236, 234, 760}

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 {\left (a+b x^2\right )^{7/6}}{x^8} \, dx\)

\(\Big \downarrow \) 247

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

\(\Big \downarrow \) 247

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

\(\Big \downarrow \) 264

\(\displaystyle \frac {1}{3} b \left (\frac {1}{15} b \left (-\frac {8 b \int \frac {1}{x^2 \left (b x^2+a\right )^{5/6}}dx}{9 a}-\frac {\sqrt [6]{a+b x^2}}{3 a x^3}\right )-\frac {\sqrt [6]{a+b x^2}}{5 x^5}\right )-\frac {\left (a+b x^2\right )^{7/6}}{7 x^7}\)

\(\Big \downarrow \) 264

\(\displaystyle \frac {1}{3} b \left (\frac {1}{15} b \left (-\frac {8 b \left (-\frac {2 b \int \frac {1}{\left (b x^2+a\right )^{5/6}}dx}{3 a}-\frac {\sqrt [6]{a+b x^2}}{a x}\right )}{9 a}-\frac {\sqrt [6]{a+b x^2}}{3 a x^3}\right )-\frac {\sqrt [6]{a+b x^2}}{5 x^5}\right )-\frac {\left (a+b x^2\right )^{7/6}}{7 x^7}\)

\(\Big \downarrow \) 236

\(\displaystyle \frac {1}{3} b \left (\frac {1}{15} b \left (-\frac {8 b \left (-\frac {2 b \int \frac {1}{\left (1-\frac {b x^2}{b x^2+a}\right )^{2/3}}d\frac {x}{\sqrt {b x^2+a}}}{3 a \sqrt [3]{\frac {a}{a+b x^2}} \sqrt [3]{a+b x^2}}-\frac {\sqrt [6]{a+b x^2}}{a x}\right )}{9 a}-\frac {\sqrt [6]{a+b x^2}}{3 a x^3}\right )-\frac {\sqrt [6]{a+b x^2}}{5 x^5}\right )-\frac {\left (a+b x^2\right )^{7/6}}{7 x^7}\)

\(\Big \downarrow \) 234

\(\displaystyle \frac {1}{3} b \left (\frac {1}{15} b \left (-\frac {8 b \left (\frac {\sqrt {-\frac {b x^2}{a+b x^2}} \sqrt [6]{a+b x^2} \int \frac {1}{\sqrt {\frac {x^3}{\left (b x^2+a\right )^{3/2}}-1}}d\sqrt [3]{1-\frac {b x^2}{b x^2+a}}}{a x \sqrt [3]{\frac {a}{a+b x^2}}}-\frac {\sqrt [6]{a+b x^2}}{a x}\right )}{9 a}-\frac {\sqrt [6]{a+b x^2}}{3 a x^3}\right )-\frac {\sqrt [6]{a+b x^2}}{5 x^5}\right )-\frac {\left (a+b x^2\right )^{7/6}}{7 x^7}\)

\(\Big \downarrow \) 760

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

Input:

Int[(a + b*x^2)^(7/6)/x^8,x]
 

Output:

-1/7*(a + b*x^2)^(7/6)/x^7 + (b*(-1/5*(a + b*x^2)^(1/6)/x^5 + (b*(-1/3*(a 
+ b*x^2)^(1/6)/(a*x^3) - (8*b*(-((a + b*x^2)^(1/6)/(a*x)) - (2*Sqrt[2 - Sq 
rt[3]]*Sqrt[-((b*x^2)/(a + b*x^2))]*(a + b*x^2)^(1/6)*(1 - (1 - (b*x^2)/(a 
 + b*x^2))^(1/3))*Sqrt[(1 + x^2/(a + b*x^2) + (1 - (b*x^2)/(a + b*x^2))^(1 
/3))/(1 - Sqrt[3] - (1 - (b*x^2)/(a + b*x^2))^(1/3))^2]*EllipticF[ArcSin[( 
1 + Sqrt[3] - (1 - (b*x^2)/(a + b*x^2))^(1/3))/(1 - Sqrt[3] - (1 - (b*x^2) 
/(a + b*x^2))^(1/3))], -7 + 4*Sqrt[3]])/(3^(1/4)*a*x*(a/(a + b*x^2))^(1/3) 
*Sqrt[-1 + x^3/(a + b*x^2)^(3/2)]*Sqrt[-((1 - (1 - (b*x^2)/(a + b*x^2))^(1 
/3))/(1 - Sqrt[3] - (1 - (b*x^2)/(a + b*x^2))^(1/3))^2)])))/(9*a)))/15))/3
 

Defintions of rubi rules used

rule 234

Int[((a_) + (b_.)*(x_)^2)^(-2/3), x_Symbol] :> Simp[3*(Sqrt[b*x^2]/(2*b*x)) 
   Subst[Int[1/Sqrt[-a + x^3], x], x, (a + b*x^2)^(1/3)], x] /; FreeQ[{a, b 
}, x]
 

rule 236

Int[((a_) + (b_.)*(x_)^2)^(-5/6), x_Symbol] :> Simp[1/((a/(a + b*x^2))^(1/3 
)*(a + b*x^2)^(1/3))   Subst[Int[1/(1 - b*x^2)^(2/3), x], x, x/Sqrt[a + b*x 
^2]], x] /; FreeQ[{a, b}, x]
 

rule 247

Int[((c_.)*(x_))^(m_.)*((a_) + (b_.)*(x_)^2)^(p_), x_Symbol] :> Simp[(c*x)^ 
(m + 1)*((a + b*x^2)^p/(c*(m + 1))), x] - Simp[2*b*(p/(c^2*(m + 1)))   Int[ 
(c*x)^(m + 2)*(a + b*x^2)^(p - 1), x], x] /; FreeQ[{a, b, c}, x] && GtQ[p, 
0] && LtQ[m, -1] &&  !ILtQ[(m + 2*p + 3)/2, 0] && IntBinomialQ[a, b, c, 2, 
m, p, x]
 

rule 264

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

rule 760

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 
] && NegQ[a]
 

Maple [F]

\[\int \frac {\left (b \,x^{2}+a \right )^{\frac {7}{6}}}{x^{8}}d x\]

Input:

int((b*x^2+a)^(7/6)/x^8,x)
 

Output:

int((b*x^2+a)^(7/6)/x^8,x)
 

Fricas [F]

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

integrate((b*x^2+a)^(7/6)/x^8,x, algorithm="fricas")
 

Output:

integral((b*x^2 + a)^(7/6)/x^8, x)
 

Sympy [C] (verification not implemented)

Result contains complex when optimal does not.

Time = 1.04 (sec) , antiderivative size = 34, normalized size of antiderivative = 0.10 \[ \int \frac {\left (a+b x^2\right )^{7/6}}{x^8} \, dx=- \frac {a^{\frac {7}{6}} {{}_{2}F_{1}\left (\begin {matrix} - \frac {7}{2}, - \frac {7}{6} \\ - \frac {5}{2} \end {matrix}\middle | {\frac {b x^{2} e^{i \pi }}{a}} \right )}}{7 x^{7}} \] Input:

integrate((b*x**2+a)**(7/6)/x**8,x)
 

Output:

-a**(7/6)*hyper((-7/2, -7/6), (-5/2,), b*x**2*exp_polar(I*pi)/a)/(7*x**7)
 

Maxima [F]

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

integrate((b*x^2+a)^(7/6)/x^8,x, algorithm="maxima")
 

Output:

integrate((b*x^2 + a)^(7/6)/x^8, x)
 

Giac [F]

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

integrate((b*x^2+a)^(7/6)/x^8,x, algorithm="giac")
                                                                                    
                                                                                    
 

Output:

integrate((b*x^2 + a)^(7/6)/x^8, x)
 

Mupad [F(-1)]

Timed out. \[ \int \frac {\left (a+b x^2\right )^{7/6}}{x^8} \, dx=\int \frac {{\left (b\,x^2+a\right )}^{7/6}}{x^8} \,d x \] Input:

int((a + b*x^2)^(7/6)/x^8,x)
 

Output:

int((a + b*x^2)^(7/6)/x^8, x)
 

Reduce [F]

\[ \int \frac {\left (a+b x^2\right )^{7/6}}{x^8} \, dx=\frac {-7 \left (b \,x^{2}+a \right )^{\frac {5}{6}} a -8 \left (b \,x^{2}+a \right )^{\frac {5}{6}} b \,x^{2}-\left (b \,x^{2}+a \right )^{\frac {2}{3}} \left (\int \frac {\left (b \,x^{2}+a \right )^{\frac {7}{6}}}{b^{2} x^{12}+2 a b \,x^{10}+a^{2} x^{8}}d x \right ) a^{2} x^{7}}{48 \left (b \,x^{2}+a \right )^{\frac {2}{3}} x^{7}} \] Input:

int((b*x^2+a)^(7/6)/x^8,x)
 

Output:

( - 7*(a + b*x**2)**(5/6)*a - 8*(a + b*x**2)**(5/6)*b*x**2 - (a + b*x**2)* 
*(2/3)*int((a + b*x**2)**(7/6)/(a**2*x**8 + 2*a*b*x**10 + b**2*x**12),x)*a 
**2*x**7)/(48*(a + b*x**2)**(2/3)*x**7)