Stability Anomalies of Some Jacobian-Free Iterative Methods of High Order of Convergence

In this manuscript, we design two classes of parametric iterative schemes to solve nonlinear problems that do not need to evaluate Jacobian matrices and need to solve three linear systems per iteration with the same divided difference operator as the coefficient matrix. The stability performance of the classes is analyzed on a quadratic polynomial system, and it is shown that for many values of the parameter, only convergence to the roots of the problem exists. Finally, we check the performance of these methods on some test problems to confirm the theoretical results.