% 畫蜂房（六角形網格）
n=6;a=0:2*pi/n:2*pi;y0=sin(a);x0=cos(a);
p=5;q=5;x1=[];y1=[];x2=[];y2=[];
for i=0:p
    for j=0:q
        x=x0+i*1.5;y=y0+(j+mod(i,2)/2)*sqrt(3);
        plot(x,y),hold on
    end
end
%x=[x1,x2];y=[y1,y2];plot(x,y),axis equal
axis equal