Aloha :)
Deine Lösungen sind richtig, aber es gibt noch eine dritte Lösung. Vermutlich bist du an die Aufgabe mit Lagrange-Multiplikatoren rangegangen. Das ist das Standard-Verfahren, leider unnötig kompliziert und fummelig zu handhaben. Das ist quasi dafür gemacht, dass man Lösungen übersieht. Besser machst du es mit Determinaten:$$f(x,y,z)=(x-1)^2+y^2+z^2\quad;\quad g(x,y,z)=x+y-z^2+1=0$$Wir bilden die beiden Gradienten und schreiben sie als Spalten in eine \(3\times2\)-Matrix. Da wir davon die Determinante nicht bestimmen können, streichen wir einfach eine Zeile, und bestimmen 2 Determinanten von 2 quadratischen Matrizen:
$$0\stackrel{!}{=}\left|\begin{array}{c}2(x-1) & 1\\\not2\not y & \not1 \\2 z & -2z\end{array}\right|=\left|\begin{array}{c}2(x-1) & 1\\2 z & -2z\end{array}\right|=\left|\begin{array}{c}2x-1 & 1\\0 & -2z\end{array}\right|=-2z(2x-1)$$$$\qquad\Rightarrow\quad z=0\quad\lor\quad x=\frac{1}{2}$$$$0\stackrel{!}{=}\left|\begin{array}{c}2(x-1) & 1\\2 y & 1 \\\not2 \not z & \not-\not2\not z\end{array}\right|=\left|\begin{array}{c}2(x-1) & 1\\2y & 1\end{array}\right|=2(x-1)-2y=2(x-y)-2$$$$\qquad\Rightarrow\quad x-y=1\quad\Rightarrow\quad y=x-1$$
Du hast nur den Fall \(x=\frac{1}{2}\) und \(y=-\frac{1}{2}\) weiter betrachtet. Setzt man diese Werte in die Nebenbedinung ein, erhält man die fehlenden Werte für \(z\):$$0=x+y-z^2+1=\frac{1}{2}-\frac{1}{2}-z^2+1=-z^2+1\quad\Leftrightarrow\quad z=\pm1$$Das sind deine beiden kritischen Punkte:$$P_1\left(\frac{1}{2}\,;\,-\frac{1}{2}\,;\,1\right)\quad;\quad P_2\left(\frac{1}{2}\,;\,-\frac{1}{2}\,;\,-1\right)$$Übersehen hast du bei der Lagrange-Fummelei die Lösung \(z=0\):$$0\stackrel{!}{=}g(x;x-1;0)=x+(x-1)-0^2+1=2x\quad\Rightarrow\quad x=0\;;\;y=x-1=-1$$Die dritte Lösung ist daher:$$P_3(0\,;\,-1\,;\,0)$$