Selbstverständlich kannst du hier den Satz über implizite Funktionen anwenden.
Betrachte dazu
$$F(x,y,u)= \begin{pmatrix} f_1(x,y,u) \\ f_2(x,y,u)\end{pmatrix} = \begin{pmatrix} \sin(x+u) -e^y + 1 \\ x^2+y+e^u - 1\end{pmatrix}$$
Für \((x_0,y_0,u_0) = (0,0,0)\) gilt offensichtlich
$$F(0,0,0) = \begin{pmatrix} 0\\ 0\end{pmatrix}$$
Jetzt musst du nur noch die Jakobi-Matrix an der Stelle \((0,0,0)\) betrachten:
$$J_F(0,0,0) = \begin{pmatrix} 1 & {\color{blue}{-1}} & {\color{blue}{1}} \\ 0 & {\color{blue}{1}} & {\color{blue}{1}}\end{pmatrix}$$
Die blau hervorgehobene Submatrix hat eine von null verschiedene Determinante. Also gibt es laut Satz über implizite Funktionen in einer Umgebung von \(x_0 = 0\) eine Funktion
\(g(x) = \begin{pmatrix} y(x)\\ u(x) \end{pmatrix}\) mit \(g(0) = \begin{pmatrix} 0\\ 0 \end{pmatrix}\) und \(F(x,g(x)) = 0\).
Damit hat das System unendlich viele Lösungen.