THETA = 0.45;
DIM = 60;
NUM_ITERATIONS = 1000;
SHOW_EVERY = 100;

x = -3 + 2 * randi(2, [DIM, DIM]);

figure;
for it = 1:NUM_ITERATIONS
  x = comb_gibbs_step(x, THETA);
  if mod(it, SHOW_EVERY) == 0
    subplot(4, 3, it/SHOW_EVERY);
    imagesc(x);
    colormap gray;
    axis off;
  end
end
