Recover the history of the row operations that have been performed. This function combines the transformation matricies into a single transformation matrix representing all row operations or may optionally print all the individual operations which have been performed.
1 2 3 4 5 6 7 
x 
a matrix A, joined with a vector of constants, b, that has been passed to

all 
logical; print individual tranformation matricies? 
... 
additional arguments 
the tranformation matrix or a list of individual transformation matricies
echelon
, gaussianElimination
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26  A < matrix(c(2, 1, 1,
3, 1, 2,
2, 1, 2), 3, 3, byrow=TRUE)
b < c(8, 11, 3)
# using row operations to reduce below diagonal to 0
Abt < Ab < cbind(A, b)
Abt < rowadd(Abt, 1, 2, 3/2)
Abt < rowadd(Abt, 1, 3, 1)
Abt < rowadd(Abt, 2, 3, 4)
Abt
# build T matrix and multiply by original form
(T < buildTmat(Abt))
T %*% Ab # same as Abt
# print all transformation matricies
buildTmat(Abt, TRUE)
# invert transformation matrix to reverse operations
inv(T) %*% Abt
# gaussian elimination
(soln < gaussianElimination(A, b))
T < buildTmat(soln)
inv(T) %*% soln

Questions? Problems? Suggestions? Tweet to @rdrrHQ or email at ian@mutexlabs.com.
All documentation is copyright its authors; we didn't write any of that.