% die_mar alternates between fair and loaded die rolls

A  = [[.95 .05]; [.1 .9]];     % transition probabilities among hidden states
E  = [[1 1 1 1 1 1]/6; [1 1 1 1 1 5]/10]; % emission probabilities

Pi = [];
x  = [];

Pi(1) = rando([1 1]/2);
x(1)  = rando(E(Pi(1),:));     % use row Pi(1) of E to get emission probs

for i=1:74,
  Pi(i+1) = rando(A(Pi(i),:));
  x(i+1)  = rando(E(Pi(i+1),:));
end
show(Pi,'FL');
show(x,'123456');