(* Définition de la matrice de connectivité ; code R. Mosseri *)

Trouveppl[l_]:=Min[(Reap[Do[If[l[[j]][[1]]==0,Sow[l[[j]][[2]]]],{j,1,Length[l]}]])[[2]]];

MatConFluxPBC[v_,p_,npt_,f_,th_]:=Module[{pnorm,l,res},l=Length[p];pnorm=qqchose*p;
res=Table[0,{x,l},{y,l}];Do[res[[j,jj]]=-Exp[I*Pi*f*(-Sin[th]*(p[[jj,1]]-p[[j,1]])+Cos[th]*(p[[jj,2]]-p[[j,2]]))*(Cos[th]*(p[[jj,1]]+p[[j,1]])+Sin[th](p[[jj,2]]+p[[j,2]]))];res[[jj,j]]=-Exp[I*Pi*f*(-Sin[th]*(p[[j,1]]-p[[jj,1]])+Cos[th]*(p[[j,2]]-p[[jj,2]]))*(Cos[th]*(p[[jj,1]]+p[[j,1]])+Sin[th](p[[jj,2]]+p[[j,2]]))],{jj,1,Length[v[[j]]]},{j,1,l}];res];

(*Construction pour k=4 ; npt = 28, taille = 28*Sqrt[2] *)

Trouveppl[{{0,1},{0,2}}]