\(\int (b d+2 c d x)^2 (a+b x+c x^2)^p \, dx\) [339]

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

Optimal result

Integrand size = 24, antiderivative size = 90 \[ \int (b d+2 c d x)^2 \left (a+b x+c x^2\right )^p \, dx=-\frac {2 d^2 (b+2 c x)^3 \left (\frac {1}{4} \left (4 a-\frac {b^2}{c}\right )+\frac {(b+2 c x)^2}{4 c}\right )^{1+p} \operatorname {Hypergeometric2F1}\left (1,\frac {5}{2}+p,\frac {5}{2},\frac {(b+2 c x)^2}{b^2-4 a c}\right )}{3 \left (b^2-4 a c\right )} \] Output:

-2*d^2*(2*c*x+b)^3*(a-1/4*b^2/c+1/4*(2*c*x+b)^2/c)^(p+1)*hypergeom([1, 5/2 
+p],[5/2],(2*c*x+b)^2/(-4*a*c+b^2))/(-12*a*c+3*b^2)
 

Mathematica [A] (verified)

Time = 0.15 (sec) , antiderivative size = 92, normalized size of antiderivative = 1.02 \[ \int (b d+2 c d x)^2 \left (a+b x+c x^2\right )^p \, dx=\frac {2^{-1-2 p} d^2 (b+2 c x)^3 (a+x (b+c x))^p \left (\frac {c (a+x (b+c x))}{-b^2+4 a c}\right )^{-p} \operatorname {Hypergeometric2F1}\left (\frac {3}{2},-p,\frac {5}{2},\frac {(b+2 c x)^2}{b^2-4 a c}\right )}{3 c} \] Input:

Integrate[(b*d + 2*c*d*x)^2*(a + b*x + c*x^2)^p,x]
 

Output:

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

Rubi [A] (verified)

Time = 0.22 (sec) , antiderivative size = 128, normalized size of antiderivative = 1.42, number of steps used = 4, number of rules used = 3, \(\frac {\text {number of rules}}{\text {integrand size}}\) = 0.125, Rules used = {1118, 279, 278}

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

\(\Big \downarrow \) 1118

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

\(\Big \downarrow \) 279

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

\(\Big \downarrow \) 278

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

Input:

Int[(b*d + 2*c*d*x)^2*(a + b*x + c*x^2)^p,x]
 

Output:

((b*d + 2*c*d*x)^3*((4*a - b^2/c)/4 + (b*d + 2*c*d*x)^2/(4*c*d^2))^p*Hyper 
geometric2F1[3/2, -p, 5/2, (b*d + 2*c*d*x)^2/((b^2 - 4*a*c)*d^2)])/(6*c*d* 
(1 - (b*d + 2*c*d*x)^2/((b^2 - 4*a*c)*d^2))^p)
 

Defintions of rubi rules used

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

rule 279
Int[((c_.)*(x_))^(m_.)*((a_) + (b_.)*(x_)^2)^(p_), x_Symbol] :> Simp[a^IntP 
art[p]*((a + b*x^2)^FracPart[p]/(1 + b*(x^2/a))^FracPart[p])   Int[(c*x)^m* 
(1 + b*(x^2/a))^p, x], x] /; FreeQ[{a, b, c, m, p}, x] &&  !IGtQ[p, 0] && 
!(ILtQ[p, 0] || GtQ[a, 0])
 

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

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

Input:

int((2*c*d*x+b*d)^2*(c*x^2+b*x+a)^p,x)
 

Output:

int((2*c*d*x+b*d)^2*(c*x^2+b*x+a)^p,x)
 

Fricas [F]

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

integrate((2*c*d*x+b*d)^2*(c*x^2+b*x+a)^p,x, algorithm="fricas")
 

Output:

integral((4*c^2*d^2*x^2 + 4*b*c*d^2*x + b^2*d^2)*(c*x^2 + b*x + a)^p, x)
 

Sympy [F]

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

integrate((2*c*d*x+b*d)**2*(c*x**2+b*x+a)**p,x)
 

Output:

d**2*(Integral(b**2*(a + b*x + c*x**2)**p, x) + Integral(4*c**2*x**2*(a + 
b*x + c*x**2)**p, x) + Integral(4*b*c*x*(a + b*x + c*x**2)**p, x))
 

Maxima [F]

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

integrate((2*c*d*x+b*d)^2*(c*x^2+b*x+a)^p,x, algorithm="maxima")
 

Output:

integrate((2*c*d*x + b*d)^2*(c*x^2 + b*x + a)^p, x)
 

Giac [F]

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

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

Output:

integrate((2*c*d*x + b*d)^2*(c*x^2 + b*x + a)^p, x)
 

Mupad [F(-1)]

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

int((b*d + 2*c*d*x)^2*(a + b*x + c*x^2)^p,x)
 

Output:

int((b*d + 2*c*d*x)^2*(a + b*x + c*x^2)^p, x)
 

Reduce [F]

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

int((2*c*d*x+b*d)^2*(c*x^2+b*x+a)^p,x)
 

Output:

(d**2*( - 8*(a + b*x + c*x**2)**p*a**2*c + 4*(a + b*x + c*x**2)**p*a*b**2* 
p + 4*(a + b*x + c*x**2)**p*a*b**2 + 8*(a + b*x + c*x**2)**p*a*b*c*p*x + 4 
*(a + b*x + c*x**2)**p*b**3*p*x + 3*(a + b*x + c*x**2)**p*b**3*x + 12*(a + 
 b*x + c*x**2)**p*b**2*c*p*x**2 + 6*(a + b*x + c*x**2)**p*b**2*c*x**2 + 8* 
(a + b*x + c*x**2)**p*b*c**2*p*x**3 + 4*(a + b*x + c*x**2)**p*b*c**2*x**3 
+ 64*int(((a + b*x + c*x**2)**p*x)/(4*a*p**2 + 8*a*p + 3*a + 4*b*p**2*x + 
8*b*p*x + 3*b*x + 4*c*p**2*x**2 + 8*c*p*x**2 + 3*c*x**2),x)*a**2*c**2*p**3 
 + 128*int(((a + b*x + c*x**2)**p*x)/(4*a*p**2 + 8*a*p + 3*a + 4*b*p**2*x 
+ 8*b*p*x + 3*b*x + 4*c*p**2*x**2 + 8*c*p*x**2 + 3*c*x**2),x)*a**2*c**2*p* 
*2 + 48*int(((a + b*x + c*x**2)**p*x)/(4*a*p**2 + 8*a*p + 3*a + 4*b*p**2*x 
 + 8*b*p*x + 3*b*x + 4*c*p**2*x**2 + 8*c*p*x**2 + 3*c*x**2),x)*a**2*c**2*p 
 - 32*int(((a + b*x + c*x**2)**p*x)/(4*a*p**2 + 8*a*p + 3*a + 4*b*p**2*x + 
 8*b*p*x + 3*b*x + 4*c*p**2*x**2 + 8*c*p*x**2 + 3*c*x**2),x)*a*b**2*c*p**3 
 - 64*int(((a + b*x + c*x**2)**p*x)/(4*a*p**2 + 8*a*p + 3*a + 4*b*p**2*x + 
 8*b*p*x + 3*b*x + 4*c*p**2*x**2 + 8*c*p*x**2 + 3*c*x**2),x)*a*b**2*c*p**2 
 - 24*int(((a + b*x + c*x**2)**p*x)/(4*a*p**2 + 8*a*p + 3*a + 4*b*p**2*x + 
 8*b*p*x + 3*b*x + 4*c*p**2*x**2 + 8*c*p*x**2 + 3*c*x**2),x)*a*b**2*c*p + 
4*int(((a + b*x + c*x**2)**p*x)/(4*a*p**2 + 8*a*p + 3*a + 4*b*p**2*x + 8*b 
*p*x + 3*b*x + 4*c*p**2*x**2 + 8*c*p*x**2 + 3*c*x**2),x)*b**4*p**3 + 8*int 
(((a + b*x + c*x**2)**p*x)/(4*a*p**2 + 8*a*p + 3*a + 4*b*p**2*x + 8*b*p...