[EN] In this work we introduce a new operator of divided differences that preserves the convergence order when it is used for approximating the Jacobian matrix in iterative method for solving nonlinear systems. We obtain ...