clear;theta=input('Enter Theta=');d=2*(pi/360)*theta;m=length(d);n=2^m; p=zeros(1,n);for i=0:n-1   b=dec2bin(i,m);   for j=1:m      p(i+1)=p(i+1)+((b(j)-'0')*2-1)*d(j);   end   g(i+1,:)=char(graycode(i,m)+'0');endpolar(p,ones(1,n),'o');for i=1:n   text(cos(p(i)),sin(p(i)),g(i,:));end