a=matrix([[0,0,1,-3,2],[2,-1,4,2,1],[4,-2,9,1,4],[2,-1,5,-1,5]])
a.rank()
a.row_space()
a.column_space()
a.left_nullity()
a*x=0
a.det()