a=matrix([[1,2,0,2,1],[-1,-2,1,1,0],[1,2,-3,-7,-2]])
a.echelon_form()