In this scheme, the tangential stiffness matrix is formed and decomposed at each iteration within a particular step as shown in the figure below. The NR method has a high convergence rate and its rate of convergence is quadratic. However, since the tangential stiffness is formed and decomposed at each iteration, which can be prohibitively expensive for large models, it may be advantageous to use another iterative method.