A = matrix([[1,2],[3,8]])
print A
U = matrix([[1,2],[0,2]])
print U
L = matrix([[1,0],[3,1]])
L*U
print A == L*U