wolffd@0: function transmat = mk_leftright_transmat(Q, p)
wolffd@0: % MK_LEFTRIGHT_TRANSMAT Q = num states, p = prob on (i,i), 1-p on (i,i+1)
wolffd@0: % function transmat = mk_leftright_transmat(Q, p)
wolffd@0: 
wolffd@0: transmat = p*diag(ones(Q,1)) + (1-p)*diag(ones(Q-1,1),1);
wolffd@0: transmat(Q,Q)=1;