function hexa=pavagehexa(n,m)
% pavage hexagonal
[x,y]=hexagone;
r1=sqrt(3)/2;
r2=3/2;
cla
hold on
for i=1:n
   for j=1:m
      if rem(j,2)
         X=x+2*i*r1;
      else
         X=x+(2*i+1)*r1;
      end
      hexa(i,j)=fill(X,y+j*r2,'c','buttondownfcn','clic','userdata',0);
   end
end
axis equal
axis off

