var('a,b,c,d,e,f')
A=matrix([[a,b],[b,c]])
B=matrix([[d,e],[e,f]])
w=(A+B)+(A-B)
if w[0,1]==w[1,0]:
    print('Symmetric')
else:
    print ('not Symmetric')

Strang-1.6-54(b)-U (last edited 2010-12-17 11:44:02 by xxx)