Constructing a High-Order Globally Convergent Iterative Method for Calculating the Matrix Sign Function

Figure 2

Basins of attractions for (17) in left and its density plot on the same domain in right (shading is done based on the number of cycles to achieve the convergence).