function m=m16(n)
q=zeros(16,1);
q(1)=n;
for i=1:15
 q(i+1)=floor(q(i)/2);
end 
m=zeros(4,4);
for i=0:3
 for j=1:4
  m(i+1,j)=(-1)^q(4*i+j);
 end
end
end
