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