functions in matrix.i - u

             unit(n)  
          or unit(n, m)  
 
     returns n-by-n (or n-by-m) unit matrix, i.e.- matrix with diagonal  
     elements all 1.0, off diagonal elements 0.0  

interpreted function, defined at i0/matrix.i   line 21