   function   X=usMat(x,n)
%  function   X=usMat(x,n);
%  Use the isometry version
global  indsu ;

X=zeros(n);
X(indsu)=x/sqrt(2);  
X=X+X';
