In this paper we present a new method for determining simultaneously all the simple roots of a quaternionic polynomial. The proposed algorithm is a two-step iterative Weierstrass-like method and has cubic order of convergence. We also illustrate a variation of the method which combines the new scheme with a recently proposed deflation procedure for the case of polynomials with spherical roots.