"Testing eigenvec"

clear

A = ((1,2,3),(2,6,4),(3,4,5))
Q = eigenvec(A)
D = dot(transpose(Q),A,Q)
B = dot(Q,D,transpose(Q))
check(infixform(B) == "((1,2,3),(2,6,4),(3,4,5))")

A = ((1,0,0,0),(0,-1,2,0),(0,2,-1,0),(0,0,0,1))
Q = eigenvec(A)
check(infixform(Q) == "((1,0,0,0),(0,0.707107,-0.707107,0),(0,0.707107,0.707107,0),(0,0,0,1))")
D = dot(transpose(Q),A,Q)
B = dot(Q,D,transpose(Q))
check(infixform(B) == "((1,0,0,0),(0,-1,2,0),(0,2,-1,0),(0,0,0,1))")
