function c=source(len)
a=randint(1,len*2);
b=reshape(a,2,len);
c=zeros(1,len);
for j=1:len
    if b(1,j)==0
        if b(2,j)==0
            c(j)=0;
        else c(j)=1;
        end
    else if b(2,j)==0
            c(j)=3;
        else c(j)=2;
        end
    end
end