a = var('a')
b=var('b')
c=var('c')
d=var('d')
A = matrix([[a,b],[c,d]])
print A 
print "let A be invertible, then inverse of A^T would be :"
print ""
B = A.transpose()
C = B.inverse()
print C
print "===================="
print " given if A is also symmetric, i.e. A=A^T"
A = matrix([[a,b],[b,d]])
print A
D= A.inverse()
E= D.transpose()
print E
print "===================="
A = matrix([[2,1],[1,1]])
B = A.transpose()
C = B.inverse()
print C
print ""
D= A.inverse()
E= D.transpose()
print E

Strang-1.7-1.22-U (last edited 2010-12-17 19:27:28 by ch.sandeep)