The Jacobian of
F
(
x
,
y
,
r
) is shown here.
Solving
F
by Newton's method requires at each iteration step solving the linear system and updating values as
x
:=
x
+
D
x
,
y
:=
y
+
D
y
,
r
:=
r
+
D
r
.