B=matrix([[1,3],[1,1]])
C=B.transpose()
A=B+C
print A
D=A.transpose()
print D
print 'A is symmetric'
K=B-C
print K
E=K.transpose()
print E
print 'K is skew-symmetric'
bool(B==A/2 +K/2)

Strang-1.6-14 (last edited 2010-12-18 11:26:46 by sai kumar)