L=matrix([[1,0,0],[1,1,0],[1,1,1]])
U=matrix([[1,1,1],[0,1,1,],[0,0,1]])
b=matrix([[4],[5],[6]])
print '---L---'
print L
print '----U----'
print U
print '----b----'
print b
print '-----L.inverse----'
x1=L.inverse()
print '--------'
print x1
c= x1*b
print '---c----'
print c
x2=U.inverse()
print '---U.inverse---'
print x2
x3=x2*c
print '---x----'
print x3
A=L*U
print '---A---'print A

Strang-1.5-33-U (last edited 2010-12-18 06:13:23 by jejiramchand)