[_linear]
Book solution method
Linear ODE
Mathematica ✓
cpu = 2.76097 (sec), leaf count = 54
Maple ✓
cpu = 0.111 (sec), leaf count = 41
DSolve[y'[x] == Sec[x]^2 + Csx[x]*Sec[x]*y[x],y[x],x]
Mathematica raw output
{{y[x] -> E^Integrate[Csx[K[1]]*Sec[K[1]], {K[1], 1, x}]*(C[1] + Integrate[Sec[K
[2]]^2/E^Integrate[Csx[K[1]]*Sec[K[1]], {K[1], 1, K[2]}], {K[2], 1, x}])}}
Maple raw input
dsolve(diff(y(x),x) = sec(x)^2+y(x)*sec(x)*Csx(x), y(x),'implicit')
Maple raw output
y(x) = (Int(2/(cos(2*x)+1)*exp(-Int(Csx(x)/cos(x),x)),x)+_C1)*exp(Int(Csx(x)/cos
(x),x))