<<VilGebra`

A=(({{8, 0, -16, -4}, {-3, -3, 3, 2}, {0, 3, 2, -2}, {1, 3, 3, 2}}));
B=Join[A,IdentityMatrix[4],2];
RowReduceComplete[B]