function Pii = getmargX(P,m,n,ii)
% marginalizes out over j\notin ii
% yields correct results for ii=1:n and ii=[]

lii = length(ii);
Pii = zeros(m^lii,1);
dimi = repmat([m],[1 lii]);
for xi = 1:m^lii
  X = ind2coord(xi,dimi);
  p = getPX(P,m,n,X,ii);
  Pii(xi) = p;
end
return