%
%  macm 202 -- 31 mar 03 -- djm
%
%  w12shuff8.m:  

Ncards = 52;  Nh = Ncards/2;

Pn = zeros(Ncards,Ncards);

for k=1:Nh
    Pn(2*k-1,   k) = 1;
    Pn(2*k  ,Nh+k) = 1;
end
Pn = sparse(Pn);

figure(1); clf

k=0;
shuff = eye(Ncards);
spy(shuff);
title(['\bf after shuffles # ' num2str(k)])

for k=1:8
    pause(1)
    shuff = Pn*shuff;
    spy(shuff)
    title(['\bf after shuffles # ' num2str(k)])
end