; A := matrix(3,3,[2,1,2,2,2,-2,3,1,1]); I3 := diag(1,1,1); rref( augment(A-I3, matrix(3,1,[0,0,0])) ); rref( augment( A - 4*I3, matrix(3,1,[0,0,0]) ) );