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