A=([2,3],[10,9])
l=A[1][0]/A[0][0]
print 'coefficient = ',l
B=([A[0][0]*l,A[0][1]*l],[10-A[0][0]*l,9-A[0][1]*l])
print 'Upper triangular system : ' , B
print 'The pivots are : ' , B[0][0] , 'and' , B[1][1]

Strang-1.3-4-U (last edited 2010-12-18 10:26:55 by Shubham)