A=matrix([[1,1,1],[-1,2,1],[2,1,-3]])
b=vector([4,8,-9])
C=b.transpose()
print "matrix A is:"
print A
P=A.det()
print "print det of A:",p
P=A.augment(C)
print P
P.swap_columns(3,0)
P1=P[:,:3]
print "matrix P1 is:"
print P1
p1=P1.det()
print "print det of P1:"
print p1
P=A.augment(C)
P.swap_columns(3,1)
P2=P[:,:3]
print "matrix P2 is:",P2
p2=P2.det()
print "print det of P2:"
print p2
P=A.augment(C)
P.swap_columns(3,2)
P3=P[:,:3]
print "matrix P3 "
print P3
p3=P3.det()
print "print det of P3:",p3
x=p1/p
y=p2/p
z=p3/p
print x,y,z

Strang-6.6-3 (last edited 2010-12-17 10:53:28 by anoopkumarreddy)