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

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

Optimal result

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

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

Mathematica [A] (verified)

Time = 3.26 (sec) , antiderivative size = 181, normalized size of antiderivative = 1.16 \[ \int \frac {\cot ^2\left (d \left (a+b \log \left (c x^n\right )\right )\right )}{x^2} \, dx=\frac {e^{2 i d \left (a+b \log \left (c x^n\right )\right )} \operatorname {Hypergeometric2F1}\left (1,1+\frac {i}{2 b d n},2+\frac {i}{2 b d n},e^{2 i d \left (a+b \log \left (c x^n\right )\right )}\right )+(i+2 b d n) \left (b d n-\cot \left (d \left (a+b \log \left (c x^n\right )\right )\right )-i \operatorname {Hypergeometric2F1}\left (1,\frac {i}{2 b d n},1+\frac {i}{2 b d n},e^{2 i d \left (a+b \log \left (c x^n\right )\right )}\right )\right )}{b d n (i+2 b d n) x} \] Input:

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

Output:

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

Rubi [A] (verified)

Time = 0.45 (sec) , antiderivative size = 211, normalized size of antiderivative = 1.35, number of steps used = 7, number of rules used = 6, \(\frac {\text {number of rules}}{\text {integrand size}}\) = 0.316, Rules used = {5009, 5007, 1004, 27, 959, 888}

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

\(\Big \downarrow \) 5009

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

\(\Big \downarrow \) 5007

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

\(\Big \downarrow \) 1004

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

\(\Big \downarrow \) 27

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

\(\Big \downarrow \) 959

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

\(\Big \downarrow \) 888

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

Input:

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

Output:

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

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 888
Int[((c_.)*(x_))^(m_.)*((a_) + (b_.)*(x_)^(n_))^(p_), x_Symbol] :> Simp[a^p 
*((c*x)^(m + 1)/(c*(m + 1)))*Hypergeometric2F1[-p, (m + 1)/n, (m + 1)/n + 1 
, (-b)*(x^n/a)], x] /; FreeQ[{a, b, c, m, n, p}, x] &&  !IGtQ[p, 0] && (ILt 
Q[p, 0] || GtQ[a, 0])
 

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

rule 1004
Int[((e_.)*(x_))^(m_.)*((a_) + (b_.)*(x_)^(n_))^(p_)*((c_) + (d_.)*(x_)^(n_ 
))^(q_), x_Symbol] :> Simp[(-(c*b - a*d))*(e*x)^(m + 1)*(a + b*x^n)^(p + 1) 
*((c + d*x^n)^(q - 1)/(a*b*e*n*(p + 1))), x] + Simp[1/(a*b*n*(p + 1))   Int 
[(e*x)^m*(a + b*x^n)^(p + 1)*(c + d*x^n)^(q - 2)*Simp[c*(c*b*n*(p + 1) + (c 
*b - a*d)*(m + 1)) + d*(c*b*n*(p + 1) + (c*b - a*d)*(m + n*(q - 1) + 1))*x^ 
n, x], x], x] /; FreeQ[{a, b, c, d, e, m, n}, x] && NeQ[b*c - a*d, 0] && Lt 
Q[p, -1] && GtQ[q, 1] && IntBinomialQ[a, b, c, d, e, m, n, p, q, x]
 

rule 5007
Int[Cot[((a_.) + Log[x_]*(b_.))*(d_.)]^(p_.)*((e_.)*(x_))^(m_.), x_Symbol] 
:> Int[(e*x)^m*((-I - I*E^(2*I*a*d)*x^(2*I*b*d))/(1 - E^(2*I*a*d)*x^(2*I*b* 
d)))^p, x] /; FreeQ[{a, b, d, e, m, p}, x]
 

rule 5009
Int[Cot[((a_.) + Log[(c_.)*(x_)^(n_.)]*(b_.))*(d_.)]^(p_.)*((e_.)*(x_))^(m_ 
.), x_Symbol] :> Simp[(e*x)^(m + 1)/(e*n*(c*x^n)^((m + 1)/n))   Subst[Int[x 
^((m + 1)/n - 1)*Cot[d*(a + b*Log[x])]^p, x], x, c*x^n], x] /; FreeQ[{a, b, 
 c, d, e, m, n, p}, x] && (NeQ[c, 1] || NeQ[n, 1])
 
Maple [F]

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

Input:

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

Output:

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

Fricas [F]

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

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

Output:

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

Sympy [F]

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

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

Output:

Integral(cot(a*d + b*d*log(c*x**n))**2/x**2, x)
 

Maxima [F]

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

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

Output:

-((b*d*cos(2*b*d*log(c))^2 + b*d*sin(2*b*d*log(c))^2)*n*cos(2*b*d*log(x^n) 
 + 2*a*d)^2 + (b*d*cos(2*b*d*log(c))^2 + b*d*sin(2*b*d*log(c))^2)*n*sin(2* 
b*d*log(x^n) + 2*a*d)^2 + b*d*n - 2*(b*d*n*cos(2*b*d*log(c)) + sin(2*b*d*l 
og(c)))*cos(2*b*d*log(x^n) + 2*a*d) - (2*b^2*d^2*n^2*x*cos(2*b*d*log(c))*c 
os(2*b*d*log(x^n) + 2*a*d) - 2*b^2*d^2*n^2*x*sin(2*b*d*log(c))*sin(2*b*d*l 
og(x^n) + 2*a*d) - b^2*d^2*n^2*x - (b^2*d^2*cos(2*b*d*log(c))^2 + b^2*d^2* 
sin(2*b*d*log(c))^2)*n^2*x*cos(2*b*d*log(x^n) + 2*a*d)^2 - (b^2*d^2*cos(2* 
b*d*log(c))^2 + b^2*d^2*sin(2*b*d*log(c))^2)*n^2*x*sin(2*b*d*log(x^n) + 2* 
a*d)^2)*integrate((cos(b*d*log(x^n) + a*d)*sin(b*d*log(c)) + cos(b*d*log(c 
))*sin(b*d*log(x^n) + a*d))/(2*b^2*d^2*n^2*x^2*cos(b*d*log(c))*cos(b*d*log 
(x^n) + a*d) - 2*b^2*d^2*n^2*x^2*sin(b*d*log(c))*sin(b*d*log(x^n) + a*d) + 
 b^2*d^2*n^2*x^2 + (b^2*d^2*cos(b*d*log(c))^2 + b^2*d^2*sin(b*d*log(c))^2) 
*n^2*x^2*cos(b*d*log(x^n) + a*d)^2 + (b^2*d^2*cos(b*d*log(c))^2 + b^2*d^2* 
sin(b*d*log(c))^2)*n^2*x^2*sin(b*d*log(x^n) + a*d)^2), x) + (2*b^2*d^2*n^2 
*x*cos(2*b*d*log(c))*cos(2*b*d*log(x^n) + 2*a*d) - 2*b^2*d^2*n^2*x*sin(2*b 
*d*log(c))*sin(2*b*d*log(x^n) + 2*a*d) - b^2*d^2*n^2*x - (b^2*d^2*cos(2*b* 
d*log(c))^2 + b^2*d^2*sin(2*b*d*log(c))^2)*n^2*x*cos(2*b*d*log(x^n) + 2*a* 
d)^2 - (b^2*d^2*cos(2*b*d*log(c))^2 + b^2*d^2*sin(2*b*d*log(c))^2)*n^2*x*s 
in(2*b*d*log(x^n) + 2*a*d)^2)*integrate(-(cos(b*d*log(x^n) + a*d)*sin(b*d* 
log(c)) + cos(b*d*log(c))*sin(b*d*log(x^n) + a*d))/(2*b^2*d^2*n^2*x^2*c...
 

Giac [F(-1)]

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

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

Output:

Timed out
 

Mupad [F(-1)]

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

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

Output:

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

Reduce [F]

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

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

Output:

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