answered 8 years ago
This is sufficient for vectors:
V = matrix(4,1,[1,2,2,1]) V.transpose() / V.norm()^2