[toy_network, toy_factors] = ConstructToyNetwork(1, 0.1);

for i =1:10
  GibbsTrans(repmat(1, 1, length(toy_network.card)), toy_network, toy_factors)
endfor

disp("****************************************************************")

for i = 1:10
  GibbsTrans(repmat(2, 1, length(toy_network.card)), toy_network, toy_factors)
endfor
