Que faites vous?

let magique v=

let N=ref (vect_length(v.(0)))
and k1=ref 0 and S1=ref 0 and S2=ref 0
and test=ref true and k3=ref 0
and V=(make_vect ((vect_length(v.(0)))*(vect_length(v.(0)))) 0) in

while !k1 < !N & !test do
for k2=0 to !N-1 do
V.(!k3)<-v.(!k1).(k2);
S1:= !S1 + v.(!k1).(k2);
S2:= !S2 + v.(k2).(!k1);
for k4=0 to !k3-1 do
if V.(!k3)=V.(k4) then test:=false;
done;
if V.(!k3)<1 or V.(!k3)>(!N)*(!N) then test:=false;
k3:=!k3+1;
done ;
if !S1<>(!N*(!N*(!N)+1)/2) or !S2<>(!N*(!N*(!N)+1)/2) then test:=false ;
k1:=(!k1)+1;S1:=0; S2:=0;
done ;

for k2=0 to !N-1 do
S1:= !S1 + v.(k2).(k2);
S2:= !S2 + v.(!N-k2-1).(!N-k2-1);
done;

if !S1<>(!N*(!N*(!N)+1)/2) or !S2<>(!N*(!N*(!N)+1)/2) then test:=false ;
!test;;

magique [|[|4;9;2|];[|3;5;7|];[|8;1;6|]|];;

let create n=
let V=(make_matrix n n 0) and x=ref (n/2+1) and y=ref (n/2) in
V.(!x).(!y)<-1;

for k=2 to n*n do
if V.( (!x+1) mod n ).( (!y+1) mod n ) = 0 then
begin
x:=(!x+1) mod n;
y:=(!y+1) mod n;
V.(!x).(!y)<-k;
end
else begin
x:=(!x+2) mod n;
y:=(!y) mod n;
V.(!x).(!y)<-k;
end;done;

V;;

create 5;;