var('a,b,c')
E=matrix([[1,0,0],[a,1,0],[b,0,1]])
F=matrix([[1,0,0],[0,1,0],[1,c,1]])
EF=E*F
FE=F*E
EE=E*E
print "EF:\n\n",EF
print "FE:\n\n",FE
print "EE:\n\n",EE

Strang-1.4.34-U (last edited 2010-12-18 11:42:22 by B lavanya)