#To find the incidence matrix and its null-space
A  = matrix([[1,0,1],[-1,1,0],[0,-1,-1]])
print A
print 'The determinant of the matrix is '
print det(A)
print 'This means infinitely many solutions'
print 'One of which is'
print vector([2,2,-2])
print 'Its null-space consists of linearly-independent elements of the form'
print '(x,x,-x)'
b = A.transpose()
print 'The transpose is '
print b
print 'Its null-space consists of linearly-independent elements of the form'
print '(x,x,x)'

* Solution by:

Strang-2.5-2-U (last edited 2010-12-18 10:44:28 by ArnabBasu)