l=matrix([[1,0,0,0],[2,1,0,0],[2,1,1,0],[3,2,4,1]])
u=matrix([[1,2,0,1,2,1],[0,0,2,2,0,0],[0,0,0,0,0,1],[0,0,0,0,0,0]])
A=l*u
A.rank()
A.left_kernel()

Strang-2.8.(a)-U (last edited 2010-12-18 05:57:30 by mamatha)