A=matrix([[1,1,3],[0,4,6],[1,5,8]])
B=matrix([[1,1,3],[0,4,6],[0,0,1]])
C=matrix([[1,1,3],[0,4,6],[1,5,9]])
A.add_multiple_of_row(2,0,-1)
A.add_multiple_of_row(2,1,-1)
print A
print A[0,0]*A[1,1]*A[2,2]
print det(B)
print det(C)
print det(A*B)
At=A.transpose()
print det(At*A)
ct=C.transpose()
print det(ct)

Strang-4.2 16-U (last edited 2010-12-18 08:23:25 by shwetaradha)