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