def transform_to_identity(a):
a = ge.transform_to_row_echleon(a)
n = a.shape[0]
for i in xrange(1, n):
print i
for j in xrange(0, i):
factor = a[j][i] / a[i][i]
for k in xrange(0, n + 1):
a[j][k] -= factor * a[i][k]
a = ge.row_scaling(a, False)
return a
def transform_to_identity(a): a = ge.transform_to_row_echleon(a) n = a.shape[0] for i in xrange(1,n): print i for j in xrange(0,i): factor = a[j][i]/a[i][i] for k in xrange(0,n+1): a[j][k] -= factor*a[i][k] a = ge.row_scaling(a, False) return a
def solve(a, b):
print "Coefficients:\n", a
print "Values:\n", b
p = mo.augmented_matrix(a, b)
print 'Augmented:\n', p
p = ge.row_scaling(p)
print 'Row scaled:\n', p
# p = partial_pivoting(p)
# print "Partially pivoted:\n", p
p = transform_to_identity(p)
print 'Transformed to identity form:\n', p
q = ge.backward_substitution(p)
print "On backward substitution:\n", q
def solve(a,b): print "Coefficients:\n", a print "Values:\n", b p = mo.augmented_matrix(a,b) print 'Augmented:\n', p p = ge.row_scaling(p) print 'Row scaled:\n', p # p = partial_pivoting(p) # print "Partially pivoted:\n", p p = transform_to_identity(p) print 'Transformed to identity form:\n', p q = ge.backward_substitution(p) print "On backward substitution:\n", q
