A=matrix([[1,1,1,1],[1,1,1,1],[1,1,1,1],[1,1,1,1]])
C=matrix([[0,1,0,1],[1,0,1,0],[0,1,0,1],[1,0,1,0]])
print"A "
show(A)
print"C"
show(C)
print"The eigen values of A"
show(A.eigenvalues())
print"The eigen values of C"
show(C.eigenvalues())
print"The rank of A"
show(A.rank())
print("The rank of C")
show(C.rank())
print"Eigen vectors corresponds to A"
show(A.eigenvectors_right())
print"Eigen vectors corresponds to C"
show(C.eigenvectors_right())
print "Eigen vector corresponds to non-zero eigen value 4 in A."
show(A.eigenvectors_right()[0][1][0])
print "Eigen vector corresponds to non-zero eigen value 2 in C."
show(C.eigenvectors_right()[0][1][0])
print "Eigen vector corresponds to non-zero eigen value -2 in C."
show(C.eigenvectors_right()[1][1][0])

Strang-5.1-16-U (last edited 2010-12-17 11:09:25 by syamkrishnanck)