Hallo :-)
Du musst zunächst die Jacobi-Matrix \(J_f\) von \(f\) bestimmen:
$$ J_f((x,y))=\begin{pmatrix}0&2y\\1&1\end{pmatrix}\in \R^{2,2} $$
Und dann kannst du schon loslegen:
Startvektor: \(v_0=\begin{pmatrix}0\\2\end{pmatrix}\)
Iterationen: Für \(k=0,1,...\) bis Abbruchkriterium (kann man sich beliebig festlegen; hier: nur zwei Iterationen) erfüllt ist mache:
1.) Löse das LGS \(J_f(v_k)\cdot d_k=-f(v_k) \).
2.) Setze \(v_{k+1}=v_k+d_k\).
Du kannst natürlich auch 1.) und 2.) zu einem Schritt zusammenfassen, indem du die Inverse von \(J_f\) bestimmst. Dann hast du nämlich \(v_{k+1}=v_k-(J_f(v_k))^{-1}\cdot f(v_k)\). Der Vorteil ist, dass du nicht in jeder Iteration ein LGS lösen musst. Das ist insbesondere bei so kleinen Gleichungssystemen ganz nützlich. Bei sehr großen Gleichungssystemen ist das keine gute Idee, weil es da sehr auffällig aufwendig wird, die Inverse von \(J_f\) explizit zu berechnen; in etwa kubischer Aufwand zur Matrixgröße. Und da greift man lieber auf das lösen des LGS pro Iteration zurück.