reset()
A=matrix([[3,4,2],[0,1,2],[0,0,0]])
show(A)
lambda1,lambda2,lambda3 = A.eigenvalues()
print "eigen values of matrix A ",(lambda1,lambda2,lambda3)
print "lambda1 + labda2 + lambda3 = ",lambda1 + lambda2 + lambda3
print "lambda1 * labda2 * lambda3 = ",lambda1 * lambda2 * lambda3
print "trace of matrix A =",A.trace()
lambda1*lambda2*lambda3
print "determinant of matrix A",A.determinant()
print '-------------------------------------'
B=matrix([[0,0,2],[0,2,0],[2,0,0]])
show(B)
lambda1,lambda2,lambda3 = B.eigenvalues()
print "eigen values of matrix A ",(lambda1,lambda2,lambda3)
print "lambda1 + labda2 + lambda3 = ",lambda1 + lambda2 + lambda3
lambda1,lambda2,lambda3
lambda1 + lambda2 + lambda3
print "lambda1 * labda2 * lambda3 = ",lambda1 * lambda2 * lambda3
print "trace of matrix B =",B.trace()
print "determinant of matrix B=",B.determinant()

Strang-5.1-7-U (last edited 2010-12-17 09:15:48 by tonylijo)