4.1.21 y(x)=4xcsc(x)sec2(x)2y(x)cot(2x)

ODE
y(x)=4xcsc(x)sec2(x)2y(x)cot(2x) ODE Classification

[_linear]

Book solution method
Linear ODE

Mathematica
cpu = 0.0606009 (sec), leaf count = 81

{{y(x)csc(x)sec(x)(c1+4iLi2(ieix)4iLi2(ieix)+4xlog(1ieix)4xlog(1+ieix))}}

Maple
cpu = 1.581 (sec), leaf count = 99

{y(x)=16(sin(2x)(idilog(1ieix)idilog(1+ieix)xln(1ieix)+xln(1+ieix))e4ix(e4ix1)2_C1/16)(cot(2x))2+1} Mathematica raw input

DSolve[y'[x] == 4*x*Csc[x]*Sec[x]^2 - 2*Cot[2*x]*y[x],y[x],x]

Mathematica raw output

{{y[x] -> Csc[x]*(C[1] + 4*x*Log[1 - I*E^(I*x)] - 4*x*Log[1 + I*E^(I*x)] + (4*I)
*PolyLog[2, (-I)*E^(I*x)] - (4*I)*PolyLog[2, I*E^(I*x)])*Sec[x]}}

Maple raw input

dsolve(diff(y(x),x) = 4*csc(x)*x*sec(x)^2-2*y(x)*cot(2*x), y(x),'implicit')

Maple raw output

y(x) = -16*(sin(2*x)*(I*dilog(1-I*exp(I*x))-I*dilog(1+I*exp(I*x))-x*ln(1-I*exp(I
*x))+x*ln(1+I*exp(I*x)))*(-exp(4*I*x)/(exp(4*I*x)-1)^2)^(1/2)-1/16*_C1)*(cot(2*x
)^2+1)^(1/2)