F=matrix(RDF,[[1,0,0,0],[2,1,0,0],[0,0,1,0],[0,0,0,1]])
print "F"
print F
G=matrix(RDF,[[1,0,0,0],[0,1,0,0],[0,2,1,0],[0,0,0,1]])
print "G"
print G
H = matrix(RDF,[[1,0,0,0],[0,1,0,0],[0,0,1,0],[0,0,2,1]])
print "H"
print H
print "F*G*H"
print F*G*H
print "H*G*F"
print H*G*F

Strang-1.5-4-U (last edited 2010-12-17 10:31:51 by 172)