function A = MakeStandard(k)
 
% Creates the standard matrix A satisfying A^2 = J
% where A is k^2-by-k^2.
 
for i=1:k
    A(1,i)=1;
end 
for i=2:k
    A(i, (k*(i-1)+1):k*i) = 1; 
end
for i=2:k 
    A((k*(i-1)+1):k*i, 1:k^2) = A(1:k,1:k^2);
end