Path: csiph.com!v102.xanadu-bbs.net!xanadu-bbs.net!feeder.erje.net!eu.feeder.erje.net!eternal-september.org!feeder.eternal-september.org!news.eternal-september.org!.POSTED!not-for-mail
From: Joe Riel <joer@san.rr.com>
Newsgroups: comp.soft-sys.math.maple
Subject: Re: Solving Equation Set
Date: Mon, 09 Jun 2014 11:17:46 -0700
Organization: A noiseless patient Spider
Lines: 47
Message-ID: <87ppiiymk5.fsf@san.rr.com>
References: <QrmdnbmNiOaMtQjOnZ2dnUVZ_rCdnZ2d@megapath.net>
Mime-Version: 1.0
Content-Type: text/plain; charset=us-ascii
Injection-Info: mx05.eternal-september.org; posting-host="15591ad2607da309a0d1a78a1d632bc7"; logging-data="31966"; mail-complaints-to="abuse@eternal-september.org";	posting-account="U2FsdGVkX1+J1KZ9gZSaJGh48jZva2qY"
User-Agent: Gnus/5.13 (Gnus v5.13) Emacs/23.4 (gnu/linux)
Cancel-Lock: sha1:AXdrfI6njdh+wlzeG//NC6LCbPU= sha1:E3VFb92BjpvjtZS/kdq9+pQwBDg=
Xref: csiph.com comp.soft-sys.math.maple:889

"Thomas D. Dean" <tomdean@speakeasy.org> writes:

> From an email about another CAS, Maxima
>
> ##
> ## In the following sequence, I'm trying to prove that the length of the
> ## cross product of two orthogonal unit vectors is =1.
> ##
> ## The normal of a vector is the length of the vector
> ## |A| = sqrt(a1^2+a2^2+a3^2)
> ## if |A| = 1 then
> ## |A|^2 = sqrt(a1^2+a2^2+a3^2)^2 = 1^2
> ## so that
> ## a1^2+a2^2+a3^2 = 1
> ##
> ## N(A)=|A|^2=a1^2+a2^2+a3^2 see above
> ## N(B)=|B|^2=b1^2+b2^2+b3^2 see above
> ## N(AxB) is the norm of AxB.
> ## A.B=a1*b1+a2*b2+a3*b3=0.
> ##
> restart;
> A:=<a1,a2,a3>;
> B:=<b1,b2,b3>;
> tmp:=LinearAlgebra[CrossProduct](A,B);
> eq1:=tmp . tmp = NAxB;  ## NAxB is the cross product
> eq2:=A . A = 1;         ## A is a unit vector
> eq3:=B . B = 1;         ## B is a unit vector
> eq4:=A . B = 0;         ## A is perpendicular ot B
>
> solve({eq1,eq2,eq3,eq4},NAxB);
>
> sol1:=lhs(eq1)+lhs(eq4)^2 = rhs(eq1)+rhs(eq4);
> sol2:=eliminate([sol1,eq2,eq3],b3);
> sol3:=eliminate(sol2[2][2],a3);
> sol4:=solve(sol3[2][2],NAxB);
>
> I can peek and poke to get a solution.  There should be an easier way.
> How do I do this?
>
> Tom Dean

simplify(tmp.tmp, [A.A=1,B.B=1,A.B=0]) assuming real;

            1

-- 
Joe Riel