def identity( n ) = vector( n, n, \r, c -> if r == c then 1 else 0 )

println( identity(3) )