\(\int \frac {(e x)^m (a+b x)}{(a c-b c x)^2} \, dx\) [319]

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

Optimal result

Integrand size = 22, antiderivative size = 71 \[ \int \frac {(e x)^m (a+b x)}{(a c-b c x)^2} \, dx=-\frac {(e x)^{1+m}}{c^2 e m (a-b x)}+\frac {(1+2 m) (e x)^{1+m} \operatorname {Hypergeometric2F1}\left (2,1+m,2+m,\frac {b x}{a}\right )}{a c^2 e m (1+m)} \] Output:

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

Mathematica [A] (verified)

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

Integrate[((e*x)^m*(a + b*x))/(a*c - b*c*x)^2,x]
 

Output:

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

Rubi [A] (verified)

Time = 0.17 (sec) , antiderivative size = 66, normalized size of antiderivative = 0.93, number of steps used = 3, number of rules used = 3, \(\frac {\text {number of rules}}{\text {integrand size}}\) = 0.136, Rules used = {87, 27, 74}

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

\(\Big \downarrow \) 87

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

\(\Big \downarrow \) 27

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

\(\Big \downarrow \) 74

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

Input:

Int[((e*x)^m*(a + b*x))/(a*c - b*c*x)^2,x]
 

Output:

(2*(e*x)^(1 + m))/(c^2*e*(a - b*x)) - ((1 + 2*m)*(e*x)^(1 + m)*Hypergeomet 
ric2F1[1, 1 + m, 2 + m, (b*x)/a])/(a*c^2*e*(1 + m))
 

Defintions of rubi rules used

rule 27
Int[(a_)*(Fx_), x_Symbol] :> Simp[a   Int[Fx, x], x] /; FreeQ[a, x] &&  !Ma 
tchQ[Fx, (b_)*(Gx_) /; FreeQ[b, x]]
 

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

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

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

Input:

int((e*x)^m*(b*x+a)/(-b*c*x+a*c)^2,x)
 

Output:

int((e*x)^m*(b*x+a)/(-b*c*x+a*c)^2,x)
 

Fricas [F]

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

integrate((e*x)^m*(b*x+a)/(-b*c*x+a*c)^2,x, algorithm="fricas")
 

Output:

integral((b*x + a)*(e*x)^m/(b^2*c^2*x^2 - 2*a*b*c^2*x + a^2*c^2), x)
 

Sympy [C] (verification not implemented)

Result contains complex when optimal does not.

Time = 2.74 (sec) , antiderivative size = 790, normalized size of antiderivative = 11.13 \[ \int \frac {(e x)^m (a+b x)}{(a c-b c x)^2} \, dx =\text {Too large to display} \] Input:

integrate((e*x)**m*(b*x+a)/(-b*c*x+a*c)**2,x)
 

Output:

a*(a*e**m*m**2*x**(m + 1)*lerchphi(b*x*exp_polar(2*I*pi)/a, 1, m + 1)*gamm 
a(m + 1)/(-a**3*c**2*gamma(m + 2) + a**2*b*c**2*x*gamma(m + 2)) + a*e**m*m 
*x**(m + 1)*lerchphi(b*x*exp_polar(2*I*pi)/a, 1, m + 1)*gamma(m + 1)/(-a** 
3*c**2*gamma(m + 2) + a**2*b*c**2*x*gamma(m + 2)) - a*e**m*m*x**(m + 1)*ga 
mma(m + 1)/(-a**3*c**2*gamma(m + 2) + a**2*b*c**2*x*gamma(m + 2)) - a*e**m 
*x**(m + 1)*gamma(m + 1)/(-a**3*c**2*gamma(m + 2) + a**2*b*c**2*x*gamma(m 
+ 2)) - b*e**m*m**2*x*x**(m + 1)*lerchphi(b*x*exp_polar(2*I*pi)/a, 1, m + 
1)*gamma(m + 1)/(-a**3*c**2*gamma(m + 2) + a**2*b*c**2*x*gamma(m + 2)) - b 
*e**m*m*x*x**(m + 1)*lerchphi(b*x*exp_polar(2*I*pi)/a, 1, m + 1)*gamma(m + 
 1)/(-a**3*c**2*gamma(m + 2) + a**2*b*c**2*x*gamma(m + 2))) + b*(a*e**m*m* 
*2*x**(m + 2)*lerchphi(b*x*exp_polar(2*I*pi)/a, 1, m + 2)*gamma(m + 2)/(-a 
**3*c**2*gamma(m + 3) + a**2*b*c**2*x*gamma(m + 3)) + 3*a*e**m*m*x**(m + 2 
)*lerchphi(b*x*exp_polar(2*I*pi)/a, 1, m + 2)*gamma(m + 2)/(-a**3*c**2*gam 
ma(m + 3) + a**2*b*c**2*x*gamma(m + 3)) - a*e**m*m*x**(m + 2)*gamma(m + 2) 
/(-a**3*c**2*gamma(m + 3) + a**2*b*c**2*x*gamma(m + 3)) + 2*a*e**m*x**(m + 
 2)*lerchphi(b*x*exp_polar(2*I*pi)/a, 1, m + 2)*gamma(m + 2)/(-a**3*c**2*g 
amma(m + 3) + a**2*b*c**2*x*gamma(m + 3)) - 2*a*e**m*x**(m + 2)*gamma(m + 
2)/(-a**3*c**2*gamma(m + 3) + a**2*b*c**2*x*gamma(m + 3)) - b*e**m*m**2*x* 
x**(m + 2)*lerchphi(b*x*exp_polar(2*I*pi)/a, 1, m + 2)*gamma(m + 2)/(-a**3 
*c**2*gamma(m + 3) + a**2*b*c**2*x*gamma(m + 3)) - 3*b*e**m*m*x*x**(m +...
                                                                                    
                                                                                    
 

Maxima [F]

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

integrate((e*x)^m*(b*x+a)/(-b*c*x+a*c)^2,x, algorithm="maxima")
 

Output:

integrate((b*x + a)*(e*x)^m/(b*c*x - a*c)^2, x)
 

Giac [F]

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

integrate((e*x)^m*(b*x+a)/(-b*c*x+a*c)^2,x, algorithm="giac")
 

Output:

integrate((b*x + a)*(e*x)^m/(b*c*x - a*c)^2, x)
 

Mupad [F(-1)]

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

int(((e*x)^m*(a + b*x))/(a*c - b*c*x)^2,x)
 

Output:

int(((e*x)^m*(a + b*x))/(a*c - b*c*x)^2, x)
 

Reduce [F]

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

int((e*x)^m*(b*x+a)/(-b*c*x+a*c)^2,x)
 

Output:

(e**m*( - 2*x**m*a*m - x**m*a - x**m*b*m*x + x**m*b*x + 2*int(x**m/(a**2*m 
*x - a**2*x - 2*a*b*m*x**2 + 2*a*b*x**2 + b**2*m*x**3 - b**2*x**3),x)*a**3 
*m**3 - int(x**m/(a**2*m*x - a**2*x - 2*a*b*m*x**2 + 2*a*b*x**2 + b**2*m*x 
**3 - b**2*x**3),x)*a**3*m**2 - int(x**m/(a**2*m*x - a**2*x - 2*a*b*m*x**2 
 + 2*a*b*x**2 + b**2*m*x**3 - b**2*x**3),x)*a**3*m - 2*int(x**m/(a**2*m*x 
- a**2*x - 2*a*b*m*x**2 + 2*a*b*x**2 + b**2*m*x**3 - b**2*x**3),x)*a**2*b* 
m**3*x + int(x**m/(a**2*m*x - a**2*x - 2*a*b*m*x**2 + 2*a*b*x**2 + b**2*m* 
x**3 - b**2*x**3),x)*a**2*b*m**2*x + int(x**m/(a**2*m*x - a**2*x - 2*a*b*m 
*x**2 + 2*a*b*x**2 + b**2*m*x**3 - b**2*x**3),x)*a**2*b*m*x))/(b*c**2*m*(a 
*m - a - b*m*x + b*x))