>	a1:=sin(delta):a:=cos(delta):Z:=matrix(3,3,[a,-a1,0,a1,a,0,0,0,1]);

>	b1:=sin(beta):b:=cos(beta):Y:=matrix(3,3,[b,0,b1,0,1,0,-b1,0,b]);

>	c1:=sin(alpha):c:=cos(alpha):X:=matrix(3,3,[1,0,0,0,c,-c1,0,c1,c]);

>	XYZ=multiply(X,Y,Z);

>	XZY=multiply(X,Z,Y);

>	YXZ=multiply(Y,X,Z);

>	YZX=multiply(Y,Z,X);

>	ZYX=multiply(Z,Y,X);

>	ZXY=multiply(Z,X,Y);

>