function P = treefillprob(m,n,T,psi)
global Ztree
P = zeros(m^n,1);
dims = repmat([m],[1 n]);
[nodesi,nodesj]=find(T);
ii = find(nodesi<nodesj);
nodesi = nodesi(ii);
nodesj = nodesj(ii);
ne = length(nodesi);
if size(psi,3) < ne
  psi = repmat(psi,[1 1 ne]);
end
for xind = 1:m^n
  x = ind2coord(xind,dims);
  p = 1;
  for e=1:length(nodesi)
    i = nodesi(e);
    j = nodesj(e);
    p = p * psi(x(i),x(j),e);
  end
  P(xind) = p;
end
Ztree = sum(P);
P = P/Ztree;
return