g=matrix([[1,0,0],[0,1,0],[0,1,1]])
f=matrix([[1,0,0],[0,1,0],[1,0,1]])
e=matrix([[1,0,0],[-2,1,0],[0,0,1]])
p=~g
q=~f
r=~e
d=g*f*e
#print d
l=r*q*p
#print l
m=d*l
n=l*d
print m
print n 

# L is the upper triangular matrix and multiplication of (g*f*e) is also a upper triangular matrix # Multiplication of d=(g*f*e) and l=(r*q*p) l*d=d*l they both are same which is identity matrix

Strang-1.5-1-U (last edited 2010-12-17 04:28:19 by sree123)