# Returns an identity matrix.

ident = function( n )
{
    return diag( fill( scalar (n); 1.0 ) );
};