ord=[3 4 6 2^10]for i=1:4subplot(2,2,i)theta=linspace(pi/ord(i),2*pi+pi/ord(i),ord(i)+1);plot(cos(theta),sin(theta));xlim(1.5*[-1,1]);ylim(1.5*[-1,1]);axis equal;end