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

Strang-1.6-54(c)-U (last edited 2010-12-17 11:46:51 by xxx)