p=eye(21,21);
p1=eye(21,21);
pm=[];
c=zeros(252,672);
a1=1:21;a2=22:42;a3=43:63;a4=64:84;a5=85:105;a6=106:126;a7=127:147;a8=148:168;a9=169:189;a10=190:210;a11=211:231;a12=232:252;a13=253:273;a14=274:294;a15=295:315;a16=316:336;a17=337:357;a18=358:378;
a19=379:399;a20=400:420;a21=421:441;a22=442:462;a23=463:483;a24=484:504;a25=505:525;a26=526:546;a27=547:567;a28=568:588;a29=589:609;a30=610:630;a31=631:651;a32=652:672;
for i=1:20
p=[p(:,2:21) p(:,1)];
pm{i}={p};
end
c(1:21,a1)=p1;
c(1:21,a4)=pm{5}{1};
c(1:21,a6)=pm{18}{1};
c(1:21,a7)=pm{16}{1};
c(1:21,a11)=pm{3}{1};
c(1:21,a12)=pm{6}{1};
c(1:21,a13)=pm{10}{1};
c(1:21,a16)&