wolffd@0: function bnet = mk_markov_chain_bnet(N, Q)
wolffd@0: 
wolffd@0: dag = zeros(N);
wolffd@0: dag(1,2)=1; dag(2,3)=1;
wolffd@0: ns = Q*ones(1,N); 
wolffd@0: bnet = mk_bnet(dag, ns);
wolffd@0: for i=1:N
wolffd@0:   bnet.CPD{i} = tabular_CPD(bnet, i);
wolffd@0: end