Returns at least 3 orthogonal latin square arrays of dimension v. The first and second arrays are the
rows and columns of a v x v square and the third is a Latin square of size v x v.
If v is prime or prime power in the set 4, 8, 16, 32, 64, 128, 9, 27, 81, 25, 49 there
are v-1 MOLS and orthogLS(v)
will return a total of v+1 arrays. If v = 10 there are two MOLS
and orthogLS(10)
will give a total of 4 arrays ofsize 10 x 10.
NB If S1 and T1 are mutually orthogonal Latin squares of order n1 and
S2 and T2 are mutually orthogonal Latin squares of order n2 then the product squares S1xS2
and T1xT2 are orthogonal to each other and have order n1n2 (not yet implemented)