如何用matlab画圆，如图所示的圆阵列？

如图所示，要求黑底白点，白色圆点的阵列，要求有代码。

hold on
flag=0:2*pi/360:2*pi
r=1
x=r*sin(flag)
y=r*cos(flag)
for i=0:5
for j=0:5
x1=x+i*4
y1=y+j*4
fill(x1,y1,[1 1 1])
end
end
set(gca,'color',[0 0 0])