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