function P = customproc_markov(m,n) A = pnormdim(ones(m) - eye(m),1); C = ones(m,1)/m; P = markfillprob(A,C,n); return