function P = markfillprob(A,C,n)
m = size(A,1);

dims= repmat(m,[1 n]);

P = zeros(m^n,1);

for xind = 1:m^(n)
  x = ind2coord(xind,dims);
  P(xind) = P(xind) + Pmark(x,A,C);
end

return

function p = Pmark(s,A,C)
p = C(s(1));
for i=2:length(s)
  p = p * A(s(i),s(i-1));
end
return