reset()
A = matrix([[1,1,1,1],[1,1,1,1],[1,1,1,1],[1,1,1,1]])
print "matrix A:"
show(A)
print "matrix B = A - I"
B = A - identity_matrix(4)
show(B)
print "eigen values of B",B.eigenvalues()
print "Determinet of B",B.determinant()

Strang-5.1-18-U (last edited 2010-12-17 09:25:22 by tonylijo)