a=matrix([[1,2],[2,4]])
b=matrix([[2,1,4],[3,1.5,6]])
c=a*b
c.rank()
var('bc,c')
a=matrix([[1,2],[2,4]])
m=matrix([[1,b],[c,bc]])
t=a*m
t.rank()