; A := matrix(3,4,[1,-2,-1,b1, -4,5,2,b2, -4,7,4,b3]); pivot(A,1,1); pivot(%,2,2); pivot(%,3,3); mulrow(%,2,-1/3); mulrow(%,3,3/2); rref(A);