% GENCOMB generate combinations
%    [X] = gencomb(n,d)
%       X : n by d^n matrix of combinations
%       n : number of components
%       d : number of states of each component
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

function [X] = gencomb(n,d)
X = ones(n,d^n);
for i = 1:d^n
   str = dec2base(i,d,n);
   for j = 1:n
      X(j,i) = str2num(str(j));
   end
end