function B = bbmat(n)
% bbmat(n) creates an n-by-n matrix.
for i = 1:n
  for j = 1:n
    B(i,j) = max(1,3*(min(i,j)-1));
  end
end
