If \(trt = p^n\) is a prime power, then \(trt\)-1 latin squares of order \(trt\)
are constructed.
The elements of the squares are numbered 1,…,\(trt\).
These squares are mutually orthogonal, i.e. if any two of them are superimposed, the resulting
array will contain each ordered pair \((i,j)\), \(i\),\(j\) in {1,…,
\(trt\)} exactly once.
The squares are in standard order, i.e. the first row is always equal to (1,…,\(trt\)).
A primitive polynomial may be constructed automatically using the internal function GF
.