r=1; %圆半径
n=6; %n边形,n>=3
%画圆
tt = (0:0.01:1)*2*pi;
x=r*sin(tt);
y=r*cos(tt);
plot(x,y);
axis equal
hold on;
%内接圆
theta = (0:1/n:1)*2*pi;
[x,y]=pol2cart(theta,r);
plot(x,y);
%外接圆
theta = pi/n*ones(1,n+1)+(0:1/n:1)*2*pi;
r1 = r/cos(2*pi/2/n)*ones(1,n+1);
[x,y]=pol2cart(theta,r1);
plot(x,y);