a=matrix([[0,1,0],[0,0,1],[1,0,0]])
a.transpose()==a.inverse()  /*returns boolean true checkn if a is orthogonal*/
var('t')
t=a*a*a     /*cube of a*/
t
t.inverse()==t.transpose()  /* returns boolean true if a^3 is orthogonal*/

Strang-3.5-3.7-U (last edited 2010-12-18 08:19:34 by manishadeepthi)