\(\int (a+b x) (A+B x) (d+e x)^m \, dx\) [232]

Optimal result
Mathematica [A] (verified)
Rubi [A] (verified)
Maple [B] (verified)
Fricas [B] (verification not implemented)
Sympy [B] (verification not implemented)
Maxima [A] (verification not implemented)
Giac [B] (verification not implemented)
Mupad [B] (verification not implemented)
Reduce [B] (verification not implemented)

Optimal result

Integrand size = 18, antiderivative size = 90 \[ \int (a+b x) (A+B x) (d+e x)^m \, dx=\frac {(b d-a e) (B d-A e) (d+e x)^{1+m}}{e^3 (1+m)}-\frac {(2 b B d-A b e-a B e) (d+e x)^{2+m}}{e^3 (2+m)}+\frac {b B (d+e x)^{3+m}}{e^3 (3+m)} \] Output:

(-a*e+b*d)*(-A*e+B*d)*(e*x+d)^(1+m)/e^3/(1+m)-(-A*b*e-B*a*e+2*B*b*d)*(e*x+ 
d)^(2+m)/e^3/(2+m)+b*B*(e*x+d)^(3+m)/e^3/(3+m)
 

Mathematica [A] (verified)

Time = 0.10 (sec) , antiderivative size = 79, normalized size of antiderivative = 0.88 \[ \int (a+b x) (A+B x) (d+e x)^m \, dx=\frac {(d+e x)^{1+m} \left (\frac {(b d-a e) (B d-A e)}{1+m}-\frac {(2 b B d-A b e-a B e) (d+e x)}{2+m}+\frac {b B (d+e x)^2}{3+m}\right )}{e^3} \] Input:

Integrate[(a + b*x)*(A + B*x)*(d + e*x)^m,x]
 

Output:

((d + e*x)^(1 + m)*(((b*d - a*e)*(B*d - A*e))/(1 + m) - ((2*b*B*d - A*b*e 
- a*B*e)*(d + e*x))/(2 + m) + (b*B*(d + e*x)^2)/(3 + m)))/e^3
 

Rubi [A] (verified)

Time = 0.24 (sec) , antiderivative size = 90, normalized size of antiderivative = 1.00, number of steps used = 2, number of rules used = 2, \(\frac {\text {number of rules}}{\text {integrand size}}\) = 0.111, Rules used = {86, 2009}

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) (A+B x) (d+e x)^m \, dx\)

\(\Big \downarrow \) 86

\(\displaystyle \int \left (\frac {(a e-b d) (A e-B d) (d+e x)^m}{e^2}+\frac {(d+e x)^{m+1} (a B e+A b e-2 b B d)}{e^2}+\frac {b B (d+e x)^{m+2}}{e^2}\right )dx\)

\(\Big \downarrow \) 2009

\(\displaystyle \frac {(b d-a e) (B d-A e) (d+e x)^{m+1}}{e^3 (m+1)}-\frac {(d+e x)^{m+2} (-a B e-A b e+2 b B d)}{e^3 (m+2)}+\frac {b B (d+e x)^{m+3}}{e^3 (m+3)}\)

Input:

Int[(a + b*x)*(A + B*x)*(d + e*x)^m,x]
 

Output:

((b*d - a*e)*(B*d - A*e)*(d + e*x)^(1 + m))/(e^3*(1 + m)) - ((2*b*B*d - A* 
b*e - a*B*e)*(d + e*x)^(2 + m))/(e^3*(2 + m)) + (b*B*(d + e*x)^(3 + m))/(e 
^3*(3 + m))
 

Defintions of rubi rules used

rule 86
Int[((a_.) + (b_.)*(x_))*((c_) + (d_.)*(x_))^(n_.)*((e_.) + (f_.)*(x_))^(p_ 
.), x_] :> Int[ExpandIntegrand[(a + b*x)*(c + d*x)^n*(e + f*x)^p, x], x] /; 
 FreeQ[{a, b, c, d, e, f, n}, x] && ((ILtQ[n, 0] && ILtQ[p, 0]) || EqQ[p, 1 
] || (IGtQ[p, 0] && ( !IntegerQ[n] || LeQ[9*p + 5*(n + 2), 0] || GeQ[n + p 
+ 1, 0] || (GeQ[n + p + 2, 0] && RationalQ[a, b, c, d, e, f]))))
 

rule 2009
Int[u_, x_Symbol] :> Simp[IntSum[u, x], x] /; SumQ[u]
 
Maple [B] (verified)

Leaf count of result is larger than twice the leaf count of optimal. \(188\) vs. \(2(90)=180\).

Time = 0.23 (sec) , antiderivative size = 189, normalized size of antiderivative = 2.10

method result size
gosper \(\frac {\left (e x +d \right )^{1+m} \left (B b \,e^{2} m^{2} x^{2}+A b \,e^{2} m^{2} x +B a \,e^{2} m^{2} x +3 B b \,e^{2} m \,x^{2}+A a \,e^{2} m^{2}+4 A b \,e^{2} m x +4 B a \,e^{2} m x -2 B b d e m x +2 b B \,x^{2} e^{2}+5 A a \,e^{2} m -A b d e m +3 A x b \,e^{2}-B a d e m +3 B x a \,e^{2}-2 B x b d e +6 A a \,e^{2}-3 A b d e -3 B a d e +2 b B \,d^{2}\right )}{e^{3} \left (m^{3}+6 m^{2}+11 m +6\right )}\) \(189\)
orering \(\frac {\left (e x +d \right ) \left (B b \,e^{2} m^{2} x^{2}+A b \,e^{2} m^{2} x +B a \,e^{2} m^{2} x +3 B b \,e^{2} m \,x^{2}+A a \,e^{2} m^{2}+4 A b \,e^{2} m x +4 B a \,e^{2} m x -2 B b d e m x +2 b B \,x^{2} e^{2}+5 A a \,e^{2} m -A b d e m +3 A x b \,e^{2}-B a d e m +3 B x a \,e^{2}-2 B x b d e +6 A a \,e^{2}-3 A b d e -3 B a d e +2 b B \,d^{2}\right ) \left (e x +d \right )^{m}}{e^{3} \left (m^{3}+6 m^{2}+11 m +6\right )}\) \(192\)
norman \(\frac {b B \,x^{3} {\mathrm e}^{m \ln \left (e x +d \right )}}{3+m}+\frac {d \left (A a \,e^{2} m^{2}+5 A a \,e^{2} m -A b d e m -B a d e m +6 A a \,e^{2}-3 A b d e -3 B a d e +2 b B \,d^{2}\right ) {\mathrm e}^{m \ln \left (e x +d \right )}}{e^{3} \left (m^{3}+6 m^{2}+11 m +6\right )}+\frac {\left (A b e m +B a e m +B b d m +3 A b e +3 B a e \right ) x^{2} {\mathrm e}^{m \ln \left (e x +d \right )}}{e \left (m^{2}+5 m +6\right )}+\frac {\left (A a \,e^{2} m^{2}+A b d e \,m^{2}+B a d e \,m^{2}+5 A a \,e^{2} m +3 A b d e m +3 B a d e m -2 B b \,d^{2} m +6 A a \,e^{2}\right ) x \,{\mathrm e}^{m \ln \left (e x +d \right )}}{e^{2} \left (m^{3}+6 m^{2}+11 m +6\right )}\) \(253\)
risch \(\frac {\left (B b \,e^{3} m^{2} x^{3}+A b \,e^{3} m^{2} x^{2}+B a \,e^{3} m^{2} x^{2}+B b d \,e^{2} m^{2} x^{2}+3 B b \,e^{3} m \,x^{3}+A a \,e^{3} m^{2} x +A b d \,e^{2} m^{2} x +4 A b \,e^{3} m \,x^{2}+B a d \,e^{2} m^{2} x +4 B a \,e^{3} m \,x^{2}+B b d \,e^{2} m \,x^{2}+2 b B \,x^{3} e^{3}+A a d \,e^{2} m^{2}+5 A a \,e^{3} m x +3 A b d \,e^{2} m x +3 A b \,e^{3} x^{2}+3 B a d \,e^{2} m x +3 B a \,e^{3} x^{2}-2 B b \,d^{2} e m x +5 A a d \,e^{2} m +6 A a \,e^{3} x -A b \,d^{2} e m -B a \,d^{2} e m +6 A a d \,e^{2}-3 A b \,d^{2} e -3 B a \,d^{2} e +2 b B \,d^{3}\right ) \left (e x +d \right )^{m}}{\left (2+m \right ) \left (3+m \right ) \left (1+m \right ) e^{3}}\) \(298\)
parallelrisch \(\frac {2 B \left (e x +d \right )^{m} b \,d^{4}+B \,x^{3} \left (e x +d \right )^{m} b d \,e^{3} m^{2}+A \,x^{2} \left (e x +d \right )^{m} b d \,e^{3} m^{2}+3 B \,x^{3} \left (e x +d \right )^{m} b d \,e^{3} m +B \,x^{2} \left (e x +d \right )^{m} a d \,e^{3} m^{2}+B \,x^{2} \left (e x +d \right )^{m} b \,d^{2} e^{2} m^{2}+4 A \,x^{2} \left (e x +d \right )^{m} b d \,e^{3} m +A x \left (e x +d \right )^{m} a d \,e^{3} m^{2}+A x \left (e x +d \right )^{m} b \,d^{2} e^{2} m^{2}+4 B \,x^{2} \left (e x +d \right )^{m} a d \,e^{3} m +B \,x^{2} \left (e x +d \right )^{m} b \,d^{2} e^{2} m +B x \left (e x +d \right )^{m} a \,d^{2} e^{2} m^{2}+5 A x \left (e x +d \right )^{m} a d \,e^{3} m +3 A x \left (e x +d \right )^{m} b \,d^{2} e^{2} m +3 B x \left (e x +d \right )^{m} a \,d^{2} e^{2} m -2 B x \left (e x +d \right )^{m} b \,d^{3} e m +2 B \,x^{3} \left (e x +d \right )^{m} b d \,e^{3}+3 A \,x^{2} \left (e x +d \right )^{m} b d \,e^{3}+A \left (e x +d \right )^{m} a \,d^{2} e^{2} m^{2}+3 B \,x^{2} \left (e x +d \right )^{m} a d \,e^{3}+6 A \left (e x +d \right )^{m} a \,d^{2} e^{2}-3 A \left (e x +d \right )^{m} b \,d^{3} e -3 B \left (e x +d \right )^{m} a \,d^{3} e +6 A x \left (e x +d \right )^{m} a d \,e^{3}+5 A \left (e x +d \right )^{m} a \,d^{2} e^{2} m -A \left (e x +d \right )^{m} b \,d^{3} e m -B \left (e x +d \right )^{m} a \,d^{3} e m}{\left (3+m \right ) \left (2+m \right ) \left (1+m \right ) e^{3} d}\) \(513\)

Input:

int((b*x+a)*(B*x+A)*(e*x+d)^m,x,method=_RETURNVERBOSE)
 

Output:

1/e^3*(e*x+d)^(1+m)/(m^3+6*m^2+11*m+6)*(B*b*e^2*m^2*x^2+A*b*e^2*m^2*x+B*a* 
e^2*m^2*x+3*B*b*e^2*m*x^2+A*a*e^2*m^2+4*A*b*e^2*m*x+4*B*a*e^2*m*x-2*B*b*d* 
e*m*x+2*B*b*e^2*x^2+5*A*a*e^2*m-A*b*d*e*m+3*A*b*e^2*x-B*a*d*e*m+3*B*a*e^2* 
x-2*B*b*d*e*x+6*A*a*e^2-3*A*b*d*e-3*B*a*d*e+2*B*b*d^2)
 

Fricas [B] (verification not implemented)

Leaf count of result is larger than twice the leaf count of optimal. 256 vs. \(2 (90) = 180\).

Time = 0.09 (sec) , antiderivative size = 256, normalized size of antiderivative = 2.84 \[ \int (a+b x) (A+B x) (d+e x)^m \, dx=\frac {{\left (A a d e^{2} m^{2} + 2 \, B b d^{3} + 6 \, A a d e^{2} - 3 \, {\left (B a + A b\right )} d^{2} e + {\left (B b e^{3} m^{2} + 3 \, B b e^{3} m + 2 \, B b e^{3}\right )} x^{3} + {\left (3 \, {\left (B a + A b\right )} e^{3} + {\left (B b d e^{2} + {\left (B a + A b\right )} e^{3}\right )} m^{2} + {\left (B b d e^{2} + 4 \, {\left (B a + A b\right )} e^{3}\right )} m\right )} x^{2} + {\left (5 \, A a d e^{2} - {\left (B a + A b\right )} d^{2} e\right )} m + {\left (6 \, A a e^{3} + {\left (A a e^{3} + {\left (B a + A b\right )} d e^{2}\right )} m^{2} - {\left (2 \, B b d^{2} e - 5 \, A a e^{3} - 3 \, {\left (B a + A b\right )} d e^{2}\right )} m\right )} x\right )} {\left (e x + d\right )}^{m}}{e^{3} m^{3} + 6 \, e^{3} m^{2} + 11 \, e^{3} m + 6 \, e^{3}} \] Input:

integrate((b*x+a)*(B*x+A)*(e*x+d)^m,x, algorithm="fricas")
 

Output:

(A*a*d*e^2*m^2 + 2*B*b*d^3 + 6*A*a*d*e^2 - 3*(B*a + A*b)*d^2*e + (B*b*e^3* 
m^2 + 3*B*b*e^3*m + 2*B*b*e^3)*x^3 + (3*(B*a + A*b)*e^3 + (B*b*d*e^2 + (B* 
a + A*b)*e^3)*m^2 + (B*b*d*e^2 + 4*(B*a + A*b)*e^3)*m)*x^2 + (5*A*a*d*e^2 
- (B*a + A*b)*d^2*e)*m + (6*A*a*e^3 + (A*a*e^3 + (B*a + A*b)*d*e^2)*m^2 - 
(2*B*b*d^2*e - 5*A*a*e^3 - 3*(B*a + A*b)*d*e^2)*m)*x)*(e*x + d)^m/(e^3*m^3 
 + 6*e^3*m^2 + 11*e^3*m + 6*e^3)
 

Sympy [B] (verification not implemented)

Leaf count of result is larger than twice the leaf count of optimal. 1982 vs. \(2 (78) = 156\).

Time = 0.65 (sec) , antiderivative size = 1982, normalized size of antiderivative = 22.02 \[ \int (a+b x) (A+B x) (d+e x)^m \, dx=\text {Too large to display} \] Input:

integrate((b*x+a)*(B*x+A)*(e*x+d)**m,x)
 

Output:

Piecewise((d**m*(A*a*x + A*b*x**2/2 + B*a*x**2/2 + B*b*x**3/3), Eq(e, 0)), 
 (-A*a*e**2/(2*d**2*e**3 + 4*d*e**4*x + 2*e**5*x**2) - A*b*d*e/(2*d**2*e** 
3 + 4*d*e**4*x + 2*e**5*x**2) - 2*A*b*e**2*x/(2*d**2*e**3 + 4*d*e**4*x + 2 
*e**5*x**2) - B*a*d*e/(2*d**2*e**3 + 4*d*e**4*x + 2*e**5*x**2) - 2*B*a*e** 
2*x/(2*d**2*e**3 + 4*d*e**4*x + 2*e**5*x**2) + 2*B*b*d**2*log(d/e + x)/(2* 
d**2*e**3 + 4*d*e**4*x + 2*e**5*x**2) + 3*B*b*d**2/(2*d**2*e**3 + 4*d*e**4 
*x + 2*e**5*x**2) + 4*B*b*d*e*x*log(d/e + x)/(2*d**2*e**3 + 4*d*e**4*x + 2 
*e**5*x**2) + 4*B*b*d*e*x/(2*d**2*e**3 + 4*d*e**4*x + 2*e**5*x**2) + 2*B*b 
*e**2*x**2*log(d/e + x)/(2*d**2*e**3 + 4*d*e**4*x + 2*e**5*x**2), Eq(m, -3 
)), (-A*a*e**2/(d*e**3 + e**4*x) + A*b*d*e*log(d/e + x)/(d*e**3 + e**4*x) 
+ A*b*d*e/(d*e**3 + e**4*x) + A*b*e**2*x*log(d/e + x)/(d*e**3 + e**4*x) + 
B*a*d*e*log(d/e + x)/(d*e**3 + e**4*x) + B*a*d*e/(d*e**3 + e**4*x) + B*a*e 
**2*x*log(d/e + x)/(d*e**3 + e**4*x) - 2*B*b*d**2*log(d/e + x)/(d*e**3 + e 
**4*x) - 2*B*b*d**2/(d*e**3 + e**4*x) - 2*B*b*d*e*x*log(d/e + x)/(d*e**3 + 
 e**4*x) + B*b*e**2*x**2/(d*e**3 + e**4*x), Eq(m, -2)), (A*a*log(d/e + x)/ 
e - A*b*d*log(d/e + x)/e**2 + A*b*x/e - B*a*d*log(d/e + x)/e**2 + B*a*x/e 
+ B*b*d**2*log(d/e + x)/e**3 - B*b*d*x/e**2 + B*b*x**2/(2*e), Eq(m, -1)), 
(A*a*d*e**2*m**2*(d + e*x)**m/(e**3*m**3 + 6*e**3*m**2 + 11*e**3*m + 6*e** 
3) + 5*A*a*d*e**2*m*(d + e*x)**m/(e**3*m**3 + 6*e**3*m**2 + 11*e**3*m + 6* 
e**3) + 6*A*a*d*e**2*(d + e*x)**m/(e**3*m**3 + 6*e**3*m**2 + 11*e**3*m ...
 

Maxima [A] (verification not implemented)

Time = 0.04 (sec) , antiderivative size = 179, normalized size of antiderivative = 1.99 \[ \int (a+b x) (A+B x) (d+e x)^m \, dx=\frac {{\left (e^{2} {\left (m + 1\right )} x^{2} + d e m x - d^{2}\right )} {\left (e x + d\right )}^{m} B a}{{\left (m^{2} + 3 \, m + 2\right )} e^{2}} + \frac {{\left (e^{2} {\left (m + 1\right )} x^{2} + d e m x - d^{2}\right )} {\left (e x + d\right )}^{m} A b}{{\left (m^{2} + 3 \, m + 2\right )} e^{2}} + \frac {{\left (e x + d\right )}^{m + 1} A a}{e {\left (m + 1\right )}} + \frac {{\left ({\left (m^{2} + 3 \, m + 2\right )} e^{3} x^{3} + {\left (m^{2} + m\right )} d e^{2} x^{2} - 2 \, d^{2} e m x + 2 \, d^{3}\right )} {\left (e x + d\right )}^{m} B b}{{\left (m^{3} + 6 \, m^{2} + 11 \, m + 6\right )} e^{3}} \] Input:

integrate((b*x+a)*(B*x+A)*(e*x+d)^m,x, algorithm="maxima")
 

Output:

(e^2*(m + 1)*x^2 + d*e*m*x - d^2)*(e*x + d)^m*B*a/((m^2 + 3*m + 2)*e^2) + 
(e^2*(m + 1)*x^2 + d*e*m*x - d^2)*(e*x + d)^m*A*b/((m^2 + 3*m + 2)*e^2) + 
(e*x + d)^(m + 1)*A*a/(e*(m + 1)) + ((m^2 + 3*m + 2)*e^3*x^3 + (m^2 + m)*d 
*e^2*x^2 - 2*d^2*e*m*x + 2*d^3)*(e*x + d)^m*B*b/((m^3 + 6*m^2 + 11*m + 6)* 
e^3)
 

Giac [B] (verification not implemented)

Leaf count of result is larger than twice the leaf count of optimal. 490 vs. \(2 (90) = 180\).

Time = 0.13 (sec) , antiderivative size = 490, normalized size of antiderivative = 5.44 \[ \int (a+b x) (A+B x) (d+e x)^m \, dx=\frac {{\left (e x + d\right )}^{m} B b e^{3} m^{2} x^{3} + {\left (e x + d\right )}^{m} B b d e^{2} m^{2} x^{2} + {\left (e x + d\right )}^{m} B a e^{3} m^{2} x^{2} + {\left (e x + d\right )}^{m} A b e^{3} m^{2} x^{2} + 3 \, {\left (e x + d\right )}^{m} B b e^{3} m x^{3} + {\left (e x + d\right )}^{m} B a d e^{2} m^{2} x + {\left (e x + d\right )}^{m} A b d e^{2} m^{2} x + {\left (e x + d\right )}^{m} A a e^{3} m^{2} x + {\left (e x + d\right )}^{m} B b d e^{2} m x^{2} + 4 \, {\left (e x + d\right )}^{m} B a e^{3} m x^{2} + 4 \, {\left (e x + d\right )}^{m} A b e^{3} m x^{2} + 2 \, {\left (e x + d\right )}^{m} B b e^{3} x^{3} + {\left (e x + d\right )}^{m} A a d e^{2} m^{2} - 2 \, {\left (e x + d\right )}^{m} B b d^{2} e m x + 3 \, {\left (e x + d\right )}^{m} B a d e^{2} m x + 3 \, {\left (e x + d\right )}^{m} A b d e^{2} m x + 5 \, {\left (e x + d\right )}^{m} A a e^{3} m x + 3 \, {\left (e x + d\right )}^{m} B a e^{3} x^{2} + 3 \, {\left (e x + d\right )}^{m} A b e^{3} x^{2} - {\left (e x + d\right )}^{m} B a d^{2} e m - {\left (e x + d\right )}^{m} A b d^{2} e m + 5 \, {\left (e x + d\right )}^{m} A a d e^{2} m + 6 \, {\left (e x + d\right )}^{m} A a e^{3} x + 2 \, {\left (e x + d\right )}^{m} B b d^{3} - 3 \, {\left (e x + d\right )}^{m} B a d^{2} e - 3 \, {\left (e x + d\right )}^{m} A b d^{2} e + 6 \, {\left (e x + d\right )}^{m} A a d e^{2}}{e^{3} m^{3} + 6 \, e^{3} m^{2} + 11 \, e^{3} m + 6 \, e^{3}} \] Input:

integrate((b*x+a)*(B*x+A)*(e*x+d)^m,x, algorithm="giac")
 

Output:

((e*x + d)^m*B*b*e^3*m^2*x^3 + (e*x + d)^m*B*b*d*e^2*m^2*x^2 + (e*x + d)^m 
*B*a*e^3*m^2*x^2 + (e*x + d)^m*A*b*e^3*m^2*x^2 + 3*(e*x + d)^m*B*b*e^3*m*x 
^3 + (e*x + d)^m*B*a*d*e^2*m^2*x + (e*x + d)^m*A*b*d*e^2*m^2*x + (e*x + d) 
^m*A*a*e^3*m^2*x + (e*x + d)^m*B*b*d*e^2*m*x^2 + 4*(e*x + d)^m*B*a*e^3*m*x 
^2 + 4*(e*x + d)^m*A*b*e^3*m*x^2 + 2*(e*x + d)^m*B*b*e^3*x^3 + (e*x + d)^m 
*A*a*d*e^2*m^2 - 2*(e*x + d)^m*B*b*d^2*e*m*x + 3*(e*x + d)^m*B*a*d*e^2*m*x 
 + 3*(e*x + d)^m*A*b*d*e^2*m*x + 5*(e*x + d)^m*A*a*e^3*m*x + 3*(e*x + d)^m 
*B*a*e^3*x^2 + 3*(e*x + d)^m*A*b*e^3*x^2 - (e*x + d)^m*B*a*d^2*e*m - (e*x 
+ d)^m*A*b*d^2*e*m + 5*(e*x + d)^m*A*a*d*e^2*m + 6*(e*x + d)^m*A*a*e^3*x + 
 2*(e*x + d)^m*B*b*d^3 - 3*(e*x + d)^m*B*a*d^2*e - 3*(e*x + d)^m*A*b*d^2*e 
 + 6*(e*x + d)^m*A*a*d*e^2)/(e^3*m^3 + 6*e^3*m^2 + 11*e^3*m + 6*e^3)
 

Mupad [B] (verification not implemented)

Time = 1.21 (sec) , antiderivative size = 259, normalized size of antiderivative = 2.88 \[ \int (a+b x) (A+B x) (d+e x)^m \, dx={\left (d+e\,x\right )}^m\,\left (\frac {x\,\left (6\,A\,a\,e^3+5\,A\,a\,e^3\,m+A\,a\,e^3\,m^2+3\,A\,b\,d\,e^2\,m+3\,B\,a\,d\,e^2\,m-2\,B\,b\,d^2\,e\,m+A\,b\,d\,e^2\,m^2+B\,a\,d\,e^2\,m^2\right )}{e^3\,\left (m^3+6\,m^2+11\,m+6\right )}+\frac {d\,\left (6\,A\,a\,e^2+2\,B\,b\,d^2+5\,A\,a\,e^2\,m+A\,a\,e^2\,m^2-3\,A\,b\,d\,e-3\,B\,a\,d\,e-A\,b\,d\,e\,m-B\,a\,d\,e\,m\right )}{e^3\,\left (m^3+6\,m^2+11\,m+6\right )}+\frac {B\,b\,x^3\,\left (m^2+3\,m+2\right )}{m^3+6\,m^2+11\,m+6}+\frac {x^2\,\left (m+1\right )\,\left (3\,A\,b\,e+3\,B\,a\,e+A\,b\,e\,m+B\,a\,e\,m+B\,b\,d\,m\right )}{e\,\left (m^3+6\,m^2+11\,m+6\right )}\right ) \] Input:

int((A + B*x)*(a + b*x)*(d + e*x)^m,x)
 

Output:

(d + e*x)^m*((x*(6*A*a*e^3 + 5*A*a*e^3*m + A*a*e^3*m^2 + 3*A*b*d*e^2*m + 3 
*B*a*d*e^2*m - 2*B*b*d^2*e*m + A*b*d*e^2*m^2 + B*a*d*e^2*m^2))/(e^3*(11*m 
+ 6*m^2 + m^3 + 6)) + (d*(6*A*a*e^2 + 2*B*b*d^2 + 5*A*a*e^2*m + A*a*e^2*m^ 
2 - 3*A*b*d*e - 3*B*a*d*e - A*b*d*e*m - B*a*d*e*m))/(e^3*(11*m + 6*m^2 + m 
^3 + 6)) + (B*b*x^3*(3*m + m^2 + 2))/(11*m + 6*m^2 + m^3 + 6) + (x^2*(m + 
1)*(3*A*b*e + 3*B*a*e + A*b*e*m + B*a*e*m + B*b*d*m))/(e*(11*m + 6*m^2 + m 
^3 + 6)))
 

Reduce [B] (verification not implemented)

Time = 0.16 (sec) , antiderivative size = 241, normalized size of antiderivative = 2.68 \[ \int (a+b x) (A+B x) (d+e x)^m \, dx=\frac {\left (e x +d \right )^{m} \left (b^{2} e^{3} m^{2} x^{3}+2 a b \,e^{3} m^{2} x^{2}+b^{2} d \,e^{2} m^{2} x^{2}+3 b^{2} e^{3} m \,x^{3}+a^{2} e^{3} m^{2} x +2 a b d \,e^{2} m^{2} x +8 a b \,e^{3} m \,x^{2}+b^{2} d \,e^{2} m \,x^{2}+2 b^{2} e^{3} x^{3}+a^{2} d \,e^{2} m^{2}+5 a^{2} e^{3} m x +6 a b d \,e^{2} m x +6 a b \,e^{3} x^{2}-2 b^{2} d^{2} e m x +5 a^{2} d \,e^{2} m +6 a^{2} e^{3} x -2 a b \,d^{2} e m +6 a^{2} d \,e^{2}-6 a b \,d^{2} e +2 b^{2} d^{3}\right )}{e^{3} \left (m^{3}+6 m^{2}+11 m +6\right )} \] Input:

int((b*x+a)*(B*x+A)*(e*x+d)^m,x)
 

Output:

((d + e*x)**m*(a**2*d*e**2*m**2 + 5*a**2*d*e**2*m + 6*a**2*d*e**2 + a**2*e 
**3*m**2*x + 5*a**2*e**3*m*x + 6*a**2*e**3*x - 2*a*b*d**2*e*m - 6*a*b*d**2 
*e + 2*a*b*d*e**2*m**2*x + 6*a*b*d*e**2*m*x + 2*a*b*e**3*m**2*x**2 + 8*a*b 
*e**3*m*x**2 + 6*a*b*e**3*x**2 + 2*b**2*d**3 - 2*b**2*d**2*e*m*x + b**2*d* 
e**2*m**2*x**2 + b**2*d*e**2*m*x**2 + b**2*e**3*m**2*x**3 + 3*b**2*e**3*m* 
x**3 + 2*b**2*e**3*x**3))/(e**3*(m**3 + 6*m**2 + 11*m + 6))