dsolve(diff(y(x),x) = y(x)*(x-y(x))/x/(x-y(x)^3),y(x))