evl:= V -> convert(evalm(V), list);
`&.`:= (U,V) -> linalg[dotprod](U, V, orthogonal);
len:= V -> sqrt(V &. V);
unit:= V -> evl(V/sqrt(V &. V));
`&x`:= (U,V) -> convert(linalg[crossprod](U,V),list);