Abstract
A general class of multi-step iterative methods for finding approximate real or complex solutions of nonlinear systems is presented. The well-known technique of undetermined coefficients is used to construct the first method of the class while the higher order schemes will be attained by a frozen Jacobian. The point of attraction theory will be taken into account to prove the convergence behavior of the main proposed iterative method. Then, it will be observed that an m-step method converges with 2m-order. A discussion of the computational efficiency index alongside numerical comparisons with the existing methods will be given. Finally, we illustrate the application of the new schemes in solving nonlinear partial differential equations.