L=matrix([[1,0,0],[4,1,0],[1,0,1]])
print L
U=matrix([[2,2,4],[0,1,3],[0,0,1]])
print U
A=L*U
print A
B=vector([0,0,1])
print B
x=A.solve_right(B)
print x
c=L.solve_right(B)
print c
x=U.solve_right(c)
print x

Strang-1-14-U (last edited 2010-12-18 09:07:12 by prakasam team)