The Jacobian matrix has the following form. It is updated either for each iteration or other modified schemes may be used to achieve convergence with less computational effort.
