A=matrix([[1,2],[2,4]])
P=matrix([[1,3],[3,6]])
P^-1*A*P