[_linear]
Book solution method
Linear ODE
Mathematica ✓
cpu = 0.0606009 (sec), leaf count = 81
Maple ✓
cpu = 1.581 (sec), leaf count = 99
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)