\(\int (A+B x) (a+b x+c x^2)^p \, dx\) [243]

Optimal result
Mathematica [C] (warning: unable to verify)
Rubi [A] (verified)
Maple [F]
Fricas [F]
Sympy [F]
Maxima [F]
Giac [F]
Mupad [F(-1)]
Reduce [F]

Optimal result

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

1/2*B*(c*x^2+b*x+a)^(p+1)/c/(p+1)-(-2*A*c+B*b)*(2*c*x+b)*(c*x^2+b*x+a)^p*h 
ypergeom([1/2, -p],[3/2],(2*c*x+b)^2/(-4*a*c+b^2))/(2^(2*p+2))/c^2/((-c*(c 
*x^2+b*x+a)/(-4*a*c+b^2))^p)
 

Mathematica [C] (warning: unable to verify)

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

Time = 0.66 (sec) , antiderivative size = 268, normalized size of antiderivative = 2.20 \[ \int (A+B x) \left (a+b x+c x^2\right )^p \, dx=\frac {1}{2} (a+x (b+c x))^p \left (B x^2 \left (\frac {b-\sqrt {b^2-4 a c}+2 c x}{b-\sqrt {b^2-4 a c}}\right )^{-p} \left (\frac {b+\sqrt {b^2-4 a c}+2 c x}{b+\sqrt {b^2-4 a c}}\right )^{-p} \operatorname {AppellF1}\left (2,-p,-p,3,-\frac {2 c x}{b+\sqrt {b^2-4 a c}},\frac {2 c x}{-b+\sqrt {b^2-4 a c}}\right )+\frac {2^p A \left (b-\sqrt {b^2-4 a c}+2 c x\right ) \left (\frac {b+\sqrt {b^2-4 a c}+2 c x}{\sqrt {b^2-4 a c}}\right )^{-p} \operatorname {Hypergeometric2F1}\left (-p,1+p,2+p,\frac {-b+\sqrt {b^2-4 a c}-2 c x}{2 \sqrt {b^2-4 a c}}\right )}{c (1+p)}\right ) \] Input:

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

Output:

((a + x*(b + c*x))^p*((B*x^2*AppellF1[2, -p, -p, 3, (-2*c*x)/(b + Sqrt[b^2 
 - 4*a*c]), (2*c*x)/(-b + Sqrt[b^2 - 4*a*c])])/(((b - Sqrt[b^2 - 4*a*c] + 
2*c*x)/(b - Sqrt[b^2 - 4*a*c]))^p*((b + Sqrt[b^2 - 4*a*c] + 2*c*x)/(b + Sq 
rt[b^2 - 4*a*c]))^p) + (2^p*A*(b - Sqrt[b^2 - 4*a*c] + 2*c*x)*Hypergeometr 
ic2F1[-p, 1 + p, 2 + p, (-b + Sqrt[b^2 - 4*a*c] - 2*c*x)/(2*Sqrt[b^2 - 4*a 
*c])])/(c*(1 + p)*((b + Sqrt[b^2 - 4*a*c] + 2*c*x)/Sqrt[b^2 - 4*a*c])^p))) 
/2
 

Rubi [A] (verified)

Time = 0.27 (sec) , antiderivative size = 158, normalized size of antiderivative = 1.30, number of steps used = 2, number of rules used = 2, \(\frac {\text {number of rules}}{\text {integrand size}}\) = 0.111, Rules used = {1160, 1096}

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

\(\Big \downarrow \) 1160

\(\displaystyle \frac {B \left (a+b x+c x^2\right )^{p+1}}{2 c (p+1)}-\frac {(b B-2 A c) \int \left (c x^2+b x+a\right )^pdx}{2 c}\)

\(\Big \downarrow \) 1096

\(\displaystyle \frac {2^p (b B-2 A c) \left (a+b x+c x^2\right )^{p+1} \left (-\frac {-\sqrt {b^2-4 a c}+b+2 c x}{\sqrt {b^2-4 a c}}\right )^{-p-1} \operatorname {Hypergeometric2F1}\left (-p,p+1,p+2,\frac {b+2 c x+\sqrt {b^2-4 a c}}{2 \sqrt {b^2-4 a c}}\right )}{c (p+1) \sqrt {b^2-4 a c}}+\frac {B \left (a+b x+c x^2\right )^{p+1}}{2 c (p+1)}\)

Input:

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

Output:

(B*(a + b*x + c*x^2)^(1 + p))/(2*c*(1 + p)) + (2^p*(b*B - 2*A*c)*(-((b - S 
qrt[b^2 - 4*a*c] + 2*c*x)/Sqrt[b^2 - 4*a*c]))^(-1 - p)*(a + b*x + c*x^2)^( 
1 + p)*Hypergeometric2F1[-p, 1 + p, 2 + p, (b + Sqrt[b^2 - 4*a*c] + 2*c*x) 
/(2*Sqrt[b^2 - 4*a*c])])/(c*Sqrt[b^2 - 4*a*c]*(1 + p))
 

Defintions of rubi rules used

rule 1096
Int[((a_.) + (b_.)*(x_) + (c_.)*(x_)^2)^(p_), x_Symbol] :> With[{q = Rt[b^2 
 - 4*a*c, 2]}, Simp[(-(a + b*x + c*x^2)^(p + 1)/(q*(p + 1)*((q - b - 2*c*x) 
/(2*q))^(p + 1)))*Hypergeometric2F1[-p, p + 1, p + 2, (b + q + 2*c*x)/(2*q) 
], x]] /; FreeQ[{a, b, c, p}, x] &&  !IntegerQ[4*p] &&  !IntegerQ[3*p]
 

rule 1160
Int[((d_.) + (e_.)*(x_))*((a_.) + (b_.)*(x_) + (c_.)*(x_)^2)^(p_), x_Symbol 
] :> Simp[e*((a + b*x + c*x^2)^(p + 1)/(2*c*(p + 1))), x] + Simp[(2*c*d - b 
*e)/(2*c)   Int[(a + b*x + c*x^2)^p, x], x] /; FreeQ[{a, b, c, d, e, p}, x] 
 && NeQ[p, -1]
 
Maple [F]

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

Input:

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

Output:

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

Fricas [F]

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

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

Output:

integral((B*x + A)*(c*x^2 + b*x + a)^p, x)
 

Sympy [F]

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

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

Output:

Integral((A + B*x)*(a + b*x + c*x**2)**p, x)
 

Maxima [F]

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

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

Output:

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

Giac [F]

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

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

Output:

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

Mupad [F(-1)]

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

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

Output:

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

Reduce [F]

\[ \int (A+B x) \left (a+b x+c x^2\right )^p \, dx=\frac {4 \left (c \,x^{2}+b x +a \right )^{p} a^{2} c p +4 \left (c \,x^{2}+b x +a \right )^{p} a^{2} c -\left (c \,x^{2}+b x +a \right )^{p} a \,b^{2}+2 \left (c \,x^{2}+b x +a \right )^{p} a b c p x +2 \left (c \,x^{2}+b x +a \right )^{p} a b c x +\left (c \,x^{2}+b x +a \right )^{p} b^{3} p x +2 \left (c \,x^{2}+b x +a \right )^{p} b^{2} c p \,x^{2}+\left (c \,x^{2}+b x +a \right )^{p} b^{2} c \,x^{2}-16 \left (\int \frac {\left (c \,x^{2}+b x +a \right )^{p} x}{2 c p \,x^{2}+2 b p x +c \,x^{2}+2 a p +b x +a}d x \right ) a^{2} c^{2} p^{3}-24 \left (\int \frac {\left (c \,x^{2}+b x +a \right )^{p} x}{2 c p \,x^{2}+2 b p x +c \,x^{2}+2 a p +b x +a}d x \right ) a^{2} c^{2} p^{2}-8 \left (\int \frac {\left (c \,x^{2}+b x +a \right )^{p} x}{2 c p \,x^{2}+2 b p x +c \,x^{2}+2 a p +b x +a}d x \right ) a^{2} c^{2} p +12 \left (\int \frac {\left (c \,x^{2}+b x +a \right )^{p} x}{2 c p \,x^{2}+2 b p x +c \,x^{2}+2 a p +b x +a}d x \right ) a \,b^{2} c \,p^{3}+18 \left (\int \frac {\left (c \,x^{2}+b x +a \right )^{p} x}{2 c p \,x^{2}+2 b p x +c \,x^{2}+2 a p +b x +a}d x \right ) a \,b^{2} c \,p^{2}+6 \left (\int \frac {\left (c \,x^{2}+b x +a \right )^{p} x}{2 c p \,x^{2}+2 b p x +c \,x^{2}+2 a p +b x +a}d x \right ) a \,b^{2} c p -2 \left (\int \frac {\left (c \,x^{2}+b x +a \right )^{p} x}{2 c p \,x^{2}+2 b p x +c \,x^{2}+2 a p +b x +a}d x \right ) b^{4} p^{3}-3 \left (\int \frac {\left (c \,x^{2}+b x +a \right )^{p} x}{2 c p \,x^{2}+2 b p x +c \,x^{2}+2 a p +b x +a}d x \right ) b^{4} p^{2}-\left (\int \frac {\left (c \,x^{2}+b x +a \right )^{p} x}{2 c p \,x^{2}+2 b p x +c \,x^{2}+2 a p +b x +a}d x \right ) b^{4} p}{2 b c \left (2 p^{2}+3 p +1\right )} \] Input:

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

Output:

(4*(a + b*x + c*x**2)**p*a**2*c*p + 4*(a + b*x + c*x**2)**p*a**2*c - (a + 
b*x + c*x**2)**p*a*b**2 + 2*(a + b*x + c*x**2)**p*a*b*c*p*x + 2*(a + b*x + 
 c*x**2)**p*a*b*c*x + (a + b*x + c*x**2)**p*b**3*p*x + 2*(a + b*x + c*x**2 
)**p*b**2*c*p*x**2 + (a + b*x + c*x**2)**p*b**2*c*x**2 - 16*int(((a + b*x 
+ c*x**2)**p*x)/(2*a*p + a + 2*b*p*x + b*x + 2*c*p*x**2 + c*x**2),x)*a**2* 
c**2*p**3 - 24*int(((a + b*x + c*x**2)**p*x)/(2*a*p + a + 2*b*p*x + b*x + 
2*c*p*x**2 + c*x**2),x)*a**2*c**2*p**2 - 8*int(((a + b*x + c*x**2)**p*x)/( 
2*a*p + a + 2*b*p*x + b*x + 2*c*p*x**2 + c*x**2),x)*a**2*c**2*p + 12*int(( 
(a + b*x + c*x**2)**p*x)/(2*a*p + a + 2*b*p*x + b*x + 2*c*p*x**2 + c*x**2) 
,x)*a*b**2*c*p**3 + 18*int(((a + b*x + c*x**2)**p*x)/(2*a*p + a + 2*b*p*x 
+ b*x + 2*c*p*x**2 + c*x**2),x)*a*b**2*c*p**2 + 6*int(((a + b*x + c*x**2)* 
*p*x)/(2*a*p + a + 2*b*p*x + b*x + 2*c*p*x**2 + c*x**2),x)*a*b**2*c*p - 2* 
int(((a + b*x + c*x**2)**p*x)/(2*a*p + a + 2*b*p*x + b*x + 2*c*p*x**2 + c* 
x**2),x)*b**4*p**3 - 3*int(((a + b*x + c*x**2)**p*x)/(2*a*p + a + 2*b*p*x 
+ b*x + 2*c*p*x**2 + c*x**2),x)*b**4*p**2 - int(((a + b*x + c*x**2)**p*x)/ 
(2*a*p + a + 2*b*p*x + b*x + 2*c*p*x**2 + c*x**2),x)*b**4*p)/(2*b*c*(2*p** 
2 + 3*p + 1))