Transforms the given matrix to reduced row echelon form.

rrechelon.py

#!/usr/bin/env python3

import sys
from fractions import Fraction

mult = lambda r,f: [f * x for x in r]
add = lambda r1,r2: [x + y for x,y in zip(r1,r2)]

def print_matrix(m):
    for r in m: print('\t'.join(str(x) for x in r))
    print()

m = []
for l in sys.stdin:
    if l != '\n': m.append(list(map(Fraction.from_float,map(float,l.split()))))

print_matrix(m)
for i in range(len(m) - 1):
    for j in range(i+1, len(m)):
        m[j] = add(m[j], mult(m[i], - m[j][i] / m[i][i]))
    print_matrix(m)

for i in reversed(range(1, len(m))):
    for j in range(i):
        m[j] = add(m[j], mult(m[i], - m[j][i] / m[i][i]))
    print_matrix(m)

for i in range(len(m)):
    m[i] = mult(m[i], 1 / m[i][i])
print_matrix(m)