DSolve[Derivative[1][y][x] == 1/(x*y[x]*(1 + x + x*y[x]^2)),y[x],x]