4.12.49 (y(x)2+x)y(x)+y(x)=a+bx

ODE
(y(x)2+x)y(x)+y(x)=a+bx ODE Classification

[_exact, _rational]

Book solution method
Exact equation

Mathematica
cpu = 0.218048 (sec), leaf count = 420

{{y(x)23(9(2ax+bx2+2c1)2+16x3+6ax+3bx2+6c1)2/32 22/3x29(2ax+bx2+2c1)2+16x3+6ax+3bx2+6c13},{y(x)i(3+i)(36a2x2+36abx3+72ac1x+9b2x4+36bc1x2+36c12+16x3+6ax+3bx2+6c1)2/3+223(1+i3)x2 22/39(2ax+bx2+2c1)2+16x3+6ax+3bx2+6c13},{y(x)xi3x239(2ax+bx2+2c1)2+16x3+6ax+3bx2+6c13i(3i)9(2ax+bx2+2c1)2+16x3+6ax+3bx2+6c132 22/3}}

Maple
cpu = 0.014 (sec), leaf count = 25

{bx22ax+xy(x)+(y(x))33+_C1=0} Mathematica raw input

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