rc=4;dy=2*rc;dx=rc*sqrt(3);
A=pi/3*[1:7];
hold on;
for yk=[0:dy:100,0:-dy:-100];
yfun=inline(['sqrt(3)*x/3+',num2str(yk)]);
for xk=[0:dx:100,0:-dx:-100];
xp=xk;
yp=yfun(xp);
if -50<xp & xp<50 & -50<yp & yp<50;
plot([xp+i*yp]+rc*exp(i*A)*2/sqrt(3),'k','linewidth',2);
plot(xp,yp,'.','markersize',6);
end
end
end
axis square
plot([0,0],[-50,50],'r:');
plot([-50,50],[0,0],'r:');
plot(50*sqrt(2)*exp(i*[pi/2*[1:5]+pi/4]),'b');