A = matrix([[0,9,-12],[-9,0,20],[12,-20,0]])
print A
print "Is symmetric?", A.is_symmetric() 
print "Is skew_symmetric?", A.is_skew_symmetric() 
print A.transpose()==A.inverse()
if(A.transpose()==A.inverse()):
    print this is orthogonal
else:
    print "this is not orthogonal"

Kreyszig-7.3-prob 5 (last edited 2010-12-17 05:20:41 by sai9876543210abcd)