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