A = matrix([[2,1,0],[1,2,1],[0,1,2]])
L = matrix([[1,0,0],[1/2,1,0],[0,2/3,1]])
U = matrix([[2,1,0],[0,3/2,1],[0,0,4/3]])
print "L*U is",L*U

Strang-2.6-1 (last edited 2010-12-18 10:42:44 by Siva Jyoshna Poreddy)