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