dot := (u,v)-> sum(u[i]*v[i],i=1..3); alias(cp=crossprod): u := vector(3): v:=vector(3):w:=vector(3):z:=vector(3): lh := cp( cp(u,v), cp(w,z) ): rh := dot( u, cp(w,z))*v - dot(v, cp(w,z))*u: NewF := lh = rh: simplify(evalm(lh-rh)[1]); simplify(evalm(lh-rh)[2]); simplify(evalm(lh-rh)[3]);