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