Function rquat() returns a quaternionic vector,
roct() returns an octonionic vector, and romat() a
quaternionic matrix.
Functions rquat() and roct() give a quick “get
you going” random onion to play with. Function romat() gives
a simple onionmat, although arguably matrix(roct(4),2,2) is
as convenient.