wolffd@0: function M = softeye(K, p)
wolffd@0: % SOFTEYE Make a stochastic matrix with p on the diagonal, and the remaining mass distributed uniformly
wolffd@0: % M = softeye(K, p)
wolffd@0: %
wolffd@0: % M is a K x K matrix.
wolffd@0: 
wolffd@0: M = p*eye(K);
wolffd@0: q = 1-p;
wolffd@0: for i=1:K
wolffd@0:   M(i, [1:i-1  i+1:K]) = q/(K-1);
wolffd@0: end