Creates all fundamentally different panmagic squares of order 4.
Usage
panmagic.4(vals = 2^(0:3))
Arguments
vals
a length four vector giving the values which are combined
in each of the $2^4$ possible ways. Thus
vals=2^sample(0:3) always gives a normal square
(0-15 in binary).