[_rational]
Book solution method
Exact equation, integrating factor
Mathematica ✓
cpu = 0.0252563 (sec), leaf count = 16
Maple ✓
cpu = 0.046 (sec), leaf count = 31
DSolve[(x + x^2 + y[x]^2)*y'[x] == y[x],y[x],x]
Mathematica raw output
Solve[ArcTan[x/y[x]] + C[1] == y[x], y[x]]
Maple raw input
dsolve((x+x^2+y(x)^2)*diff(y(x),x) = y(x), y(x),'implicit')
Maple raw output
_C1-exp(-2*I*y(x))*(I*x+y(x))/(2*I*y(x)+2*x) = 0