function A1 = mcstate_blowup(A,k)
% turn A(x|y) into A1(x^k|y) where x^k \in Omega^k

m = size(A,1);
A1 = zeros(m^k,m);

for y = 1:m
  for xi=1:m^k
    x = i2c(xi,m,k);
    p = A(x(1),y);
    for i=2:k
      p = p * A(x(i),x(i-1));
    end
    A1(xi,y) = p;
  end
end

return


function i=c2i(c,m,n)
dims = repmat(m,[1 n]);
i = coord2ind(c,dims);
return

function c=i2c(i,m,n)
dims = repmat(m,[1 n]);
c = ind2coord(i,dims);
return