a=matrix([[1,2],[3,6]])
b=matrix([[0,0],[1,2]])
c=matrix([[1,1,0,0],[0,1,0,1]])
a.bases()
b.bases()
c.bases()