\(\int \frac {(a+b x^3+c x^6)^p}{x} \, dx\) [249]

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

Optimal result

Integrand size = 18, antiderivative size = 157 \[ \int \frac {\left (a+b x^3+c x^6\right )^p}{x} \, dx=\frac {2^{-1+2 p} \left (\frac {b-\sqrt {b^2-4 a c}+2 c x^3}{c x^3}\right )^{-p} \left (\frac {b+\sqrt {b^2-4 a c}+2 c x^3}{c x^3}\right )^{-p} \left (a+b x^3+c x^6\right )^p \operatorname {AppellF1}\left (-2 p,-p,-p,1-2 p,-\frac {b-\sqrt {b^2-4 a c}}{2 c x^3},-\frac {b+\sqrt {b^2-4 a c}}{2 c x^3}\right )}{3 p} \] Output:

1/3*2^(-1+2*p)*(c*x^6+b*x^3+a)^p*AppellF1(-2*p,-p,-p,1-2*p,-1/2*(b-(-4*a*c 
+b^2)^(1/2))/c/x^3,-1/2*(b+(-4*a*c+b^2)^(1/2))/c/x^3)/p/(((b-(-4*a*c+b^2)^ 
(1/2)+2*c*x^3)/c/x^3)^p)/(((b+(-4*a*c+b^2)^(1/2)+2*c*x^3)/c/x^3)^p)
 

Mathematica [A] (verified)

Time = 0.35 (sec) , antiderivative size = 157, normalized size of antiderivative = 1.00 \[ \int \frac {\left (a+b x^3+c x^6\right )^p}{x} \, dx=\frac {2^{-1+2 p} \left (\frac {b-\sqrt {b^2-4 a c}+2 c x^3}{c x^3}\right )^{-p} \left (\frac {b+\sqrt {b^2-4 a c}+2 c x^3}{c x^3}\right )^{-p} \left (a+b x^3+c x^6\right )^p \operatorname {AppellF1}\left (-2 p,-p,-p,1-2 p,-\frac {b+\sqrt {b^2-4 a c}}{2 c x^3},\frac {-b+\sqrt {b^2-4 a c}}{2 c x^3}\right )}{3 p} \] Input:

Integrate[(a + b*x^3 + c*x^6)^p/x,x]
 

Output:

(2^(-1 + 2*p)*(a + b*x^3 + c*x^6)^p*AppellF1[-2*p, -p, -p, 1 - 2*p, -1/2*( 
b + Sqrt[b^2 - 4*a*c])/(c*x^3), (-b + Sqrt[b^2 - 4*a*c])/(2*c*x^3)])/(3*p* 
((b - Sqrt[b^2 - 4*a*c] + 2*c*x^3)/(c*x^3))^p*((b + Sqrt[b^2 - 4*a*c] + 2* 
c*x^3)/(c*x^3))^p)
 

Rubi [A] (verified)

Time = 0.32 (sec) , antiderivative size = 157, normalized size of antiderivative = 1.00, number of steps used = 4, number of rules used = 3, \(\frac {\text {number of rules}}{\text {integrand size}}\) = 0.167, Rules used = {1693, 1178, 150}

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^3+c x^6\right )^p}{x} \, dx\)

\(\Big \downarrow \) 1693

\(\displaystyle \frac {1}{3} \int \frac {\left (c x^6+b x^3+a\right )^p}{x^3}dx^3\)

\(\Big \downarrow \) 1178

\(\displaystyle -\frac {1}{3} 4^p \left (\frac {1}{x^3}\right )^{2 p} \left (\frac {-\sqrt {b^2-4 a c}+b+2 c x^3}{c x^3}\right )^{-p} \left (\frac {\sqrt {b^2-4 a c}+b+2 c x^3}{c x^3}\right )^{-p} \left (a+b x^3+c x^6\right )^p \int \left (\frac {b-\sqrt {b^2-4 a c}}{2 c x^3}+1\right )^p \left (\frac {b+\sqrt {b^2-4 a c}}{2 c x^3}+1\right )^p \left (\frac {1}{x^3}\right )^{-2 p-1}d\frac {1}{x^3}\)

\(\Big \downarrow \) 150

\(\displaystyle \frac {2^{2 p-1} \left (\frac {-\sqrt {b^2-4 a c}+b+2 c x^3}{c x^3}\right )^{-p} \left (\frac {\sqrt {b^2-4 a c}+b+2 c x^3}{c x^3}\right )^{-p} \left (a+b x^3+c x^6\right )^p \operatorname {AppellF1}\left (-2 p,-p,-p,1-2 p,-\frac {b-\sqrt {b^2-4 a c}}{2 c x^3},-\frac {b+\sqrt {b^2-4 a c}}{2 c x^3}\right )}{3 p}\)

Input:

Int[(a + b*x^3 + c*x^6)^p/x,x]
 

Output:

(2^(-1 + 2*p)*(a + b*x^3 + c*x^6)^p*AppellF1[-2*p, -p, -p, 1 - 2*p, -1/2*( 
b - Sqrt[b^2 - 4*a*c])/(c*x^3), -1/2*(b + Sqrt[b^2 - 4*a*c])/(c*x^3)])/(3* 
p*((b - Sqrt[b^2 - 4*a*c] + 2*c*x^3)/(c*x^3))^p*((b + Sqrt[b^2 - 4*a*c] + 
2*c*x^3)/(c*x^3))^p)
 

Defintions of rubi rules used

rule 150
Int[((b_.)*(x_))^(m_)*((c_) + (d_.)*(x_))^(n_)*((e_) + (f_.)*(x_))^(p_), x_ 
] :> Simp[c^n*e^p*((b*x)^(m + 1)/(b*(m + 1)))*AppellF1[m + 1, -n, -p, m + 2 
, (-d)*(x/c), (-f)*(x/e)], x] /; FreeQ[{b, c, d, e, f, m, n, p}, x] &&  !In 
tegerQ[m] &&  !IntegerQ[n] && GtQ[c, 0] && (IntegerQ[p] || GtQ[e, 0])
 

rule 1178
Int[((d_.) + (e_.)*(x_))^(m_)*((a_.) + (b_.)*(x_) + (c_.)*(x_)^2)^(p_), x_S 
ymbol] :> With[{q = Rt[b^2 - 4*a*c, 2]}, Simp[(-(1/(d + e*x))^(2*p))*((a + 
b*x + c*x^2)^p/(e*(e*((b - q + 2*c*x)/(2*c*(d + e*x))))^p*(e*((b + q + 2*c* 
x)/(2*c*(d + e*x))))^p))   Subst[Int[x^(-m - 2*(p + 1))*Simp[1 - (d - e*((b 
 - q)/(2*c)))*x, x]^p*Simp[1 - (d - e*((b + q)/(2*c)))*x, x]^p, x], x, 1/(d 
 + e*x)], x]] /; FreeQ[{a, b, c, d, e, p}, x] && ILtQ[m, 0]
 

rule 1693
Int[(x_)^(m_.)*((a_) + (c_.)*(x_)^(n2_.) + (b_.)*(x_)^(n_))^(p_.), x_Symbol 
] :> Simp[1/n   Subst[Int[x^(Simplify[(m + 1)/n] - 1)*(a + b*x + c*x^2)^p, 
x], x, x^n], x] /; FreeQ[{a, b, c, m, n, p}, x] && EqQ[n2, 2*n] && IntegerQ 
[Simplify[(m + 1)/n]]
 
Maple [F]

\[\int \frac {\left (c \,x^{6}+b \,x^{3}+a \right )^{p}}{x}d x\]

Input:

int((c*x^6+b*x^3+a)^p/x,x)
 

Output:

int((c*x^6+b*x^3+a)^p/x,x)
 

Fricas [F]

\[ \int \frac {\left (a+b x^3+c x^6\right )^p}{x} \, dx=\int { \frac {{\left (c x^{6} + b x^{3} + a\right )}^{p}}{x} \,d x } \] Input:

integrate((c*x^6+b*x^3+a)^p/x,x, algorithm="fricas")
 

Output:

integral((c*x^6 + b*x^3 + a)^p/x, x)
 

Sympy [F(-1)]

Timed out. \[ \int \frac {\left (a+b x^3+c x^6\right )^p}{x} \, dx=\text {Timed out} \] Input:

integrate((c*x**6+b*x**3+a)**p/x,x)
 

Output:

Timed out
 

Maxima [F]

\[ \int \frac {\left (a+b x^3+c x^6\right )^p}{x} \, dx=\int { \frac {{\left (c x^{6} + b x^{3} + a\right )}^{p}}{x} \,d x } \] Input:

integrate((c*x^6+b*x^3+a)^p/x,x, algorithm="maxima")
 

Output:

integrate((c*x^6 + b*x^3 + a)^p/x, x)
 

Giac [F]

\[ \int \frac {\left (a+b x^3+c x^6\right )^p}{x} \, dx=\int { \frac {{\left (c x^{6} + b x^{3} + a\right )}^{p}}{x} \,d x } \] Input:

integrate((c*x^6+b*x^3+a)^p/x,x, algorithm="giac")
 

Output:

integrate((c*x^6 + b*x^3 + a)^p/x, x)
 

Mupad [F(-1)]

Timed out. \[ \int \frac {\left (a+b x^3+c x^6\right )^p}{x} \, dx=\int \frac {{\left (c\,x^6+b\,x^3+a\right )}^p}{x} \,d x \] Input:

int((a + b*x^3 + c*x^6)^p/x,x)
 

Output:

int((a + b*x^3 + c*x^6)^p/x, x)
 

Reduce [F]

\[ \int \frac {\left (a+b x^3+c x^6\right )^p}{x} \, dx=\frac {\left (c \,x^{6}+b \,x^{3}+a \right )^{p}+3 \left (\int \frac {\left (c \,x^{6}+b \,x^{3}+a \right )^{p}}{c \,x^{7}+b \,x^{4}+a x}d x \right ) a p -3 \left (\int \frac {\left (c \,x^{6}+b \,x^{3}+a \right )^{p} x^{5}}{c \,x^{6}+b \,x^{3}+a}d x \right ) c p}{3 p} \] Input:

int((c*x^6+b*x^3+a)^p/x,x)
 

Output:

((a + b*x**3 + c*x**6)**p + 3*int((a + b*x**3 + c*x**6)**p/(a*x + b*x**4 + 
 c*x**7),x)*a*p - 3*int(((a + b*x**3 + c*x**6)**p*x**5)/(a + b*x**3 + c*x* 
*6),x)*c*p)/(3*p)