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