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

Strang-1.5-31-U (last edited 2010-12-16 12:07:32 by lakshmi)