A=matrix([[5,4],[1,2]])
print A
A.eigenvectors_left()
X=matrix([[4,1],[1,-1]])
print X
print X^-1 * A * X