If \(v_1,\ldots,v_n\) is the standard basis for
\(\mathbb{R}^n\) then \(\phi_i\) is defined so
that \(\phi_i(v_j)=\delta_{ij}\). phi(n)
returns
\(\phi_n\).
If n
is a vector of strictly positive integers, then
phi(n)
returns the tensor cross product of \(\phi\)
applied to the individual elements of n
[which is a lot
easier and more obvious than it sounds].