dsolve(diff(y(x),x) = y(x)*(x-y(x))*(1+y(x))/x/(x*y(x)+x-y(x)),y(x))