Las raíces del polinomio
$$ f(x)= x^8 + 8 x^6 + 60 x^4 + 176 x^2 + 2187 x + 1942 $$
permiten satisfacer las condiciones para $x \ne y \ne z$.
Hemos de redescubrir el 2 soluciones de $x=y=z=-1$$x=y=z=-2$, pero también una solución compleja
$(x,y,z)=(2.44944 -1.66645i,-1.74091+2.72124i, 0.791465+3.15828i) $
donde las rotaciones $(y,z,x), (z,x,y)$ y los complejos conjugados $(\overline x,\overline y, \overline z)$ son también soluciones.
Si escribimos $$g(x) = -{x^2+2\over 3}$$
a continuación, el conjunto de ecuaciones es también
$$ \begin{eqnarray} y=g(x)&;& z=g(y)=g(g(x))&;& x=g(z)=g(g(g(x))) \end{eqnarray}$$
y la ecuación de $x = g(g(g(x))))$ conduce a la polinomio $f(x)$ por encima.
He calculado las raíces utilizando Pari/GP y hacerlas llegar a los primeros dígitos:
$$ \pequeño
\begin{array} {rr}
-2.00000000000 \\
-1.00000000000 \\
2.44944037937&-1.66644590342*I \\
2.44944037937&+1.66644590342*I \\
-1.74090540769&-2.72123992391*I \\
-1.74090540769&+2.72123992391*I \\
0.791465028319&-3.15828086610*I \\
0.791465028319&+3.15828086610*I
\end{array}$$
Podemos reescribir el conjunto de ecuaciones
$$ \begin{eqnarray} y = -{x^2+2\over3} & ;&z = -{y^2+2\over3}&;&x = -{z^2+2\over3} \end{eqnarray}$$
Si puedo insertar las raíces que me pasa por las sucesivas evalutations la siguiente ( y última entrada debe ser igual que el primero de cada fila):
$$ \pequeño
\begin{array} {rr|rr|rr|rr}
x=root_r&& y=g(x) && z=g(y) && x=g(z) \\
\hline
-2.00000 && -2.00000 && -2.00000 && -2.00000 \\
-1.00000 && -1.00000 && -1.00000 && -1.00000 \\
2.44944 &-1.66645i & -1.74091&+2.72124i & 0.791465&+3.15828i & 2.44944&-1.66645i \\
2.44944 &+1.66645i & -1.74091&-2.72124i & 0.791465&-3.15828i & 2.44944&+1.66645i \\
-1.74091&-2.72124i & 0.791465&-3.15828i & 2.44944&+1.66645i & -1.74091&-2.72124i \\
-1.74091&+2.72124i & 0.791465&+3.15828i & 2.44944&-1.66645i & -1.74091&+2.72124i \\
0.791465&-3.15828i & 2.44944&+1.66645i & -1.74091&-2.72124i & 0.791465&-3.15828i \\
0.791465&+3.15828i & 2.44944&-1.66645i & -1.74091&+2.72124i & 0.791465&+3.15828i
\end{array} $$
Podemos expresar el 8 posibles soluciones de conjuntos de $(x,y,z)$ por los siguientes índices en el polroots, donde empezamos a contar en 1:
$[1,1,1],[2,2,2],[3,6,8],[4,5,7],[5,7,4],[6,8,3],[7,4,5],[8,3,6]$
Las dos primeras soluciones de $x=y=z$
mientras que en el siguiente 6 de soluciones de $x \ne y \ne z$ (siendo rotaciones y complejos conjugados de cada uno de los otros)