[_rational]
Book solution method
Exact equation, integrating factor
Mathematica ✓
cpu = 0.504613 (sec), leaf count = 23
Maple ✓
cpu = 0.092 (sec), leaf count = 26
DSolve[x*(x + y[x] + 2*y[x]^3)*y'[x] == (x - y[x])*y[x],y[x],x]
Mathematica raw output
Solve[C[1] + x/y[x] == Log[x] + Log[y[x]] + y[x]^2, y[x]]
Maple raw input
dsolve(x*(x+y(x)+2*y(x)^3)*diff(y(x),x) = (x-y(x))*y(x), y(x),'implicit')
Maple raw output
-ln(x)+x/y(x)-y(x)^2-ln(y(x))+_C1 = 0