Da in der zweiten Gleichung das y alleine steht kannst du direkt nach y auflösen
3y^2 - 3 = 0 --> y = ±1
Für die dritte Gleichung wende den Satz vom Nullprodukt an
2·z - 2·x·z = 2·z·(1 - x) = 0 --> x = 1 ∨ z = 0
Mit diesen 2 Fällen gehst du in die letzte Gleichung
Fall 1: x = 1
4·1 + z^2 = 0 --> z = ± 2·i oder keine Lösung in R
Fall 2: z = 0
4·x + 0^2 = 0 → x = 0
Damit hast du alle Lösungen beisammen
[x = 0 ∧ y = 1 ∧ z = 0,
x = 0 ∧ y = -1 ∧ z = 0,
x = 1 ∧ y = 1 ∧ z = 2·i,
x = 1 ∧ y = 1 ∧ z = - 2·i,
x = 1 ∧ y = -1 ∧ z = 2·i,
x = 1 ∧ y = -1 ∧ z = - 2·i]
Wie gesagt im Bereich der reellen Zahlen entfallen die Komplexen Lösungen.