ord=[3 4 6 2^10]for i=1:4 subplot(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