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