\(\int (d+e x^3)^q (d^2-e^2 x^6)^p \, dx\) [12]

Optimal result

Integrand size = 24, antiderivative size = 92 \[ \int \left (d+e x^3\right )^q \left (d^2-e^2 x^6\right )^p \, dx=x \left (d+e x^3\right )^q \left (1-\frac {e x^3}{d}\right )^{-p} \left (1+\frac {e x^3}{d}\right )^{-p-q} \left (d^2-e^2 x^6\right )^p \operatorname {AppellF1}\left (\frac {1}{3},-p,-p-q,\frac {4}{3},\frac {e x^3}{d},-\frac {e x^3}{d}\right ) \] Output:

x*(e*x^3+d)^q*(1+e*x^3/d)^(-p-q)*(-e^2*x^6+d^2)^p*AppellF1(1/3,-p,-p-q,4/3 
,e*x^3/d,-e*x^3/d)/((1-e*x^3/d)^p)
 

Mathematica [F]

\[ \int \left (d+e x^3\right )^q \left (d^2-e^2 x^6\right )^p \, dx=\int \left (d+e x^3\right )^q \left (d^2-e^2 x^6\right )^p \, dx \] Input:

Integrate[(d + e*x^3)^q*(d^2 - e^2*x^6)^p,x]
 

Output:

Integrate[(d + e*x^3)^q*(d^2 - e^2*x^6)^p, x]
 

Rubi [A] (verified)

Time = 0.26 (sec) , antiderivative size = 92, normalized size of antiderivative = 1.00, number of steps used = 4, number of rules used = 4, \(\frac {\text {number of rules}}{\text {integrand size}}\) = 0.167, Rules used = {1396, 937, 937, 936}

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 \left (d^2-e^2 x^6\right )^p \left (d+e x^3\right )^q \, dx\)

\(\Big \downarrow \) 1396

\(\displaystyle \left (d-e x^3\right )^{-p} \left (d+e x^3\right )^{-p} \left (d^2-e^2 x^6\right )^p \int \left (d-e x^3\right )^p \left (e x^3+d\right )^{p+q}dx\)

\(\Big \downarrow \) 937

\(\displaystyle \left (d+e x^3\right )^{-p} \left (1-\frac {e x^3}{d}\right )^{-p} \left (d^2-e^2 x^6\right )^p \int \left (e x^3+d\right )^{p+q} \left (1-\frac {e x^3}{d}\right )^pdx\)

\(\Big \downarrow \) 937

\(\displaystyle \left (1-\frac {e x^3}{d}\right )^{-p} \left (d^2-e^2 x^6\right )^p \left (d+e x^3\right )^q \left (\frac {e x^3}{d}+1\right )^{-p-q} \int \left (1-\frac {e x^3}{d}\right )^p \left (\frac {e x^3}{d}+1\right )^{p+q}dx\)

\(\Big \downarrow \) 936

\(\displaystyle x \left (1-\frac {e x^3}{d}\right )^{-p} \left (d^2-e^2 x^6\right )^p \left (d+e x^3\right )^q \left (\frac {e x^3}{d}+1\right )^{-p-q} \operatorname {AppellF1}\left (\frac {1}{3},-p,-p-q,\frac {4}{3},\frac {e x^3}{d},-\frac {e x^3}{d}\right )\)

Input:

Int[(d + e*x^3)^q*(d^2 - e^2*x^6)^p,x]
 

Output:

(x*(d + e*x^3)^q*(1 + (e*x^3)/d)^(-p - q)*(d^2 - e^2*x^6)^p*AppellF1[1/3, 
-p, -p - q, 4/3, (e*x^3)/d, -((e*x^3)/d)])/(1 - (e*x^3)/d)^p
 

Defintions of rubi rules used

rule 936

Int[((a_) + (b_.)*(x_)^(n_))^(p_)*((c_) + (d_.)*(x_)^(n_))^(q_), x_Symbol] 
:> Simp[a^p*c^q*x*AppellF1[1/n, -p, -q, 1 + 1/n, (-b)*(x^n/a), (-d)*(x^n/c) 
], x] /; FreeQ[{a, b, c, d, n, p, q}, x] && NeQ[b*c - a*d, 0] && NeQ[n, -1] 
 && (IntegerQ[p] || GtQ[a, 0]) && (IntegerQ[q] || GtQ[c, 0])
 

rule 937

Int[((a_) + (b_.)*(x_)^(n_))^(p_)*((c_) + (d_.)*(x_)^(n_))^(q_), x_Symbol] 
:> Simp[a^IntPart[p]*((a + b*x^n)^FracPart[p]/(1 + b*(x^n/a))^FracPart[p]) 
  Int[(1 + b*(x^n/a))^p*(c + d*x^n)^q, x], x] /; FreeQ[{a, b, c, d, n, p, q 
}, x] && NeQ[b*c - a*d, 0] && NeQ[n, -1] &&  !(IntegerQ[p] || GtQ[a, 0])
 

rule 1396

Int[(u_.)*((a_) + (c_.)*(x_)^(n2_.))^(p_)*((d_) + (e_.)*(x_)^(n_))^(q_.), x 
_Symbol] :> Simp[(a + c*x^(2*n))^FracPart[p]/((d + e*x^n)^FracPart[p]*(a/d 
+ c*(x^n/e))^FracPart[p])   Int[u*(d + e*x^n)^(p + q)*(a/d + (c/e)*x^n)^p, 
x], x] /; FreeQ[{a, c, d, e, n, p, q}, x] && EqQ[n2, 2*n] && EqQ[c*d^2 + a* 
e^2, 0] &&  !IntegerQ[p] &&  !(EqQ[q, 1] && EqQ[n, 2])
 

Maple [F]

\[\int \left (e \,x^{3}+d \right )^{q} \left (-e^{2} x^{6}+d^{2}\right )^{p}d x\]

Input:

int((e*x^3+d)^q*(-e^2*x^6+d^2)^p,x)
 

Output:

int((e*x^3+d)^q*(-e^2*x^6+d^2)^p,x)
 

Fricas [F]

\[ \int \left (d+e x^3\right )^q \left (d^2-e^2 x^6\right )^p \, dx=\int { {\left (-e^{2} x^{6} + d^{2}\right )}^{p} {\left (e x^{3} + d\right )}^{q} \,d x } \] Input:

integrate((e*x^3+d)^q*(-e^2*x^6+d^2)^p,x, algorithm="fricas")
 

Output:

integral((-e^2*x^6 + d^2)^p*(e*x^3 + d)^q, x)
 

Sympy [F(-1)]

Timed out. \[ \int \left (d+e x^3\right )^q \left (d^2-e^2 x^6\right )^p \, dx=\text {Timed out} \] Input:

integrate((e*x**3+d)**q*(-e**2*x**6+d**2)**p,x)
 

Output:

Timed out
 

Maxima [F]

\[ \int \left (d+e x^3\right )^q \left (d^2-e^2 x^6\right )^p \, dx=\int { {\left (-e^{2} x^{6} + d^{2}\right )}^{p} {\left (e x^{3} + d\right )}^{q} \,d x } \] Input:

integrate((e*x^3+d)^q*(-e^2*x^6+d^2)^p,x, algorithm="maxima")
 

Output:

integrate((-e^2*x^6 + d^2)^p*(e*x^3 + d)^q, x)
 

Giac [F]

\[ \int \left (d+e x^3\right )^q \left (d^2-e^2 x^6\right )^p \, dx=\int { {\left (-e^{2} x^{6} + d^{2}\right )}^{p} {\left (e x^{3} + d\right )}^{q} \,d x } \] Input:

integrate((e*x^3+d)^q*(-e^2*x^6+d^2)^p,x, algorithm="giac")
 

Output:

integrate((-e^2*x^6 + d^2)^p*(e*x^3 + d)^q, x)
 

Mupad [F(-1)]

Timed out. \[ \int \left (d+e x^3\right )^q \left (d^2-e^2 x^6\right )^p \, dx=\int {\left (d^2-e^2\,x^6\right )}^p\,{\left (e\,x^3+d\right )}^q \,d x \] Input:

int((d^2 - e^2*x^6)^p*(d + e*x^3)^q,x)
 

Output:

int((d^2 - e^2*x^6)^p*(d + e*x^3)^q, x)
 

Reduce [F]

\[ \int \left (d+e x^3\right )^q \left (d^2-e^2 x^6\right )^p \, dx=\frac {\left (e \,x^{3}+d \right )^{q} \left (-e^{2} x^{6}+d^{2}\right )^{p} x -18 \left (\int \frac {\left (e \,x^{3}+d \right )^{q} \left (-e^{2} x^{6}+d^{2}\right )^{p} x^{3}}{-6 e^{2} p \,x^{6}-3 e^{2} q \,x^{6}-e^{2} x^{6}+6 d^{2} p +3 d^{2} q +d^{2}}d x \right ) d e p q -9 \left (\int \frac {\left (e \,x^{3}+d \right )^{q} \left (-e^{2} x^{6}+d^{2}\right )^{p} x^{3}}{-6 e^{2} p \,x^{6}-3 e^{2} q \,x^{6}-e^{2} x^{6}+6 d^{2} p +3 d^{2} q +d^{2}}d x \right ) d e \,q^{2}-3 \left (\int \frac {\left (e \,x^{3}+d \right )^{q} \left (-e^{2} x^{6}+d^{2}\right )^{p} x^{3}}{-6 e^{2} p \,x^{6}-3 e^{2} q \,x^{6}-e^{2} x^{6}+6 d^{2} p +3 d^{2} q +d^{2}}d x \right ) d e q +36 \left (\int \frac {\left (e \,x^{3}+d \right )^{q} \left (-e^{2} x^{6}+d^{2}\right )^{p}}{-6 e^{2} p \,x^{6}-3 e^{2} q \,x^{6}-e^{2} x^{6}+6 d^{2} p +3 d^{2} q +d^{2}}d x \right ) d^{2} p^{2}+36 \left (\int \frac {\left (e \,x^{3}+d \right )^{q} \left (-e^{2} x^{6}+d^{2}\right )^{p}}{-6 e^{2} p \,x^{6}-3 e^{2} q \,x^{6}-e^{2} x^{6}+6 d^{2} p +3 d^{2} q +d^{2}}d x \right ) d^{2} p q +6 \left (\int \frac {\left (e \,x^{3}+d \right )^{q} \left (-e^{2} x^{6}+d^{2}\right )^{p}}{-6 e^{2} p \,x^{6}-3 e^{2} q \,x^{6}-e^{2} x^{6}+6 d^{2} p +3 d^{2} q +d^{2}}d x \right ) d^{2} p +9 \left (\int \frac {\left (e \,x^{3}+d \right )^{q} \left (-e^{2} x^{6}+d^{2}\right )^{p}}{-6 e^{2} p \,x^{6}-3 e^{2} q \,x^{6}-e^{2} x^{6}+6 d^{2} p +3 d^{2} q +d^{2}}d x \right ) d^{2} q^{2}+3 \left (\int \frac {\left (e \,x^{3}+d \right )^{q} \left (-e^{2} x^{6}+d^{2}\right )^{p}}{-6 e^{2} p \,x^{6}-3 e^{2} q \,x^{6}-e^{2} x^{6}+6 d^{2} p +3 d^{2} q +d^{2}}d x \right ) d^{2} q}{6 p +3 q +1} \] Input:

int((e*x^3+d)^q*(-e^2*x^6+d^2)^p,x)
 

Output:

((d + e*x**3)**q*(d**2 - e**2*x**6)**p*x - 18*int(((d + e*x**3)**q*(d**2 - 
 e**2*x**6)**p*x**3)/(6*d**2*p + 3*d**2*q + d**2 - 6*e**2*p*x**6 - 3*e**2* 
q*x**6 - e**2*x**6),x)*d*e*p*q - 9*int(((d + e*x**3)**q*(d**2 - e**2*x**6) 
**p*x**3)/(6*d**2*p + 3*d**2*q + d**2 - 6*e**2*p*x**6 - 3*e**2*q*x**6 - e* 
*2*x**6),x)*d*e*q**2 - 3*int(((d + e*x**3)**q*(d**2 - e**2*x**6)**p*x**3)/ 
(6*d**2*p + 3*d**2*q + d**2 - 6*e**2*p*x**6 - 3*e**2*q*x**6 - e**2*x**6),x 
)*d*e*q + 36*int(((d + e*x**3)**q*(d**2 - e**2*x**6)**p)/(6*d**2*p + 3*d** 
2*q + d**2 - 6*e**2*p*x**6 - 3*e**2*q*x**6 - e**2*x**6),x)*d**2*p**2 + 36* 
int(((d + e*x**3)**q*(d**2 - e**2*x**6)**p)/(6*d**2*p + 3*d**2*q + d**2 - 
6*e**2*p*x**6 - 3*e**2*q*x**6 - e**2*x**6),x)*d**2*p*q + 6*int(((d + e*x** 
3)**q*(d**2 - e**2*x**6)**p)/(6*d**2*p + 3*d**2*q + d**2 - 6*e**2*p*x**6 - 
 3*e**2*q*x**6 - e**2*x**6),x)*d**2*p + 9*int(((d + e*x**3)**q*(d**2 - e** 
2*x**6)**p)/(6*d**2*p + 3*d**2*q + d**2 - 6*e**2*p*x**6 - 3*e**2*q*x**6 - 
e**2*x**6),x)*d**2*q**2 + 3*int(((d + e*x**3)**q*(d**2 - e**2*x**6)**p)/(6 
*d**2*p + 3*d**2*q + d**2 - 6*e**2*p*x**6 - 3*e**2*q*x**6 - e**2*x**6),x)* 
d**2*q)/(6*p + 3*q + 1)