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 + Dx, y := y + Dy, r := r + Dr.