A=matrix(QQ,([2,3,0],[4,5,-1],[2,-1,-3]))
print A
b=vector([0,3,5])
c=A.augment(b.transpose())
print
print c.echelon_form()