DSolve[Derivative[1][y][x] == -b/2 - (a*x)/2 + x*Sqrt[b^2 - 4*c + 2*a*b*x + a^2*x^2 + 4*a*y[x]],y[x],x]