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

Strang-1.6-54(a)-U (last edited 2010-12-17 11:41:56 by xxx)