b=matrix([[1,2,0],[0,1,1],[1,2,0]])
c=vector([1,0,1])
d=b.augment(c.transpose())
print d
q= d.echelon_form()
print q
q.rank()
b.pivots()                                     /*    ONE FREE VARIABLE  (R  PIVOTS AND N-R GIVE FREE VARIABLES    */
a=matrix([[1,2],[4,5],[7,8]])
b=vector([3,6,9])
c=a.augment(b.transpose())
print c
q= c.echelon_form()
print q
q.rank()
a.pivots()                               /*     ONE FREE VARIABLE      */

Strang-2.2-5-U (last edited 2010-12-18 07:40:12 by tejaswi)