#	How to do a cross-product (4x4) with 3 products

#	Form the products
for i to 3 do
	p.i := (x.i.1*a1+x.i.2*a2+x.i.3*a3+x.i.4*a4) *
		(y.i.1*b1+y.i.2*b2+y.i.3*b3+y.i.4*b4)
	od;

#	Form the linear combinations of products
eq := - (a1*b1+a2*b2+a3*b3+a4*b4);
for j to 3 do eq := eq + p.j od;

eqa := {};
eq := expand(eq);
for i to 4 do
	for j to 4 do
		eqa := eqa + {coeff( coeff(eq,a.i,1), b.j, 1 )} od od;

printlevel := 10000;
read `/u/ghgonnet/handbook/6/1/solsys3`;
gc();
prettyprint := 0;
