matthiasm@8: function M = sampleUniformInts(N, r, c)
matthiasm@8: 
matthiasm@8: % M is an rxc matrix of integers in 1..N
matthiasm@8: 
matthiasm@8: prob = normalize(ones(N,1));
matthiasm@8: M = sample_discrete(prob, r, c);