# Example 3, Page 238
lambda1 = var('lambda1')
A = matrix([[1 - lambda1, 4, 5], [0, 3/4 - lambda1, 6], [0, 0, 1/2 - lambda1]])
print A.eigenvalues()

Strang-5.1-3 (last edited 2010-12-17 03:53:38 by VenkataRamana)