samer@16: % nary - return array contain nary sequece over N columns
samer@4: %
samer@16: % nary :: M:natural, N:natural -> [[M^N,N]->1..M].
samer@4: 
samer@4: function B=nary(M,N)
samer@4: 
samer@4: if (N==1), B=(1:M)';
samer@4: else
samer@4: 	b=nary(M,N-1);
samer@4: 	m=size(b,1);
samer@4: 	B=[ kron((1:M)',ones(m,1)), repmat(b,M,1)];
samer@4: end