sol = 
      NDSolve[{y''[x] + 2*(y[x]^2 + 1)*y'[x] + y[x] == 0, 
          y[0] == 0, y'[0] == 1}, y, {x, 0, 10}]

Plot[y[x] /. sol, {x, 0, 10}]