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