[_linear]
Book solution method
Linear ODE
Mathematica ✓
cpu = 28.6501 (sec), leaf count = 155
Maple ✓
cpu = 1.476 (sec), leaf count = 125
DSolve[y'[x] == 4*x*Csc[x]*(1 - Tan[x]^2 + y[x]),y[x],x]
Mathematica raw output
{{y[x] -> E^(4*x*(Log[1 - E^(I*x)] - Log[1 + E^(I*x)]) + (4*I)*PolyLog[2, -E^(I*
x)] - (4*I)*PolyLog[2, E^(I*x)])*(C[1] + Integrate[4*E^(4*K[1]*(-Log[1 - E^(I*K[
1])] + Log[1 + E^(I*K[1])]) - (4*I)*PolyLog[2, -E^(I*K[1])] + (4*I)*PolyLog[2, E
^(I*K[1])])*Cos[2*K[1]]*Csc[K[1]]*K[1]*Sec[K[1]]^2, {K[1], 1, x}])}}
Maple raw input
dsolve(diff(y(x),x) = 4*csc(x)*x*(1-tan(x)^2+y(x)), y(x),'implicit')
Maple raw output
y(x) = exp(-4*I*(polylog(2,exp(I*x))-polylog(2,-exp(I*x))))*(1-exp(I*x))^(4*x)*(
1+exp(I*x))^(-4*x)*(4*Int(4*x*(1-exp(I*x))^(-4*x)*(1+exp(I*x))^(4*x)*exp(4*I*(po
lylog(2,exp(I*x))-polylog(2,-exp(I*x))))*(-sin(3*x)+sin(x))/(-1+cos(4*x)),x)+_C1
)