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