a := vector([a1,a2,a3]); b:= vector([b1,b2,b3]); c := vector([c1,c2,c3]); alias(cp=crossprod); axb := cp(a,b); bxc:=cp(b,c); a_xbxc := cp(a,bxc); axbx_c:=cp(axb,c); evalm(a_xbxc - axbx_c); map(simplify,"); res := ": solutions := ": solutions[16]; u :=[-2,0,-1]; v:=[x,y,z]; w:=[2,0,1]; u := [1,0,-1]; cp(u,cp(v,w)); cp(cp(u,v),w);