M=matrix([[4,4,2],[3,-1,2],[3,7,1]])
B=matrix([[0],[0],[0]])
print B
[0]
[0]
[0
Print M
[ 4  4  2]
[ 3 -1  2]
[ 3  7  1]
A=M.augment(B)
E=A.echelon_form()
print E
[ 1  5  0  0]
[ 0  8 -1  0]
[ 0  0  0  0]
//notice that set of linear equations has infinite solutions

solutions x+5*y=0
          8*y-z=0

Kreyszig-18.1.10-U (last edited 2010-12-18 11:23:01 by Aditya J)