dsolve(diff(y(x),x) = -1/2*a*x-1/2*b+x*(a^2*x^2+2*a*b*x+b^2+4*a*y(x)-4*c)^(1/2),y(x))