Function FAURE generates elements 0,..,k of a d-dimensional, base-b Faure sequence, an ingredient of some 'quasi Monte Carlo' (QMC) algorithms. See the function's help header for details.
Disclaimer: the implementation is inefficient (forgoing the opportunity to vectorize some calculations, or exploit the recursive nature of generator matrices, for example), solely for the sake of greater transparency.