function P = markovfillprob(A0,C0,n0)
global A C n m
A=A0; C=C0; n=n0;
m = size(A,1);
if size(A,3)==1
  A = repmat(A,[1 1 n-1]);
end

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

P = zeros(m^n,1);
for sind = 1:m^n
  s = ind2coord(sind,dimS);
  P(sind) = Ps(s);
end

return

function p = Ps(s)
% p = P(s,x)
global A C n
p = C(s(1));
for i=2:length(s)
  p = p * A(s(i),s(i-1),i-1);
end
return