function A = aa(n)
% aa(n) creates an n-by-n matrix
for i = 1:n
  for j = 1:n
    if(i >= j)
      A(i,j) = 1;
    else
      A(i,j) = 0;
    end
  end
end
