clear;
di=input('Enter d_i=');
dq=input('Enter d_q=');
mi=length(di);
mq=length(dq);
m=mi+mq;
ni=2^mi;
nq=2^mq;
n=ni*nq;
pi=zeros(1,ni);
pq=zeros(1,nq);
for i=0:ni-1
   b=dec2bin(i,mi);
   for j=1:mi
      pi(i+1)=pi(i+1)+((b(j)-'0')*2-1)*di(j);
   end
   gi(i+1,:)=char(graycode(i,mi)+'0');
   gni(i+1,:)=graycode(i,mi);
end
for i=0:nq-1
   b=dec2bin(i,mq);
   for j=1:mq
      pq(i+1)=pq(i+1)+((b(j)-'0')*2-1)*dq(j);
   end
   gq(i+1,:)=char(graycode(i,mq)+'0');
   gnq(i+1,:)=graycode(i,mq);
end
k=1;
for i=1:ni
   for j=1:nq
      px(k)=pi(i);
      py(k)=pq(j);
      ii=1;
      iq=1;
      while (ii+iq <= m+1)
         if (ii<=mi)
            g(k,ii+iq-1)=gi(i,ii);
            gn(k,ii+iq-1)=gni(i,ii);
            ii=ii+1;
         end
         if (iq<=mq)
            g(k,ii+iq-1)=gq(j,iq);
            gn(k,ii+iq-1)=gnq(j,iq);
            iq=iq+1;
         end
      end
      k=k+1;
   end
end

plot(px,py,'o');
for i=1:n
   text(px(i),py(i)-0.05,g(i,:));
end