function elementaryMatrix = emat(m, i, j, val)
  elementaryMatrix = eye(m,m);
  elementaryMatrix(i,j) = val;
end