Abstract
In this work, we propose a higher-order derivative free family for solving nonlinear systems. The local order of convergence of the constructed family is first determined using first-order divided difference operators for functions of several variables and direct computation by Taylor's series expansion. Computational efficiencies of the developed scheme is considered and compared with their closest competitors. Moreover, numerical experiments are performed on several large and complex non-linear systems. The results are found to be effective and comparable to other existing methods which also support the theoretical results.