\(\int \frac {(e x)^{-1+2 n}}{a+b \csc (c+d x^n)} \, dx\) [80]

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

Optimal result

Integrand size = 24, antiderivative size = 338 \[ \int \frac {(e x)^{-1+2 n}}{a+b \csc \left (c+d x^n\right )} \, dx=\frac {(e x)^{2 n}}{2 a e n}+\frac {i b x^{-n} (e x)^{2 n} \log \left (1-\frac {i a e^{i \left (c+d x^n\right )}}{b-\sqrt {-a^2+b^2}}\right )}{a \sqrt {-a^2+b^2} d e n}-\frac {i b x^{-n} (e x)^{2 n} \log \left (1-\frac {i a e^{i \left (c+d x^n\right )}}{b+\sqrt {-a^2+b^2}}\right )}{a \sqrt {-a^2+b^2} d e n}+\frac {b x^{-2 n} (e x)^{2 n} \operatorname {PolyLog}\left (2,\frac {i a e^{i \left (c+d x^n\right )}}{b-\sqrt {-a^2+b^2}}\right )}{a \sqrt {-a^2+b^2} d^2 e n}-\frac {b x^{-2 n} (e x)^{2 n} \operatorname {PolyLog}\left (2,\frac {i a e^{i \left (c+d x^n\right )}}{b+\sqrt {-a^2+b^2}}\right )}{a \sqrt {-a^2+b^2} d^2 e n} \] Output:

1/2*(e*x)^(2*n)/a/e/n+I*b*(e*x)^(2*n)*ln(1-I*a*exp(I*(c+d*x^n))/(b-(-a^2+b 
^2)^(1/2)))/a/(-a^2+b^2)^(1/2)/d/e/n/(x^n)-I*b*(e*x)^(2*n)*ln(1-I*a*exp(I* 
(c+d*x^n))/(b+(-a^2+b^2)^(1/2)))/a/(-a^2+b^2)^(1/2)/d/e/n/(x^n)+b*(e*x)^(2 
*n)*polylog(2,I*a*exp(I*(c+d*x^n))/(b-(-a^2+b^2)^(1/2)))/a/(-a^2+b^2)^(1/2 
)/d^2/e/n/(x^(2*n))-b*(e*x)^(2*n)*polylog(2,I*a*exp(I*(c+d*x^n))/(b+(-a^2+ 
b^2)^(1/2)))/a/(-a^2+b^2)^(1/2)/d^2/e/n/(x^(2*n))
 

Mathematica [B] (warning: unable to verify)

Both result and optimal contain complex but leaf count is larger than twice the leaf count of optimal. \(1003\) vs. \(2(338)=676\).

Time = 3.94 (sec) , antiderivative size = 1003, normalized size of antiderivative = 2.97 \[ \int \frac {(e x)^{-1+2 n}}{a+b \csc \left (c+d x^n\right )} \, dx =\text {Too large to display} \] Input:

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

Output:

((e*x)^(2*n)*Csc[c + d*x^n]*(1 - (2*b*((Pi*ArcTan[(a + b*Tan[(c + d*x^n)/2 
])/Sqrt[-a^2 + b^2]])/Sqrt[-a^2 + b^2] + (2*(c - ArcCos[-(b/a)])*ArcTanh[( 
(a - b)*Cot[(2*c + Pi + 2*d*x^n)/4])/Sqrt[a^2 - b^2]] + (-2*c + Pi - 2*d*x 
^n)*ArcTanh[((a + b)*Tan[(2*c + Pi + 2*d*x^n)/4])/Sqrt[a^2 - b^2]] - (ArcC 
os[-(b/a)] - (2*I)*ArcTanh[((a - b)*Cot[(2*c + Pi + 2*d*x^n)/4])/Sqrt[a^2 
- b^2]])*Log[((a + b)*(a - b - I*Sqrt[a^2 - b^2])*(1 + I*Cot[(2*c + Pi + 2 
*d*x^n)/4]))/(a*(a + b + Sqrt[a^2 - b^2]*Cot[(2*c + Pi + 2*d*x^n)/4]))] + 
(ArcCos[-(b/a)] + (2*I)*(-ArcTanh[((a - b)*Cot[(2*c + Pi + 2*d*x^n)/4])/Sq 
rt[a^2 - b^2]] + ArcTanh[((a + b)*Tan[(2*c + Pi + 2*d*x^n)/4])/Sqrt[a^2 - 
b^2]]))*Log[((-1)^(1/4)*Sqrt[a^2 - b^2])/(Sqrt[2]*Sqrt[a]*E^((I/2)*(c + d* 
x^n))*Sqrt[b + a*Sin[c + d*x^n]])] + (ArcCos[-(b/a)] + (2*I)*ArcTanh[((a - 
 b)*Cot[(2*c + Pi + 2*d*x^n)/4])/Sqrt[a^2 - b^2]] - (2*I)*ArcTanh[((a + b) 
*Tan[(2*c + Pi + 2*d*x^n)/4])/Sqrt[a^2 - b^2]])*Log[-(((-1)^(3/4)*Sqrt[a^2 
 - b^2]*E^((I/2)*(c + d*x^n)))/(Sqrt[2]*Sqrt[a]*Sqrt[b + a*Sin[c + d*x^n]] 
))] - (ArcCos[-(b/a)] + (2*I)*ArcTanh[((a - b)*Cot[(2*c + Pi + 2*d*x^n)/4] 
)/Sqrt[a^2 - b^2]])*Log[1 + (I*(I*b + Sqrt[a^2 - b^2])*(a + b + Sqrt[a^2 - 
 b^2]*Tan[(2*c - Pi + 2*d*x^n)/4]))/(a*(a + b + Sqrt[a^2 - b^2]*Cot[(2*c + 
 Pi + 2*d*x^n)/4]))] + I*(PolyLog[2, ((b - I*Sqrt[a^2 - b^2])*(a + b + Sqr 
t[a^2 - b^2]*Tan[(2*c - Pi + 2*d*x^n)/4]))/(a*(a + b + Sqrt[a^2 - b^2]*Cot 
[(2*c + Pi + 2*d*x^n)/4]))] - PolyLog[2, ((b + I*Sqrt[a^2 - b^2])*(a + ...
 

Rubi [A] (verified)

Time = 0.85 (sec) , antiderivative size = 283, normalized size of antiderivative = 0.84, number of steps used = 6, number of rules used = 5, \(\frac {\text {number of rules}}{\text {integrand size}}\) = 0.208, Rules used = {4697, 4693, 3042, 4679, 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 \frac {(e x)^{2 n-1}}{a+b \csc \left (c+d x^n\right )} \, dx\)

\(\Big \downarrow \) 4697

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

\(\Big \downarrow \) 4693

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

\(\Big \downarrow \) 3042

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

\(\Big \downarrow \) 4679

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

\(\Big \downarrow \) 2009

\(\displaystyle \frac {x^{-2 n} (e x)^{2 n} \left (\frac {b \operatorname {PolyLog}\left (2,\frac {i a e^{i \left (d x^n+c\right )}}{b-\sqrt {b^2-a^2}}\right )}{a d^2 \sqrt {b^2-a^2}}-\frac {b \operatorname {PolyLog}\left (2,\frac {i a e^{i \left (d x^n+c\right )}}{b+\sqrt {b^2-a^2}}\right )}{a d^2 \sqrt {b^2-a^2}}+\frac {i b x^n \log \left (1-\frac {i a e^{i \left (c+d x^n\right )}}{b-\sqrt {b^2-a^2}}\right )}{a d \sqrt {b^2-a^2}}-\frac {i b x^n \log \left (1-\frac {i a e^{i \left (c+d x^n\right )}}{\sqrt {b^2-a^2}+b}\right )}{a d \sqrt {b^2-a^2}}+\frac {x^{2 n}}{2 a}\right )}{e n}\)

Input:

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

Output:

((e*x)^(2*n)*(x^(2*n)/(2*a) + (I*b*x^n*Log[1 - (I*a*E^(I*(c + d*x^n)))/(b 
- Sqrt[-a^2 + b^2])])/(a*Sqrt[-a^2 + b^2]*d) - (I*b*x^n*Log[1 - (I*a*E^(I* 
(c + d*x^n)))/(b + Sqrt[-a^2 + b^2])])/(a*Sqrt[-a^2 + b^2]*d) + (b*PolyLog 
[2, (I*a*E^(I*(c + d*x^n)))/(b - Sqrt[-a^2 + b^2])])/(a*Sqrt[-a^2 + b^2]*d 
^2) - (b*PolyLog[2, (I*a*E^(I*(c + d*x^n)))/(b + Sqrt[-a^2 + b^2])])/(a*Sq 
rt[-a^2 + b^2]*d^2)))/(e*n*x^(2*n))
 

Defintions of rubi rules used

rule 2009
Int[u_, x_Symbol] :> Simp[IntSum[u, x], x] /; SumQ[u]
 

rule 3042
Int[u_, x_Symbol] :> Int[DeactivateTrig[u, x], x] /; FunctionOfTrigOfLinear 
Q[u, x]
 

rule 4679
Int[(csc[(e_.) + (f_.)*(x_)]*(b_.) + (a_))^(n_.)*((c_.) + (d_.)*(x_))^(m_.) 
, x_Symbol] :> Int[ExpandIntegrand[(c + d*x)^m, 1/(Sin[e + f*x]^n/(b + a*Si 
n[e + f*x])^n), x], x] /; FreeQ[{a, b, c, d, e, f}, x] && ILtQ[n, 0] && IGt 
Q[m, 0]
 

rule 4693
Int[((a_.) + Csc[(c_.) + (d_.)*(x_)^(n_)]*(b_.))^(p_.)*(x_)^(m_.), x_Symbol 
] :> Simp[1/n   Subst[Int[x^(Simplify[(m + 1)/n] - 1)*(a + b*Csc[c + d*x])^ 
p, x], x, x^n], x] /; FreeQ[{a, b, c, d, m, n, p}, x] && IGtQ[Simplify[(m + 
 1)/n], 0] && IntegerQ[p]
 

rule 4697
Int[((a_.) + Csc[(c_.) + (d_.)*(x_)^(n_)]*(b_.))^(p_.)*((e_)*(x_))^(m_.), x 
_Symbol] :> Simp[e^IntPart[m]*((e*x)^FracPart[m]/x^FracPart[m])   Int[x^m*( 
a + b*Csc[c + d*x^n])^p, x], x] /; FreeQ[{a, b, c, d, e, m, n, p}, x]
 
Maple [C] (warning: unable to verify)

Result contains higher order function than in optimal. Order 9 vs. order 4.

Time = 0.42 (sec) , antiderivative size = 777, normalized size of antiderivative = 2.30

method result size
risch \(\frac {x \,{\mathrm e}^{\frac {\left (-1+2 n \right ) \left (-i \operatorname {csgn}\left (i e \right ) \operatorname {csgn}\left (i x \right ) \operatorname {csgn}\left (i e x \right ) \pi +i \operatorname {csgn}\left (i e \right ) \operatorname {csgn}\left (i e x \right )^{2} \pi +i \operatorname {csgn}\left (i x \right ) \operatorname {csgn}\left (i e x \right )^{2} \pi -i \operatorname {csgn}\left (i e x \right )^{3} \pi +2 \ln \left (x \right )+2 \ln \left (e \right )\right )}{2}}}{2 a n}+\frac {\left (-d \,x^{n} \ln \left (\frac {i {\mathrm e}^{i c} b +a \,{\mathrm e}^{i \left (d \,x^{n}+2 c \right )}-\sqrt {a^{2} {\mathrm e}^{2 i c}-{\mathrm e}^{2 i c} b^{2}}}{i {\mathrm e}^{i c} b -\sqrt {a^{2} {\mathrm e}^{2 i c}-{\mathrm e}^{2 i c} b^{2}}}\right )+d \,x^{n} \ln \left (\frac {i {\mathrm e}^{i c} b +a \,{\mathrm e}^{i \left (d \,x^{n}+2 c \right )}+\sqrt {a^{2} {\mathrm e}^{2 i c}-{\mathrm e}^{2 i c} b^{2}}}{i {\mathrm e}^{i c} b +\sqrt {a^{2} {\mathrm e}^{2 i c}-{\mathrm e}^{2 i c} b^{2}}}\right )+i \operatorname {dilog}\left (\frac {i b \,{\mathrm e}^{i c}}{i {\mathrm e}^{i c} b -\sqrt {a^{2} {\mathrm e}^{2 i c}-{\mathrm e}^{2 i c} b^{2}}}+\frac {a \,{\mathrm e}^{i \left (d \,x^{n}+2 c \right )}}{i {\mathrm e}^{i c} b -\sqrt {a^{2} {\mathrm e}^{2 i c}-{\mathrm e}^{2 i c} b^{2}}}-\frac {\sqrt {a^{2} {\mathrm e}^{2 i c}-{\mathrm e}^{2 i c} b^{2}}}{i {\mathrm e}^{i c} b -\sqrt {a^{2} {\mathrm e}^{2 i c}-{\mathrm e}^{2 i c} b^{2}}}\right )-i \operatorname {dilog}\left (\frac {i b \,{\mathrm e}^{i c}}{i {\mathrm e}^{i c} b +\sqrt {a^{2} {\mathrm e}^{2 i c}-{\mathrm e}^{2 i c} b^{2}}}+\frac {a \,{\mathrm e}^{i \left (d \,x^{n}+2 c \right )}}{i {\mathrm e}^{i c} b +\sqrt {a^{2} {\mathrm e}^{2 i c}-{\mathrm e}^{2 i c} b^{2}}}+\frac {\sqrt {a^{2} {\mathrm e}^{2 i c}-{\mathrm e}^{2 i c} b^{2}}}{i {\mathrm e}^{i c} b +\sqrt {a^{2} {\mathrm e}^{2 i c}-{\mathrm e}^{2 i c} b^{2}}}\right )\right ) \sqrt {a^{2} {\mathrm e}^{2 i c}-{\mathrm e}^{2 i c} b^{2}}\, e^{2 n} b \,{\mathrm e}^{-\frac {i \left (2 \pi n \,\operatorname {csgn}\left (i e \right ) \operatorname {csgn}\left (i x \right ) \operatorname {csgn}\left (i e x \right )-2 \pi n \,\operatorname {csgn}\left (i e \right ) \operatorname {csgn}\left (i e x \right )^{2}-2 \pi n \,\operatorname {csgn}\left (i x \right ) \operatorname {csgn}\left (i e x \right )^{2}+2 \pi n \operatorname {csgn}\left (i e x \right )^{3}-\pi \,\operatorname {csgn}\left (i e \right ) \operatorname {csgn}\left (i x \right ) \operatorname {csgn}\left (i e x \right )+\pi \,\operatorname {csgn}\left (i e \right ) \operatorname {csgn}\left (i e x \right )^{2}+\pi \,\operatorname {csgn}\left (i x \right ) \operatorname {csgn}\left (i e x \right )^{2}-\pi \operatorname {csgn}\left (i e x \right )^{3}+2 c \right )}{2}}}{\left (a^{2}-b^{2}\right ) d^{2} n e a}\) \(777\)

Input:

int((e*x)^(-1+2*n)/(a+b*csc(c+d*x^n)),x,method=_RETURNVERBOSE)
 

Output:

1/2/a/n*x*exp(1/2*(-1+2*n)*(-I*csgn(I*e)*csgn(I*x)*csgn(I*e*x)*Pi+I*csgn(I 
*e)*csgn(I*e*x)^2*Pi+I*csgn(I*x)*csgn(I*e*x)^2*Pi-I*csgn(I*e*x)^3*Pi+2*ln( 
x)+2*ln(e)))+1/(a^2-b^2)*(-d*x^n*ln((I*exp(I*c)*b+a*exp(I*(d*x^n+2*c))-(a^ 
2*exp(2*I*c)-exp(2*I*c)*b^2)^(1/2))/(I*exp(I*c)*b-(a^2*exp(2*I*c)-exp(2*I* 
c)*b^2)^(1/2)))+d*x^n*ln((I*exp(I*c)*b+a*exp(I*(d*x^n+2*c))+(a^2*exp(2*I*c 
)-exp(2*I*c)*b^2)^(1/2))/(I*exp(I*c)*b+(a^2*exp(2*I*c)-exp(2*I*c)*b^2)^(1/ 
2)))+I*dilog(I/(I*exp(I*c)*b-(a^2*exp(2*I*c)-exp(2*I*c)*b^2)^(1/2))*b*exp( 
I*c)+a/(I*exp(I*c)*b-(a^2*exp(2*I*c)-exp(2*I*c)*b^2)^(1/2))*exp(I*(d*x^n+2 
*c))-1/(I*exp(I*c)*b-(a^2*exp(2*I*c)-exp(2*I*c)*b^2)^(1/2))*(a^2*exp(2*I*c 
)-exp(2*I*c)*b^2)^(1/2))-I*dilog(I/(I*exp(I*c)*b+(a^2*exp(2*I*c)-exp(2*I*c 
)*b^2)^(1/2))*b*exp(I*c)+a/(I*exp(I*c)*b+(a^2*exp(2*I*c)-exp(2*I*c)*b^2)^( 
1/2))*exp(I*(d*x^n+2*c))+1/(I*exp(I*c)*b+(a^2*exp(2*I*c)-exp(2*I*c)*b^2)^( 
1/2))*(a^2*exp(2*I*c)-exp(2*I*c)*b^2)^(1/2)))*(a^2*exp(2*I*c)-exp(2*I*c)*b 
^2)^(1/2)/d^2/n/e*(e^n)^2/a*b*exp(-1/2*I*(2*Pi*n*csgn(I*e)*csgn(I*x)*csgn( 
I*e*x)-2*Pi*n*csgn(I*e)*csgn(I*e*x)^2-2*Pi*n*csgn(I*x)*csgn(I*e*x)^2+2*Pi* 
n*csgn(I*e*x)^3-Pi*csgn(I*e)*csgn(I*x)*csgn(I*e*x)+Pi*csgn(I*e)*csgn(I*e*x 
)^2+Pi*csgn(I*x)*csgn(I*e*x)^2-Pi*csgn(I*e*x)^3+2*c))
 

Fricas [B] (verification not implemented)

Both result and optimal contain complex but leaf count of result is larger than twice the leaf count of optimal. 1235 vs. \(2 (302) = 604\).

Time = 0.28 (sec) , antiderivative size = 1235, normalized size of antiderivative = 3.65 \[ \int \frac {(e x)^{-1+2 n}}{a+b \csc \left (c+d x^n\right )} \, dx=\text {Too large to display} \] Input:

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

Output:

-1/2*(a*b*c*e^(2*n - 1)*sqrt((a^2 - b^2)/a^2)*log(2*a*cos(d*x^n + c) + 2*I 
*a*sin(d*x^n + c) + 2*a*sqrt((a^2 - b^2)/a^2) + 2*I*b) + a*b*c*e^(2*n - 1) 
*sqrt((a^2 - b^2)/a^2)*log(2*a*cos(d*x^n + c) - 2*I*a*sin(d*x^n + c) + 2*a 
*sqrt((a^2 - b^2)/a^2) - 2*I*b) - a*b*c*e^(2*n - 1)*sqrt((a^2 - b^2)/a^2)* 
log(-2*a*cos(d*x^n + c) + 2*I*a*sin(d*x^n + c) + 2*a*sqrt((a^2 - b^2)/a^2) 
 + 2*I*b) - a*b*c*e^(2*n - 1)*sqrt((a^2 - b^2)/a^2)*log(-2*a*cos(d*x^n + c 
) - 2*I*a*sin(d*x^n + c) + 2*a*sqrt((a^2 - b^2)/a^2) - 2*I*b) - (a^2 - b^2 
)*d^2*e^(2*n - 1)*x^(2*n) - I*a*b*e^(2*n - 1)*sqrt((a^2 - b^2)/a^2)*dilog( 
((a*sqrt((a^2 - b^2)/a^2) + I*b)*cos(d*x^n + c) + (I*a*sqrt((a^2 - b^2)/a^ 
2) - b)*sin(d*x^n + c) - a)/a + 1) - I*a*b*e^(2*n - 1)*sqrt((a^2 - b^2)/a^ 
2)*dilog(-((a*sqrt((a^2 - b^2)/a^2) + I*b)*cos(d*x^n + c) - (I*a*sqrt((a^2 
 - b^2)/a^2) - b)*sin(d*x^n + c) + a)/a + 1) + I*a*b*e^(2*n - 1)*sqrt((a^2 
 - b^2)/a^2)*dilog(((a*sqrt((a^2 - b^2)/a^2) - I*b)*cos(d*x^n + c) + (-I*a 
*sqrt((a^2 - b^2)/a^2) - b)*sin(d*x^n + c) - a)/a + 1) + I*a*b*e^(2*n - 1) 
*sqrt((a^2 - b^2)/a^2)*dilog(-((a*sqrt((a^2 - b^2)/a^2) - I*b)*cos(d*x^n + 
 c) - (-I*a*sqrt((a^2 - b^2)/a^2) - b)*sin(d*x^n + c) + a)/a + 1) + (a*b*d 
*e^(2*n - 1)*x^n*sqrt((a^2 - b^2)/a^2) + a*b*c*e^(2*n - 1)*sqrt((a^2 - b^2 
)/a^2))*log(-((a*sqrt((a^2 - b^2)/a^2) + I*b)*cos(d*x^n + c) + (I*a*sqrt(( 
a^2 - b^2)/a^2) - b)*sin(d*x^n + c) - a)/a) - (a*b*d*e^(2*n - 1)*x^n*sqrt( 
(a^2 - b^2)/a^2) + a*b*c*e^(2*n - 1)*sqrt((a^2 - b^2)/a^2))*log(((a*sqr...
 

Sympy [F]

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

integrate((e*x)**(-1+2*n)/(a+b*csc(c+d*x**n)),x)
 

Output:

Integral((e*x)**(2*n - 1)/(a + b*csc(c + d*x**n)), x)
 

Maxima [F]

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

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

Output:

-1/2*(4*a*b*e^(2*n + 1)*n*integrate((2*b*x^(2*n)*cos(d*x^n + c)^2 + a*x^(2 
*n)*cos(d*x^n + c)*sin(2*d*x^n + 2*c) - a*x^(2*n)*cos(2*d*x^n + 2*c)*sin(d 
*x^n + c) + 2*b*x^(2*n)*sin(d*x^n + c)^2 + a*x^(2*n)*sin(d*x^n + c))/(a^3* 
e*x*cos(2*d*x^n + 2*c)^2 + 4*a*b^2*e*x*cos(d*x^n + c)^2 + 4*a^2*b*e*x*cos( 
d*x^n + c)*sin(2*d*x^n + 2*c) + a^3*e*x*sin(2*d*x^n + 2*c)^2 + 4*a*b^2*e*x 
*sin(d*x^n + c)^2 + 4*a^2*b*e*x*sin(d*x^n + c) + a^3*e*x - 2*(2*a^2*b*e*x* 
sin(d*x^n + c) + a^3*e*x)*cos(2*d*x^n + 2*c)), x) - e^(2*n)*x^(2*n))/(a*e* 
n)
 

Giac [F]

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

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

Output:

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

Mupad [F(-1)]

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

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

Output:

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

Reduce [F]

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

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

Output:

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