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