A=matrix([[1,2],[2,4]])
B=matrix([[2,1,4],[3,1.5,6]])
var('c,b')
M=matrix([[1,b],[c,b*c]])
print '---A---'
print A
print '---B---'
print B
print '---M---'
print M
print 'Rank of AB'
print (A*B).rank()
print 'Rank of AM'
print (A*M).rank()

Strang-2.2-22-U (last edited 2010-12-18 07:55:13 by jejiramchand)