DSolve[b + a*x*y[x] + x^2*(y[x]^2 + Derivative[1][y][x]) == 0,y[x],x]