>function IdentityMatrix (n:index)
$  return setdiag(zeros(n,n),0,1);
$endfunction