A = matrix([[2, -1], [-1, 2]])
B=copy(A)
Asqr= A*A
Ainv= B.inverse()
x1=A.eigenvalues()
x2=Asqr.eigenvalues()
x3=Ainv.eigenvalues()
print 'eigen values of A are ', x1
print 'eigen values of A*A are ', x2
print 'eigen values of A inverse are ', x3
print 'eigen vectors of A',A.eigenvectors_right()
print 'eigen vectors of A*A are ',Asqr.eigenvectors_right()
print 'eigen vectors of Ainverse are ',Ainv.eigenvectors_right()

Strang-6.1-6.1A (last edited 2010-12-18 08:03:36 by vasanthi emani)