This function returns a list of n mutually orthogonal Latin squares of size p^m.

Two Latin squares, A and B, are mutually orthogonal if the set of ordered pairs $\left\{A_{i,j}\,B_{i,j}\,0\le i\,j\le p^m\right\}$ has no repeated elements.

$\mathrm{MOLS}\left(3\,1\,2\right)$

$\left[\left[\begin{array}{ccc}0& 1& 2\\ 1& 2& 0\\ 2& 0& 1\end{array}\right]\,\left[\begin{array}{ccc}0& 1& 2\\ 2& 0& 1\\ 1& 2& 0\end{array}\right]\right]$