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