a=matrix([[1,1/2,1/3],[1/2,1/3,1/4],[1/3,1/4,1/5]])#coefficient matrix
c=a^-1 #inverse if matrix a
b=matrix(QQ,[[1,0,0]]) #defining the constants matrix
d=c*b.transpose()
print d

SageDays/Kreyszig-20.3-9 (last edited 2010-08-12 18:52:19 by PraveenKumar)