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

S.K.Manoj, Student, IIIT Hyderabad