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

Strang-1.1-1-U (last edited 2010-12-17 09:23:55 by 172)