a=matrix([[1,1,1],[1,1,3],[2,5,8]])
b=matrix([[1,1,1],[0,0,2],[0,3,6]])
l=a*~b
l
p=matrix([[1,0,0],[0,0,1],[0,1,0]])
u=matrix([[1,1,1],[0,3,6],[0,0,2]])
p*a
l*u

Strang-1.5-14-U (last edited 2010-12-17 10:23:54 by lakshmi)