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