function B=permut(A,p)
if p>0
  B=ones(A);
  N=2^p;
  Ipair=1:2:N;Iimp=2:2:N;
  B(Ipair)=permut(A((1:N/2)),p-1)
  B(Iimp)=permut(A((N/2+1):N),p-1)
else
  B=A
end
