MathGroup Archive 1991

[Date Index] [Thread Index] [Author Index]

Search the Archive

Re: Question about Simplify and related

  • To: mathgroup at
  • Subject: Re: Question about Simplify and related
  • From: mbp at (Mark B. Phillips)
  • Date: Sat, 19 Jan 91 18:07:47 EST

> 	This is a question about Simplify and related objects. Below is 
> a Mma session in which a 3x3 orthogonal matrix (u) is formed.
> ...
> In fact I have gotten to the point below where one can show in a
> couple of minutes with a paper and pencil that uutrans is indeed the
> 3x3 unit matrix.
> 	Sould all this be so hard to do, or am I missing something 
> simple and obvious?

Part of the problem lies in that your simplification rules like si^2 +
ci^2 -> 1 really need to to be applied at different times to different

For example to simplify:		c3^2 + c1^2 s3^2 + s1^2 s3^2
using only the rules
  si^2 + ci^2 -> 1
one needs to first factor
s3^2 out of the 2nd two terms:		c3^2 + s3^2( c1^2 + s1^2 )
apply s1^2 + c1^2 -> 1:			c3^2 + s3^2
and then apply c3^2 + s3^2 -> 1:	1.

Fortunately in this case the rules are simple enough to be written in
a form which will completely eliminate a variable when applied.  Using
si^2 -> 1 - ci^2 works much better because it complelely eliminates
si^2 from the expression:

In[1]:= u={{c1 c2 - c3 s1 s2,c2 s1 + c3 c1 s2,s2 s3},
      {-s2 c1 - c3 s1 c2, -s2 s1 + c3 c1 c2, c2 s3},
       {s3 s1, -s3 c1, c3}}
Out[1]= {{c1 c2 - c3 s1 s2, c2 s1 + c1 c3 s2, s2 s3}, 
>    {-(c2 c3 s1) - c1 s2, c1 c2 c3 - s1 s2, c2 s3}, {s1 s3, -(c1 s3), c3}}

In[2]:= uutrans = Expand[ u . Transpose[u] ]

Out[2]= [[ lengthy output deleted ]]

In[3]:= Expand[ uutrans /. {s1^2->1-c1^2, s2^2->1-c2^2, s3^2->1-c3^2} ]

Out[3]= {{1, 0, 0}, {0, 1, 0}, {0, 0, 1}}

There are other case, however, where the simplification rules are
not so simple, in which I have no idea how to proceed.

Mark Phillips
mbp at

Mathematics and Computer Science
University of Richmond			
Richmond, Virginia  23173

phone: (804) 289-8090
fax: (804) 289-8943

  • Prev by Date: Hashing for efficient definition
  • Next by Date: Mathematica Conference?
  • Previous by thread: Question about Simplify and related
  • Next by thread: Re: Question about Simplify and related