\(\int \frac {(a d e+(c d^2+a e^2) x+c d e x^2)^p}{(d+e x)^3} \, dx\) [395]

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 = 35, antiderivative size = 90 \[ \int \frac {\left (a d e+\left (c d^2+a e^2\right ) x+c d e x^2\right )^p}{(d+e x)^3} \, dx=\frac {\left (a d e+\left (c d^2+a e^2\right ) x+c d e x^2\right )^{1+p} \operatorname {Hypergeometric2F1}\left (1,-1+2 p,-1+p,\frac {c d (d+e x)}{c d^2-a e^2}\right )}{\left (c d^2-a e^2\right ) (2-p) (d+e x)^3} \] Output:

(a*d*e+(a*e^2+c*d^2)*x+c*d*e*x^2)^(p+1)*hypergeom([1, -1+2*p],[-1+p],c*d*( 
e*x+d)/(-a*e^2+c*d^2))/(-a*e^2+c*d^2)/(2-p)/(e*x+d)^3
 

Mathematica [A] (verified)

Time = 0.08 (sec) , antiderivative size = 112, normalized size of antiderivative = 1.24 \[ \int \frac {\left (a d e+\left (c d^2+a e^2\right ) x+c d e x^2\right )^p}{(d+e x)^3} \, dx=\frac {c^2 d^2 (a e+c d x) \left (\frac {c d (d+e x)}{c d^2-a e^2}\right )^{-p} ((a e+c d x) (d+e x))^p \operatorname {Hypergeometric2F1}\left (3-p,1+p,2+p,\frac {e (a e+c d x)}{-c d^2+a e^2}\right )}{\left (c d^2-a e^2\right )^3 (1+p)} \] Input:

Integrate[(a*d*e + (c*d^2 + a*e^2)*x + c*d*e*x^2)^p/(d + e*x)^3,x]
 

Output:

(c^2*d^2*(a*e + c*d*x)*((a*e + c*d*x)*(d + e*x))^p*Hypergeometric2F1[3 - p 
, 1 + p, 2 + p, (e*(a*e + c*d*x))/(-(c*d^2) + a*e^2)])/((c*d^2 - a*e^2)^3* 
(1 + p)*((c*d*(d + e*x))/(c*d^2 - a*e^2))^p)
 

Rubi [A] (verified)

Time = 0.27 (sec) , antiderivative size = 124, normalized size of antiderivative = 1.38, number of steps used = 3, number of rules used = 3, \(\frac {\text {number of rules}}{\text {integrand size}}\) = 0.086, Rules used = {1138, 80, 79}

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 (x \left (a e^2+c d^2\right )+a d e+c d e x^2\right )^p}{(d+e x)^3} \, dx\)

\(\Big \downarrow \) 1138

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

\(\Big \downarrow \) 80

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

\(\Big \downarrow \) 79

\(\displaystyle \frac {c^2 d^2 (a e+c d x) \left (\frac {c d (d+e x)}{c d^2-a e^2}\right )^{-p} \left (x \left (a e^2+c d^2\right )+a d e+c d e x^2\right )^p \operatorname {Hypergeometric2F1}\left (3-p,p+1,p+2,-\frac {e (a e+c d x)}{c d^2-a e^2}\right )}{(p+1) \left (c d^2-a e^2\right )^3}\)

Input:

Int[(a*d*e + (c*d^2 + a*e^2)*x + c*d*e*x^2)^p/(d + e*x)^3,x]
 

Output:

(c^2*d^2*(a*e + c*d*x)*(a*d*e + (c*d^2 + a*e^2)*x + c*d*e*x^2)^p*Hypergeom 
etric2F1[3 - p, 1 + p, 2 + p, -((e*(a*e + c*d*x))/(c*d^2 - a*e^2))])/((c*d 
^2 - a*e^2)^3*(1 + p)*((c*d*(d + e*x))/(c*d^2 - a*e^2))^p)
 

Defintions of rubi rules used

rule 79
Int[((a_) + (b_.)*(x_))^(m_)*((c_) + (d_.)*(x_))^(n_), x_Symbol] :> Simp[(( 
a + b*x)^(m + 1)/(b*(m + 1)*(b/(b*c - a*d))^n))*Hypergeometric2F1[-n, m + 1 
, m + 2, (-d)*((a + b*x)/(b*c - a*d))], x] /; FreeQ[{a, b, c, d, m, n}, x] 
&&  !IntegerQ[m] &&  !IntegerQ[n] && GtQ[b/(b*c - a*d), 0] && (RationalQ[m] 
 ||  !(RationalQ[n] && GtQ[-d/(b*c - a*d), 0]))
 

rule 80
Int[((a_) + (b_.)*(x_))^(m_)*((c_) + (d_.)*(x_))^(n_), x_Symbol] :> Simp[(c 
 + d*x)^FracPart[n]/((b/(b*c - a*d))^IntPart[n]*(b*((c + d*x)/(b*c - a*d))) 
^FracPart[n])   Int[(a + b*x)^m*Simp[b*(c/(b*c - a*d)) + b*d*(x/(b*c - a*d) 
), x]^n, x], x] /; FreeQ[{a, b, c, d, m, n}, x] &&  !IntegerQ[m] &&  !Integ 
erQ[n] && (RationalQ[m] ||  !SimplerQ[n + 1, m + 1])
 

rule 1138
Int[((d_) + (e_.)*(x_))^(m_)*((a_.) + (b_.)*(x_) + (c_.)*(x_)^2)^(p_), x_Sy 
mbol] :> Simp[d^m*((a + b*x + c*x^2)^FracPart[p]/((1 + e*(x/d))^FracPart[p] 
*(a/d + (c*x)/e)^FracPart[p]))   Int[(1 + e*(x/d))^(m + p)*(a/d + (c/e)*x)^ 
p, x], x] /; FreeQ[{a, b, c, d, e, m}, x] && EqQ[c*d^2 - b*d*e + a*e^2, 0] 
&& (IntegerQ[m] || GtQ[d, 0]) &&  !(IGtQ[m, 0] && (IntegerQ[3*p] || Integer 
Q[4*p]))
 
Maple [F]

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

Input:

int((a*d*e+(a*e^2+c*d^2)*x+c*d*x^2*e)^p/(e*x+d)^3,x)
 

Output:

int((a*d*e+(a*e^2+c*d^2)*x+c*d*x^2*e)^p/(e*x+d)^3,x)
 

Fricas [F]

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

integrate((a*d*e+(a*e^2+c*d^2)*x+c*d*e*x^2)^p/(e*x+d)^3,x, algorithm="fric 
as")
 

Output:

integral((c*d*e*x^2 + a*d*e + (c*d^2 + a*e^2)*x)^p/(e^3*x^3 + 3*d*e^2*x^2 
+ 3*d^2*e*x + d^3), x)
 

Sympy [F(-1)]

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

integrate((a*d*e+(a*e**2+c*d**2)*x+c*d*e*x**2)**p/(e*x+d)**3,x)
 

Output:

Timed out
 

Maxima [F]

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

integrate((a*d*e+(a*e^2+c*d^2)*x+c*d*e*x^2)^p/(e*x+d)^3,x, algorithm="maxi 
ma")
 

Output:

integrate((c*d*e*x^2 + a*d*e + (c*d^2 + a*e^2)*x)^p/(e*x + d)^3, x)
 

Giac [F]

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

integrate((a*d*e+(a*e^2+c*d^2)*x+c*d*e*x^2)^p/(e*x+d)^3,x, algorithm="giac 
")
 

Output:

integrate((c*d*e*x^2 + a*d*e + (c*d^2 + a*e^2)*x)^p/(e*x + d)^3, x)
 

Mupad [F(-1)]

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

int((x*(a*e^2 + c*d^2) + a*d*e + c*d*e*x^2)^p/(d + e*x)^3,x)
 

Output:

int((x*(a*e^2 + c*d^2) + a*d*e + c*d*e*x^2)^p/(d + e*x)^3, x)
 

Reduce [F]

\[ \int \frac {\left (a d e+\left (c d^2+a e^2\right ) x+c d e x^2\right )^p}{(d+e x)^3} \, dx=\text {too large to display} \] Input:

int((a*d*e+(a*e^2+c*d^2)*x+c*d*e*x^2)^p/(e*x+d)^3,x)
 

Output:

((a*d*e + a*e**2*x + c*d**2*x + c*d*e*x**2)**p*a*e - int(((a*d*e + a*e**2* 
x + c*d**2*x + c*d*e*x**2)**p*x)/(a**2*d**3*e**3*p - 2*a**2*d**3*e**3 + 3* 
a**2*d**2*e**4*p*x - 6*a**2*d**2*e**4*x + 3*a**2*d*e**5*p*x**2 - 6*a**2*d* 
e**5*x**2 + a**2*e**6*p*x**3 - 2*a**2*e**6*x**3 + a*c*d**5*e*p + 4*a*c*d** 
4*e**2*p*x - 2*a*c*d**4*e**2*x + 6*a*c*d**3*e**3*p*x**2 - 6*a*c*d**3*e**3* 
x**2 + 4*a*c*d**2*e**4*p*x**3 - 6*a*c*d**2*e**4*x**3 + a*c*d*e**5*p*x**4 - 
 2*a*c*d*e**5*x**4 + c**2*d**6*p*x + 3*c**2*d**5*e*p*x**2 + 3*c**2*d**4*e* 
*2*p*x**3 + c**2*d**3*e**3*p*x**4),x)*a**2*c*d**3*e**4*p**2 + 2*int(((a*d* 
e + a*e**2*x + c*d**2*x + c*d*e*x**2)**p*x)/(a**2*d**3*e**3*p - 2*a**2*d** 
3*e**3 + 3*a**2*d**2*e**4*p*x - 6*a**2*d**2*e**4*x + 3*a**2*d*e**5*p*x**2 
- 6*a**2*d*e**5*x**2 + a**2*e**6*p*x**3 - 2*a**2*e**6*x**3 + a*c*d**5*e*p 
+ 4*a*c*d**4*e**2*p*x - 2*a*c*d**4*e**2*x + 6*a*c*d**3*e**3*p*x**2 - 6*a*c 
*d**3*e**3*x**2 + 4*a*c*d**2*e**4*p*x**3 - 6*a*c*d**2*e**4*x**3 + a*c*d*e* 
*5*p*x**4 - 2*a*c*d*e**5*x**4 + c**2*d**6*p*x + 3*c**2*d**5*e*p*x**2 + 3*c 
**2*d**4*e**2*p*x**3 + c**2*d**3*e**3*p*x**4),x)*a**2*c*d**3*e**4*p - 2*in 
t(((a*d*e + a*e**2*x + c*d**2*x + c*d*e*x**2)**p*x)/(a**2*d**3*e**3*p - 2* 
a**2*d**3*e**3 + 3*a**2*d**2*e**4*p*x - 6*a**2*d**2*e**4*x + 3*a**2*d*e**5 
*p*x**2 - 6*a**2*d*e**5*x**2 + a**2*e**6*p*x**3 - 2*a**2*e**6*x**3 + a*c*d 
**5*e*p + 4*a*c*d**4*e**2*p*x - 2*a*c*d**4*e**2*x + 6*a*c*d**3*e**3*p*x**2 
 - 6*a*c*d**3*e**3*x**2 + 4*a*c*d**2*e**4*p*x**3 - 6*a*c*d**2*e**4*x**3...