function six
close all
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');