D=matrix([[1,0,0],[5,1,0],[4,0,-1]])
print D
D.det()
(D/2).det()
(-D).det()
(D**2).det()
E=D.inverse()
print E.det()