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

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

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

/* code_ends

Strang-1.5-9-U (last edited 2010-12-17 10:43:30 by 172)