function r=randomFlipMatrix(n);
% returns matrix with components 1 at n random positions of each column 
r=zeros(156,26);
for i=1:26
    x=randperm(156);
    r(x(1:n),i)=1;
end
