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(d)-U (last edited 2010-12-17 11:50:34 by xxx)