l=matrix([[1,0],[4,1]])
l
b=matrix([[2],[11]])
b
c=l.inverse()*b
c
U=matrix([[2,4],[0,1]])
U
x=U.inverse()*c
x