MAGIC matrix function
Description
MAGIC(n) creates a n × n matrix with as elements the integers 1 through n2, with equal row and column sums.
The argument n must be larger than 2.
If n=1 then the function returns the scalar {1}; if n=2 then a non-magic square matrix is returned.