(a=matrix([[4,0,4],[0,2,0],[6,0,1]])
print a
b=inv(a)
print b