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