Reference: isNormalPanmagicSquare
isNormalPanmagicSquare(M), where M is a real matrix, returns true if and only if M is a normal panmagic (or, diabolic) square, i.e. a N×N panmagic square in which every integer from 1 to N^2 occurs (exactly once). Thus, isNormalPanmagicSquare(M) = isPanmagicSquare(M) ∧ isNormalSquare(M).