// vérification de magic 
A=testmatrix('magic',5)
v=[sum(A,'r'),sum(A,'c')'];
and(v == v(1))