function[B]=bigmat(b,r)
for i=1:r
B(:,i)=b;
end
