[_exact, _rational]
Book solution method
Exact equation
Mathematica ✓
cpu = 0.218048 (sec), leaf count = 420
Maple ✓
cpu = 0.014 (sec), leaf count = 25
DSolve[y[x] + (x + y[x]^2)*y'[x] == a + b*x,y[x],x]
Mathematica raw output
{{y[x] -> (-2*2^(2/3)*x + 2^(1/3)*(6*a*x + 3*b*x^2 + 6*C[1] + Sqrt[16*x^3 + 9*(2
*a*x + b*x^2 + 2*C[1])^2])^(2/3))/(2*(6*a*x + 3*b*x^2 + 6*C[1] + Sqrt[16*x^3 + 9
*(2*a*x + b*x^2 + 2*C[1])^2])^(1/3))}, {y[x] -> (2*2^(1/3)*(1 + I*Sqrt[3])*x + I
*(I + Sqrt[3])*(6*a*x + 3*b*x^2 + 6*C[1] + Sqrt[36*a^2*x^2 + 16*x^3 + 36*a*b*x^3
+ 9*b^2*x^4 + 72*a*x*C[1] + 36*b*x^2*C[1] + 36*C[1]^2])^(2/3))/(2*2^(2/3)*(6*a*
x + 3*b*x^2 + 6*C[1] + Sqrt[16*x^3 + 9*(2*a*x + b*x^2 + 2*C[1])^2])^(1/3))}, {y[
x] -> (x - I*Sqrt[3]*x)/(2^(1/3)*(6*a*x + 3*b*x^2 + 6*C[1] + Sqrt[16*x^3 + 9*(2*
a*x + b*x^2 + 2*C[1])^2])^(1/3)) - ((I/2)*(-I + Sqrt[3])*(6*a*x + 3*b*x^2 + 6*C[
1] + Sqrt[16*x^3 + 9*(2*a*x + b*x^2 + 2*C[1])^2])^(1/3))/2^(2/3)}}
Maple raw input
dsolve((x+y(x)^2)*diff(y(x),x)+y(x) = b*x+a, y(x),'implicit')
Maple raw output
-1/2*b*x^2-a*x+x*y(x)+1/3*y(x)^3+_C1 = 0