B = matrix ([[0,9,-12],[-9,0,20],[12,-20,0]])
print B

print B.transpose()
print "Is symmetric?",B.is_symmetric()
print "Is skew_symmetric?", B.is_skew_symmetric()

if ( B.is_invertible()==true):
    if(B.transpose()==B.inverse()):
        print "matrix is orthogonal"
    else:
        print "matrix is not orthogonal"
else:
    print"matrix is not orhtogonal"
B.eigenvalues()

A = matrix ([[3,4],[1,3]])
print A

print A.transpose()
print "Is symmetric?",A.is_symmetric()
print "Is skew_symmetric?", A.is_skew_symmetric()

if ( A.is_invertible()==true):
    if(A.transpose()==A.inverse()):
        print "matrix is orthogonal"
    else:
        print "matrix is not orthogonal"
else:
    print"matrix is not orhtogonal"



A.eigenvalues()

Kreyszig-7.3-3 (last edited 2010-12-15 13:14:15 by 172)