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