% Example for THEORETICAL Gaussian elimimination
% It is sheer luck that roundoff errors do not affect these computations

% A linear system with a 3x5 matrix
A = [0 1 2 3 1;0 4 5 6 1;0 7 8 9 1];
b = [5;11;17];

% Transformation to step normal form by row transformations
[Z,y,r,I] = GE_theoretical(A,[5;11;17]),
% Z contains the transformed matrix
% y contains the transformed right hand side vector
% r gives the rank of the matrix
% I contains the indices of columns of Z agreeing with unit vectors, see
% Def. 3.3.C



