!! OBSOLETE !!
Its functionality is now implemented in my function PERMN, available on the File Exchange
permnsub - subset of all permutations with repetition
M = permnsub(V,N,IX) returns a subset of all permutations, with
repetition, of N elements, taken from the elements in vector V. M has
the size numel(IX)-by-N. V can be an array of numbers, cells or strings,
and M will be of the same type as V.
This is an excellent implementation of a very efficient method. Together with the exhaustive documentation this submission is useful, usable and valuable for the productive work and education. Thanks, Jos!
A tiny simplification: You do not need B(:)', because B is a row vector already.
24 Jun 2013
1.1 - implemented suggestion of Jan; output is now always a row vector