r = 5; % Radius of small circle
rho = 50; % Radius of big circle
Np = 16; % Number of circle
theta = linspace(pi,-pi,Np+1);
[x,y] = pol2cart(theta,rho); % coordinates of small circle center
t = linspace(0, 2*pi);
apoint = 1:Np; % Line connection, starting point
bpoint = circshift(1:16,-4); % Line coonection, end point
for w = 1:Np
line([x(apoint(w)) x(bpoint(w))],[y(apoint(w)) y(bpoint(w))],'Color','k'); % Draw lines first
end
rectangle('Position',[-rho -rho 2*rho 2*rho],'Curvature',[1 1])% % Draw big circle
for k = 1:Np
patch(x(k)+r*cos(t), y(k)+r*sin(t), 'w','FaceAlpha',1); % Draw small circle
end
text(x(1:Np),y(1:Np),string(1:Np),'HorizontalAlignment','center'); % Add the text
axis equal
axis off
