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