function E=elemat(n,i,j,val);
% E=elemat(n,i,j,val) -- compute the elementary matrix E
%      of dimension n by n with entry val in (i,j) position

E = eye(n); E(i,j) = val;

%END


