a=matrix([[1,0,0,1,0,0],[-1,1,0,0,1,0],[0,-1,1,0,0,1]])
print '[A I]:',a
b=a.echelon_form()
print '[I A^-1]:',b

S.K.Manoj, Student, IIIT Hyderabad