clear all
close all
G=[0 1 0 1 1 0 0; 1 0 1 0 1 0 0; 0 1 1 0 0 1 0; 1 1 0 0 0 0 1];
b=0:15;
a=dec2bin(b,4);
for i=1:length(b)
    u=a(i,:);
    c(i,:)=zeros(1,length(G));
    for j=1:length(u)
        h(j,:)=G(j,:)*u(j);
        c(i,:)=mod(c(i,:)+h(j,:),2);
    end
end
c