function D = pairwisedist(m,n)
dims = repmat([m],[1 n]);
N = m^n;
D = zeros(N);
for ii = 1:N
  x = ind2coord(ii,dims);
  for jj = 1:N
    y = ind2coord(jj,dims);
    D(ii,jj) = length(find(x~=y));
  end
end
return