\(\int \frac {(a+b x)^n (c+d x)^2}{x^2} \, dx\) [546]

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

Optimal result

Integrand size = 18, antiderivative size = 87 \[ \int \frac {(a+b x)^n (c+d x)^2}{x^2} \, dx=\frac {d^2 (a+b x)^{1+n}}{b (1+n)}-\frac {c^2 (a+b x)^{1+n}}{a x}-\frac {c (2 a d+b c n) (a+b x)^{1+n} \operatorname {Hypergeometric2F1}\left (1,1+n,2+n,1+\frac {b x}{a}\right )}{a^2 (1+n)} \] Output:

d^2*(b*x+a)^(1+n)/b/(1+n)-c^2*(b*x+a)^(1+n)/a/x-c*(b*c*n+2*a*d)*(b*x+a)^(1 
+n)*hypergeom([1, 1+n],[2+n],1+b*x/a)/a^2/(1+n)
 

Mathematica [A] (verified)

Time = 0.06 (sec) , antiderivative size = 73, normalized size of antiderivative = 0.84 \[ \int \frac {(a+b x)^n (c+d x)^2}{x^2} \, dx=\frac {(a+b x)^{1+n} \left (a \left (-b c^2 (1+n)+a d^2 x\right )-b c (2 a d+b c n) x \operatorname {Hypergeometric2F1}\left (1,1+n,2+n,1+\frac {b x}{a}\right )\right )}{a^2 b (1+n) x} \] Input:

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

Output:

((a + b*x)^(1 + n)*(a*(-(b*c^2*(1 + n)) + a*d^2*x) - b*c*(2*a*d + b*c*n)*x 
*Hypergeometric2F1[1, 1 + n, 2 + n, 1 + (b*x)/a]))/(a^2*b*(1 + n)*x)
 

Rubi [A] (verified)

Time = 0.20 (sec) , antiderivative size = 93, normalized size of antiderivative = 1.07, number of steps used = 3, number of rules used = 3, \(\frac {\text {number of rules}}{\text {integrand size}}\) = 0.167, Rules used = {100, 90, 75}

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 {(c+d x)^2 (a+b x)^n}{x^2} \, dx\)

\(\Big \downarrow \) 100

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

\(\Big \downarrow \) 90

\(\displaystyle \frac {c (2 a d+b c n) \int \frac {(a+b x)^n}{x}dx+\frac {a d^2 (a+b x)^{n+1}}{b (n+1)}}{a}-\frac {c^2 (a+b x)^{n+1}}{a x}\)

\(\Big \downarrow \) 75

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

Input:

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

Output:

-((c^2*(a + b*x)^(1 + n))/(a*x)) + ((a*d^2*(a + b*x)^(1 + n))/(b*(1 + n)) 
- (c*(2*a*d + b*c*n)*(a + b*x)^(1 + n)*Hypergeometric2F1[1, 1 + n, 2 + n, 
1 + (b*x)/a])/(a*(1 + n)))/a
 

Defintions of rubi rules used

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

rule 90
Int[((a_.) + (b_.)*(x_))*((c_.) + (d_.)*(x_))^(n_.)*((e_.) + (f_.)*(x_))^(p 
_.), x_] :> Simp[b*(c + d*x)^(n + 1)*((e + f*x)^(p + 1)/(d*f*(n + p + 2))), 
 x] + Simp[(a*d*f*(n + p + 2) - b*(d*e*(n + 1) + c*f*(p + 1)))/(d*f*(n + p 
+ 2))   Int[(c + d*x)^n*(e + f*x)^p, x], x] /; FreeQ[{a, b, c, d, e, f, n, 
p}, x] && NeQ[n + p + 2, 0]
 

rule 100
Int[((a_.) + (b_.)*(x_))^2*((c_.) + (d_.)*(x_))^(n_)*((e_.) + (f_.)*(x_))^( 
p_), x_] :> Simp[(b*c - a*d)^2*(c + d*x)^(n + 1)*((e + f*x)^(p + 1)/(d^2*(d 
*e - c*f)*(n + 1))), x] - Simp[1/(d^2*(d*e - c*f)*(n + 1))   Int[(c + d*x)^ 
(n + 1)*(e + f*x)^p*Simp[a^2*d^2*f*(n + p + 2) + b^2*c*(d*e*(n + 1) + c*f*( 
p + 1)) - 2*a*b*d*(d*e*(n + 1) + c*f*(p + 1)) - b^2*d*(d*e - c*f)*(n + 1)*x 
, x], x], x] /; FreeQ[{a, b, c, d, e, f, n, p}, x] && (LtQ[n, -1] || (EqQ[n 
 + p + 3, 0] && NeQ[n, -1] && (SumSimplerQ[n, 1] ||  !SumSimplerQ[p, 1])))
 
Maple [F]

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

Input:

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

Output:

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

Fricas [F]

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

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

Output:

integral((d^2*x^2 + 2*c*d*x + c^2)*(b*x + a)^n/x^2, x)
 

Sympy [A] (verification not implemented)

Time = 2.48 (sec) , antiderivative size = 265, normalized size of antiderivative = 3.05 \[ \int \frac {(a+b x)^n (c+d x)^2}{x^2} \, dx=d^{2} \left (\begin {cases} a^{n} x & \text {for}\: b = 0 \\\frac {\begin {cases} \frac {\left (a + b x\right )^{n + 1}}{n + 1} & \text {for}\: n \neq -1 \\\log {\left (a + b x \right )} & \text {otherwise} \end {cases}}{b} & \text {otherwise} \end {cases}\right ) - \frac {2 b^{n + 1} c d n \left (\frac {a}{b} + x\right )^{n + 1} \Phi \left (1 + \frac {b x}{a}, 1, n + 1\right ) \Gamma \left (n + 1\right )}{a \Gamma \left (n + 2\right )} - \frac {2 b^{n + 1} c d \left (\frac {a}{b} + x\right )^{n + 1} \Phi \left (1 + \frac {b x}{a}, 1, n + 1\right ) \Gamma \left (n + 1\right )}{a \Gamma \left (n + 2\right )} - \frac {b^{n + 2} c^{2} n \left (\frac {a}{b} + x\right )^{n + 1} \Gamma \left (n + 1\right )}{a b x \Gamma \left (n + 2\right )} - \frac {b^{n + 2} c^{2} \left (\frac {a}{b} + x\right )^{n + 1} \Gamma \left (n + 1\right )}{a b x \Gamma \left (n + 2\right )} - \frac {b^{n + 2} c^{2} n^{2} \left (\frac {a}{b} + x\right )^{n + 1} \Phi \left (1 + \frac {b x}{a}, 1, n + 1\right ) \Gamma \left (n + 1\right )}{a^{2} \Gamma \left (n + 2\right )} - \frac {b^{n + 2} c^{2} n \left (\frac {a}{b} + x\right )^{n + 1} \Phi \left (1 + \frac {b x}{a}, 1, n + 1\right ) \Gamma \left (n + 1\right )}{a^{2} \Gamma \left (n + 2\right )} \] Input:

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

Output:

d**2*Piecewise((a**n*x, Eq(b, 0)), (Piecewise(((a + b*x)**(n + 1)/(n + 1), 
 Ne(n, -1)), (log(a + b*x), True))/b, True)) - 2*b**(n + 1)*c*d*n*(a/b + x 
)**(n + 1)*lerchphi(1 + b*x/a, 1, n + 1)*gamma(n + 1)/(a*gamma(n + 2)) - 2 
*b**(n + 1)*c*d*(a/b + x)**(n + 1)*lerchphi(1 + b*x/a, 1, n + 1)*gamma(n + 
 1)/(a*gamma(n + 2)) - b**(n + 2)*c**2*n*(a/b + x)**(n + 1)*gamma(n + 1)/( 
a*b*x*gamma(n + 2)) - b**(n + 2)*c**2*(a/b + x)**(n + 1)*gamma(n + 1)/(a*b 
*x*gamma(n + 2)) - b**(n + 2)*c**2*n**2*(a/b + x)**(n + 1)*lerchphi(1 + b* 
x/a, 1, n + 1)*gamma(n + 1)/(a**2*gamma(n + 2)) - b**(n + 2)*c**2*n*(a/b + 
 x)**(n + 1)*lerchphi(1 + b*x/a, 1, n + 1)*gamma(n + 1)/(a**2*gamma(n + 2) 
)
 

Maxima [F]

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

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

Output:

(b*x + a)^(n + 1)*d^2/(b*(n + 1)) + integrate((2*c*d*x + c^2)*(b*x + a)^n/ 
x^2, x)
 

Giac [F]

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

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

Output:

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

Mupad [F(-1)]

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

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

Output:

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

Reduce [F]

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

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

Output:

((a + b*x)**n*a*d**2*n*x - (a + b*x)**n*b*c**2*n**2 - (a + b*x)**n*b*c**2* 
n + 2*(a + b*x)**n*b*c*d*n*x + 2*(a + b*x)**n*b*c*d*x + (a + b*x)**n*b*d** 
2*n*x**2 + 2*int((a + b*x)**n/(a*x + b*x**2),x)*a*b*c*d*n**2*x + 2*int((a 
+ b*x)**n/(a*x + b*x**2),x)*a*b*c*d*n*x + int((a + b*x)**n/(a*x + b*x**2), 
x)*b**2*c**2*n**3*x + int((a + b*x)**n/(a*x + b*x**2),x)*b**2*c**2*n**2*x) 
/(b*n*x*(n + 1))