function P=permhaar(n0)
global m n
%n = n0-1;
n = n0;
m = n0;

P = ones(m^n,1);
for xi=1:m^n
  x = i2c(xi);
  P(xi) = isperm(x);
end

P = P+eps;
P=P/sum(P);
return


function b = isperm(x)
Ux = getucnt(x,unique(x));
b = (length(x)==length(Ux));
return


function Ucnt = getucnt(U,Uniq)
Ucnt = zeros(size(Uniq));
for j=1:length(Uniq)
  Ucnt(j) = length(find(U==Uniq(j)));
end
return