a=matrix([[1,-1],[2,4]])
b=a.eigenvalues()
c=b[0]+b[1]
d=b[0]*b[1]
e=bool(c == a.trace())
f=bool( d == a.det())
e and f