a=matrix([[1,2,3],[0,4,0],[0,0,5]])
cofr=a.adjoint()
print cofr
cofr=cofr.transpose()
adet=a.det()
c=identity_matrix(3)*adet
x=cofr.transpose()
print "Cofactor matrix"
print x
d=x*a
if c==d:
 print "acT=(det a)I"
else:
 print "acT!=(det a)I"
invA=a.inverse()
print "Inverse matrix"
print invA

Kreyszig-4.4.1-U (last edited 2010-12-18 05:38:28 by B lavanya)