I think it’s analogous to the identity operator I in LinearAlgebra — when you act I*x on a vector (or matrix), it infers the correct dimension of the identity from x, analogous to UnsetSpace in ApproxFun. But occasionally it’s nice to have an explicit identity matrix, e.g. via I(n) to get an n \times n Diagonal matrix.
3 Likes