reset()
A = matrix([[1,4],[2,3]])
show(A)
B = A+identity_matrix(2)
show(B)
print "eigen values of A",A.eigenvalues()

print "eigen values of B",B.eigenvalues()
print "eigen vectors of matrix A",A.eigenvectors_right()[0][1][0],A.eigenvectors_right()[1][1][0]
print "eigen vectors of matrix B",B.eigenvectors_right()[0][1][0],B.eigenvectors_right()[1][1][0]
print "so the eigen vectors of A and B are the same and their eigen value of A is one less than that in B"

Strang-5.1-19-U (last edited 2010-12-17 10:31:59 by tonylijo)