A = matrix(RDF,[[1,-1,0,0], [-1,2,-1,0], [0,-1,2,-1],[0,0,-1,2]])
print A.parent()
P, L, U = A.LU()
print "P"
print P
print "L"
print L
print "U"
print U

Strang-1.5-5 (last edited 2010-12-17 10:12:48 by 172)