In the note we consider an iterative generalisation of the rejection sampling method. In high energy physics, this sampling is frequently used for event generation, i.e. preparation of phase space points distributed according to a matrix element squared |M |2 for a scattering process. In many realistic cases |M |2 is a complicated multi-dimensional function, so, the standard von Neumann procedu...