In this paper, we deal with a new family of iterative methods for approximating the solution of nonlinear systems for non-differentiable operators. The novelty of this family is that it is a m-step generalization of the Steffensen-type method by updating the divided difference operator in the first two steps but not in the following ones. This procedure allows us to increase both the order of convergence and the efficiency index with respect to that obtained in the family that updates divided differences only in the first step. We perform a semilocal convergence study that allows us to fix the convergence domain and uniqueness for real applied problems, where the existence of a solution is not known a priori. After this study, some numerical tests are developed to apply the semilocal convergence theoretical results obtained. Finally, mediating the dynamic planes generated by the different numerical methods that compose the family, we study the symmetry of the basins of attraction generated by each solution, the shape of these basins, and the convergence to each root of a polynomial function.