a = [1 1 2 3]';
  b = [2 3 1 3]';
  n = 4;

  a_mat = repmat(a, [1 n]);
  b_mat = repmat(b', [n 1]);

  A = a_mat == b_mat;