dsolve(diff(y(x),x) = (-b^3+6*b^2*x-12*b*x^2+8*x^3-4*b*y(x)^2+8*x*y(x)^2+8*y(x)^3)/(2*x-b)^3,y(x))