clc;clear;
gc1=[0 2];
gs1=[1 0];
gc2=[2 2];
gs2=[0 1];
c1=[0 0 1 1];
c2=[0 1 0 1];
s1=[0 0 1 1];
s2=[0 1 0 1];
for m=1:4
    for n=1:4
        temp=c1(n).*gc1+c2(n).*gc2+s1(m).*gs1+s2(m).*gs2;
        x1(m,n)=temp(1);
        x2(m,n)=temp(2);       
    end
end
x1=mod(x1,4);
x2=mod(x2,4);