\(\int \frac {(a+b x+c x^2)^m (d+e x+f x^2+g x^3)^n (-2 a d+(-b d-a e+b d m+a e n) x+(2 c d m+b e m+b e n+2 a f n) x^2+(c e+b f+a g+2 c e m+b f m+c e n+2 b f n+3 a g n) x^3+(2 c f+2 b g+2 c f m+b g m+2 c f n+3 b g n) x^4+c g (3+2 m+3 n) x^5)}{x^3} \, dx\) [617]

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

Optimal result

Integrand size = 163, antiderivative size = 37 \[ \int \frac {\left (a+b x+c x^2\right )^m \left (d+e x+f x^2+g x^3\right )^n \left (-2 a d+(-b d-a e+b d m+a e n) x+(2 c d m+b e m+b e n+2 a f n) x^2+(c e+b f+a g+2 c e m+b f m+c e n+2 b f n+3 a g n) x^3+(2 c f+2 b g+2 c f m+b g m+2 c f n+3 b g n) x^4+c g (3+2 m+3 n) x^5\right )}{x^3} \, dx=\frac {\left (a+b x+c x^2\right )^{1+m} \left (d+e x+f x^2+g x^3\right )^{1+n}}{x^2} \]

[Out]

(c*x^2+b*x+a)^(1+m)*(g*x^3+f*x^2+e*x+d)^(1+n)/x^2

Rubi [F]

\[ \int \frac {\left (a+b x+c x^2\right )^m \left (d+e x+f x^2+g x^3\right )^n \left (-2 a d+(-b d-a e+b d m+a e n) x+(2 c d m+b e m+b e n+2 a f n) x^2+(c e+b f+a g+2 c e m+b f m+c e n+2 b f n+3 a g n) x^3+(2 c f+2 b g+2 c f m+b g m+2 c f n+3 b g n) x^4+c g (3+2 m+3 n) x^5\right )}{x^3} \, dx=\int \frac {\left (a+b x+c x^2\right )^m \left (d+e x+f x^2+g x^3\right )^n \left (-2 a d+(-b d-a e+b d m+a e n) x+(2 c d m+b e m+b e n+2 a f n) x^2+(c e+b f+a g+2 c e m+b f m+c e n+2 b f n+3 a g n) x^3+(2 c f+2 b g+2 c f m+b g m+2 c f n+3 b g n) x^4+c g (3+2 m+3 n) x^5\right )}{x^3} \, dx \]

[In]

Int[((a + b*x + c*x^2)^m*(d + e*x + f*x^2 + g*x^3)^n*(-2*a*d + (-(b*d) - a*e + b*d*m + a*e*n)*x + (2*c*d*m + b
*e*m + b*e*n + 2*a*f*n)*x^2 + (c*e + b*f + a*g + 2*c*e*m + b*f*m + c*e*n + 2*b*f*n + 3*a*g*n)*x^3 + (2*c*f + 2
*b*g + 2*c*f*m + b*g*m + 2*c*f*n + 3*b*g*n)*x^4 + c*g*(3 + 2*m + 3*n)*x^5))/x^3,x]

[Out]

(c*e*(1 + 2*m + n) + b*f*(1 + m + 2*n) + a*g*(1 + 3*n))*Defer[Int][(a + b*x + c*x^2)^m*(d + e*x + f*x^2 + g*x^
3)^n, x] - 2*a*d*Defer[Int][((a + b*x + c*x^2)^m*(d + e*x + f*x^2 + g*x^3)^n)/x^3, x] - (b*d*(1 - m) + a*e*(1
- n))*Defer[Int][((a + b*x + c*x^2)^m*(d + e*x + f*x^2 + g*x^3)^n)/x^2, x] + (2*c*d*m + 2*a*f*n + b*e*(m + n))
*Defer[Int][((a + b*x + c*x^2)^m*(d + e*x + f*x^2 + g*x^3)^n)/x, x] + (2*c*f*(1 + m + n) + b*g*(2 + m + 3*n))*
Defer[Int][x*(a + b*x + c*x^2)^m*(d + e*x + f*x^2 + g*x^3)^n, x] + c*g*(3 + 2*m + 3*n)*Defer[Int][x^2*(a + b*x
 + c*x^2)^m*(d + e*x + f*x^2 + g*x^3)^n, x]

Rubi steps \begin{align*} \text {integral}& = \int \left (c e \left (1+\frac {c e (2 m+n)+b f (1+m+2 n)+a (g+3 g n)}{c e}\right ) \left (a+b x+c x^2\right )^m \left (d+e x+f x^2+g x^3\right )^n-\frac {2 a d \left (a+b x+c x^2\right )^m \left (d+e x+f x^2+g x^3\right )^n}{x^3}+\frac {(-b d (1-m)-a e (1-n)) \left (a+b x+c x^2\right )^m \left (d+e x+f x^2+g x^3\right )^n}{x^2}+\frac {(2 c d m+2 a f n+b e (m+n)) \left (a+b x+c x^2\right )^m \left (d+e x+f x^2+g x^3\right )^n}{x}+(2 c f (1+m+n)+b g (2+m+3 n)) x \left (a+b x+c x^2\right )^m \left (d+e x+f x^2+g x^3\right )^n+c g (3+2 m+3 n) x^2 \left (a+b x+c x^2\right )^m \left (d+e x+f x^2+g x^3\right )^n\right ) \, dx \\ & = -\left ((2 a d) \int \frac {\left (a+b x+c x^2\right )^m \left (d+e x+f x^2+g x^3\right )^n}{x^3} \, dx\right )+(-b d (1-m)-a e (1-n)) \int \frac {\left (a+b x+c x^2\right )^m \left (d+e x+f x^2+g x^3\right )^n}{x^2} \, dx+(c g (3+2 m+3 n)) \int x^2 \left (a+b x+c x^2\right )^m \left (d+e x+f x^2+g x^3\right )^n \, dx+(2 c d m+2 a f n+b e (m+n)) \int \frac {\left (a+b x+c x^2\right )^m \left (d+e x+f x^2+g x^3\right )^n}{x} \, dx+(c e (1+2 m+n)+b f (1+m+2 n)+a g (1+3 n)) \int \left (a+b x+c x^2\right )^m \left (d+e x+f x^2+g x^3\right )^n \, dx+(2 c f (1+m+n)+b g (2+m+3 n)) \int x \left (a+b x+c x^2\right )^m \left (d+e x+f x^2+g x^3\right )^n \, dx \\ \end{align*}

Mathematica [A] (verified)

Time = 9.15 (sec) , antiderivative size = 34, normalized size of antiderivative = 0.92 \[ \int \frac {\left (a+b x+c x^2\right )^m \left (d+e x+f x^2+g x^3\right )^n \left (-2 a d+(-b d-a e+b d m+a e n) x+(2 c d m+b e m+b e n+2 a f n) x^2+(c e+b f+a g+2 c e m+b f m+c e n+2 b f n+3 a g n) x^3+(2 c f+2 b g+2 c f m+b g m+2 c f n+3 b g n) x^4+c g (3+2 m+3 n) x^5\right )}{x^3} \, dx=\frac {(a+x (b+c x))^{1+m} (d+x (e+x (f+g x)))^{1+n}}{x^2} \]

[In]

Integrate[((a + b*x + c*x^2)^m*(d + e*x + f*x^2 + g*x^3)^n*(-2*a*d + (-(b*d) - a*e + b*d*m + a*e*n)*x + (2*c*d
*m + b*e*m + b*e*n + 2*a*f*n)*x^2 + (c*e + b*f + a*g + 2*c*e*m + b*f*m + c*e*n + 2*b*f*n + 3*a*g*n)*x^3 + (2*c
*f + 2*b*g + 2*c*f*m + b*g*m + 2*c*f*n + 3*b*g*n)*x^4 + c*g*(3 + 2*m + 3*n)*x^5))/x^3,x]

[Out]

((a + x*(b + c*x))^(1 + m)*(d + x*(e + x*(f + g*x)))^(1 + n))/x^2

Maple [A] (verified)

Time = 82.80 (sec) , antiderivative size = 38, normalized size of antiderivative = 1.03

method result size
gosper \(\frac {\left (c \,x^{2}+b x +a \right )^{1+m} \left (g \,x^{3}+f \,x^{2}+e x +d \right )^{1+n}}{x^{2}}\) \(38\)
risch \(\frac {\left (c g \,x^{5}+b g \,x^{4}+c f \,x^{4}+a g \,x^{3}+b f \,x^{3}+c \,x^{3} e +a f \,x^{2}+e \,x^{2} b +x^{2} c d +a e x +b d x +a d \right ) \left (c \,x^{2}+b x +a \right )^{m} \left (g \,x^{3}+f \,x^{2}+e x +d \right )^{n}}{x^{2}}\) \(100\)
parallelrisch \(\frac {x^{5} \left (c \,x^{2}+b x +a \right )^{m} \left (g \,x^{3}+f \,x^{2}+e x +d \right )^{n} c^{2} g^{2}+x^{4} \left (c \,x^{2}+b x +a \right )^{m} \left (g \,x^{3}+f \,x^{2}+e x +d \right )^{n} b c \,g^{2}+x^{4} \left (c \,x^{2}+b x +a \right )^{m} \left (g \,x^{3}+f \,x^{2}+e x +d \right )^{n} c^{2} f g +x^{3} \left (c \,x^{2}+b x +a \right )^{m} \left (g \,x^{3}+f \,x^{2}+e x +d \right )^{n} a c \,g^{2}+x^{3} \left (c \,x^{2}+b x +a \right )^{m} \left (g \,x^{3}+f \,x^{2}+e x +d \right )^{n} b c f g +x^{3} \left (c \,x^{2}+b x +a \right )^{m} \left (g \,x^{3}+f \,x^{2}+e x +d \right )^{n} c^{2} e g +x^{2} \left (c \,x^{2}+b x +a \right )^{m} \left (g \,x^{3}+f \,x^{2}+e x +d \right )^{n} a c f g +x^{2} \left (c \,x^{2}+b x +a \right )^{m} \left (g \,x^{3}+f \,x^{2}+e x +d \right )^{n} b c e g +x^{2} \left (c \,x^{2}+b x +a \right )^{m} \left (g \,x^{3}+f \,x^{2}+e x +d \right )^{n} c^{2} d g +x \left (c \,x^{2}+b x +a \right )^{m} \left (g \,x^{3}+f \,x^{2}+e x +d \right )^{n} a c e g +x \left (c \,x^{2}+b x +a \right )^{m} \left (g \,x^{3}+f \,x^{2}+e x +d \right )^{n} b c d g +\left (c \,x^{2}+b x +a \right )^{m} \left (g \,x^{3}+f \,x^{2}+e x +d \right )^{n} a c d g}{x^{2} g c}\) \(456\)

[In]

int((c*x^2+b*x+a)^m*(g*x^3+f*x^2+e*x+d)^n*(-2*a*d+(a*e*n+b*d*m-a*e-b*d)*x+(2*a*f*n+b*e*m+b*e*n+2*c*d*m)*x^2+(3
*a*g*n+b*f*m+2*b*f*n+2*c*e*m+c*e*n+a*g+b*f+c*e)*x^3+(b*g*m+3*b*g*n+2*c*f*m+2*c*f*n+2*b*g+2*c*f)*x^4+c*g*(3+2*m
+3*n)*x^5)/x^3,x,method=_RETURNVERBOSE)

[Out]

(c*x^2+b*x+a)^(1+m)*(g*x^3+f*x^2+e*x+d)^(1+n)/x^2

Fricas [F(-1)]

Timed out. \[ \int \frac {\left (a+b x+c x^2\right )^m \left (d+e x+f x^2+g x^3\right )^n \left (-2 a d+(-b d-a e+b d m+a e n) x+(2 c d m+b e m+b e n+2 a f n) x^2+(c e+b f+a g+2 c e m+b f m+c e n+2 b f n+3 a g n) x^3+(2 c f+2 b g+2 c f m+b g m+2 c f n+3 b g n) x^4+c g (3+2 m+3 n) x^5\right )}{x^3} \, dx=\text {Timed out} \]

[In]

integrate((c*x^2+b*x+a)^m*(g*x^3+f*x^2+e*x+d)^n*(-2*a*d+(a*e*n+b*d*m-a*e-b*d)*x+(2*a*f*n+b*e*m+b*e*n+2*c*d*m)*
x^2+(3*a*g*n+b*f*m+2*b*f*n+2*c*e*m+c*e*n+a*g+b*f+c*e)*x^3+(b*g*m+3*b*g*n+2*c*f*m+2*c*f*n+2*b*g+2*c*f)*x^4+c*g*
(3+2*m+3*n)*x^5)/x^3,x, algorithm="fricas")

[Out]

Timed out

Sympy [F(-1)]

Timed out. \[ \int \frac {\left (a+b x+c x^2\right )^m \left (d+e x+f x^2+g x^3\right )^n \left (-2 a d+(-b d-a e+b d m+a e n) x+(2 c d m+b e m+b e n+2 a f n) x^2+(c e+b f+a g+2 c e m+b f m+c e n+2 b f n+3 a g n) x^3+(2 c f+2 b g+2 c f m+b g m+2 c f n+3 b g n) x^4+c g (3+2 m+3 n) x^5\right )}{x^3} \, dx=\text {Timed out} \]

[In]

integrate((c*x**2+b*x+a)**m*(g*x**3+f*x**2+e*x+d)**n*(-2*a*d+(a*e*n+b*d*m-a*e-b*d)*x+(2*a*f*n+b*e*m+b*e*n+2*c*
d*m)*x**2+(3*a*g*n+b*f*m+2*b*f*n+2*c*e*m+c*e*n+a*g+b*f+c*e)*x**3+(b*g*m+3*b*g*n+2*c*f*m+2*c*f*n+2*b*g+2*c*f)*x
**4+c*g*(3+2*m+3*n)*x**5)/x**3,x)

[Out]

Timed out

Maxima [B] (verification not implemented)

Leaf count of result is larger than twice the leaf count of optimal. 95 vs. \(2 (37) = 74\).

Time = 0.28 (sec) , antiderivative size = 95, normalized size of antiderivative = 2.57 \[ \int \frac {\left (a+b x+c x^2\right )^m \left (d+e x+f x^2+g x^3\right )^n \left (-2 a d+(-b d-a e+b d m+a e n) x+(2 c d m+b e m+b e n+2 a f n) x^2+(c e+b f+a g+2 c e m+b f m+c e n+2 b f n+3 a g n) x^3+(2 c f+2 b g+2 c f m+b g m+2 c f n+3 b g n) x^4+c g (3+2 m+3 n) x^5\right )}{x^3} \, dx=\frac {{\left (c g x^{5} + {\left (c f + b g\right )} x^{4} + {\left (c e + b f + a g\right )} x^{3} + {\left (c d + b e + a f\right )} x^{2} + a d + {\left (b d + a e\right )} x\right )} e^{\left (n \log \left (g x^{3} + f x^{2} + e x + d\right ) + m \log \left (c x^{2} + b x + a\right )\right )}}{x^{2}} \]

[In]

integrate((c*x^2+b*x+a)^m*(g*x^3+f*x^2+e*x+d)^n*(-2*a*d+(a*e*n+b*d*m-a*e-b*d)*x+(2*a*f*n+b*e*m+b*e*n+2*c*d*m)*
x^2+(3*a*g*n+b*f*m+2*b*f*n+2*c*e*m+c*e*n+a*g+b*f+c*e)*x^3+(b*g*m+3*b*g*n+2*c*f*m+2*c*f*n+2*b*g+2*c*f)*x^4+c*g*
(3+2*m+3*n)*x^5)/x^3,x, algorithm="maxima")

[Out]

(c*g*x^5 + (c*f + b*g)*x^4 + (c*e + b*f + a*g)*x^3 + (c*d + b*e + a*f)*x^2 + a*d + (b*d + a*e)*x)*e^(n*log(g*x
^3 + f*x^2 + e*x + d) + m*log(c*x^2 + b*x + a))/x^2

Giac [F(-1)]

Timed out. \[ \int \frac {\left (a+b x+c x^2\right )^m \left (d+e x+f x^2+g x^3\right )^n \left (-2 a d+(-b d-a e+b d m+a e n) x+(2 c d m+b e m+b e n+2 a f n) x^2+(c e+b f+a g+2 c e m+b f m+c e n+2 b f n+3 a g n) x^3+(2 c f+2 b g+2 c f m+b g m+2 c f n+3 b g n) x^4+c g (3+2 m+3 n) x^5\right )}{x^3} \, dx=\text {Timed out} \]

[In]

integrate((c*x^2+b*x+a)^m*(g*x^3+f*x^2+e*x+d)^n*(-2*a*d+(a*e*n+b*d*m-a*e-b*d)*x+(2*a*f*n+b*e*m+b*e*n+2*c*d*m)*
x^2+(3*a*g*n+b*f*m+2*b*f*n+2*c*e*m+c*e*n+a*g+b*f+c*e)*x^3+(b*g*m+3*b*g*n+2*c*f*m+2*c*f*n+2*b*g+2*c*f)*x^4+c*g*
(3+2*m+3*n)*x^5)/x^3,x, algorithm="giac")

[Out]

Timed out

Mupad [B] (verification not implemented)

Time = 26.02 (sec) , antiderivative size = 146, normalized size of antiderivative = 3.95 \[ \int \frac {\left (a+b x+c x^2\right )^m \left (d+e x+f x^2+g x^3\right )^n \left (-2 a d+(-b d-a e+b d m+a e n) x+(2 c d m+b e m+b e n+2 a f n) x^2+(c e+b f+a g+2 c e m+b f m+c e n+2 b f n+3 a g n) x^3+(2 c f+2 b g+2 c f m+b g m+2 c f n+3 b g n) x^4+c g (3+2 m+3 n) x^5\right )}{x^3} \, dx={\left (c\,x^2+b\,x+a\right )}^m\,{\left (g\,x^3+f\,x^2+e\,x+d\right )}^n\,\left (a\,f+b\,e+c\,d+c\,g\,x^3+a\,g\,x+b\,f\,x+c\,e\,x+b\,g\,x^2+c\,f\,x^2\right )+\frac {\left (a\,e+b\,d\right )\,{\left (c\,x^2+b\,x+a\right )}^m\,{\left (g\,x^3+f\,x^2+e\,x+d\right )}^n}{x}+\frac {a\,d\,{\left (c\,x^2+b\,x+a\right )}^m\,{\left (g\,x^3+f\,x^2+e\,x+d\right )}^n}{x^2} \]

[In]

int(((a + b*x + c*x^2)^m*(d + e*x + f*x^2 + g*x^3)^n*(x^4*(2*b*g + 2*c*f + b*g*m + 2*c*f*m + 3*b*g*n + 2*c*f*n
) - 2*a*d - x*(a*e + b*d - b*d*m - a*e*n) + x^2*(b*e*m + 2*c*d*m + 2*a*f*n + b*e*n) + x^3*(a*g + b*f + c*e + b
*f*m + 2*c*e*m + 3*a*g*n + 2*b*f*n + c*e*n) + c*g*x^5*(2*m + 3*n + 3)))/x^3,x)

[Out]

(a + b*x + c*x^2)^m*(d + e*x + f*x^2 + g*x^3)^n*(a*f + b*e + c*d + c*g*x^3 + a*g*x + b*f*x + c*e*x + b*g*x^2 +
 c*f*x^2) + ((a*e + b*d)*(a + b*x + c*x^2)^m*(d + e*x + f*x^2 + g*x^3)^n)/x + (a*d*(a + b*x + c*x^2)^m*(d + e*
x + f*x^2 + g*x^3)^n)/x^2