% Triangle équilatéral
a=1;
x=[0 a a/2 0];
y=[0 0 sqrt(3)*a/2 0];
clf
hold on

n=14;
m=10;
for i=1:n
	for j=1:m
		h=plot(x+i+rem(j,2)*a/2,y+j*sqrt(3)*a/2);
		set(h,'color',[0.5 0.5 0.5],'linewidth',3)
	end
end
for i=1:n
	for j=1:m
		plot(i+a/2+rem(j,2)*a/2,j*sqrt(3)*a/2,'.w','markersize',18)
		plot(i+a/2+rem(j,2)*a/2,j*sqrt(3)*a/2,'.k','markersize',16)
		plot(i+a/2+rem(j,2)*a/2,(j+1)*sqrt(3)*a/2,'.w','markersize',28)
		plot(i+a/2+rem(j,2)*a/2,(j+1)*sqrt(3)*a/2,'.k','markersize',20)
	end 
end
set(gcf,'color','w')

axis equal
axis off
h=plot([1.75 n+.75 n+.75 1.75 1.75],[1.5 1.5 m-.5 m-.5 1.5]);
set(h,'color',[0.5 0.5 0.5],'linewidth',3)
axis([1.75 n+.75 1.5 m-0.5])
hold off
