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

Strang-1.5-31-u-U (last edited 2010-12-16 12:39:37 by ravichandra)