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-4 (last edited 2010-12-17 09:13:09 by 172)