function P = latentfillprob(Ph,Poh,mhid,mobs,n)


dimX = repmat([mobs],[1 n]);
dimS = repmat([mhid],[1 n]);

P = zeros(mobs^n,1);

for xind = 1:mobs^(n)
  x = ind2coord(xind,dimX);
  for sind = 1:mhid^n
    s = ind2coord(sind,dimS);
    P(xind) = P(xind) + Ph(sind)*Pcond(Poh,x,s);
  end
end

return

function p = Pcond(Poh,x,s)
% p = P(x|s)
p = 1;
for i=1:length(x)
  p = p * Poh(x(i),s(i),i);
end
return