As a starting point...
N=6;
t=2*pi*(1/(2*N):1/N:1).';
x=cos(t);y=sin(t);
fill(x,y,'k')
axis square
hold on
hH=fill(x+2*max(x),y,'r');
axis square
ylim(xlim)
Those on the same row are obvious, the fillers in between need an adjustment of the x positon of only 1*x and of the height of the vertical side for y for the initial point -- then it's also tiled by 2X the X,Y position.
<fill> has a builtin replication routine that, I think, can be used to build the full array automagically, but I've not dug into just how to use it/how it works; there's not an example and didn't have time(*) to dig at the moment.
(*) Well, more like "inclination" would be more accurate... :)
