In this paper, in order to solve systems of nonlinear equations, a new class of frozen Jacobian multi-step iterative methods is presented. Our proposed algorithms are characterized by a highly convergent order and an excellent efficiency index. The theoretical analysis is presented in detail. Finally, numerical experiments are presented for showing the performance of the proposed methods, when compared with known algorithms taken from the literature.

The local convergence analysis of the multi-step seventh order method to solve nonlinear equations is presented in this paper. The point of this paper is that our proposed study requires a weak hypothesis where the Fréchet derivative of the nonlinear operator satisfies the $\psi$-continuity condition, which thereby extends the applicability of the method when both Lipschitz and Hölder conditions fail. The convergence in this study is considered under the hypotheses on the first-order derivative without involving derivatives of the higher-order. To find a subset of the original convergence domain, a strategy is devised here. As a result, the new Lipschitz constants are at least as tight as the old ones, allowing for a more precise convergence analysis in the local convergence case. Some concrete numerical examples showing the performance of the method over some existing schemes are presented in this article.