This manuscript is devoted to a derivative-free parametric iterative
step to obtain roots simultaneously for both nonlinear systems and
equations. We prove that when it is added to an arbitrary scheme, it
doubles the convergence order of the original procedure and defines a
new algorithm that obtains several solutions simultaneously. Different
numerical experiments are carried out to check the behaviour of the
iterative methods by changing the value of the parameter and the initial
guesses. Also, it is perform a numerical example where the dynamical
planes are carried out to see and compare the basins of attraction.