format rational A=[2 -1 4 -5 19; -1 3 -2 -1 -18; 0 1 5 -6 12; 3 -8 1 0 29] A = 2 -1 4 -5 19 -1 3 -2 -1 -18 0 1 5 -6 12 3 -8 1 0 29 Aog=A Aog = 2 -1 4 -5 19 -1 3 -2 -1 -18 0 1 5 -6 12 3 -8 1 0 29 A(1,:)=A(1,:)+A(2,:) A = 1 2 2 -6 1 -1 3 -2 -1 -18 0 1 5 -6 12 3 -8 1 0 29 A(2,:)=A(1,:)+A(1,:) A = 1 2 2 -6 1 2 4 4 -12 2 0 1 5 -6 12 3 -8 1 0 29 A(2,:)=A(1,:)-A(1,:) A = 1 2 2 -6 1 0 0 0 0 0 0 1 5 -6 12 3 -8 1 0 29 A(2,:)=A(1,:)-A(1,:) A = 1 2 2 -6 1 0 0 0 0 0 0 1 5 -6 12 3 -8 1 0 29 A=Aog A = 2 -1 4 -5 19 -1 3 -2 -1 -18 0 1 5 -6 12 3 -8 1 0 29 A(1,:)=A(1,:)+A(2,:) A = 1 2 2 -6 1 -1 3 -2 -1 -18 0 1 5 -6 12 3 -8 1 0 29 A(2,:)=A(2,:)+A(1,:) A = 1 2 2 -6 1 0 5 0 -7 -17 0 1 5 -6 12 3 -8 1 0 29 A(4,:)=A(4,:)-3*A(1,:) A = 1 2 2 -6 1 0 5 0 -7 -17 0 1 5 -6 12 0 -14 -5 18 26 A(2,:)=A(2,:)-4*A(3,:) A = 1 2 2 -6 1 0 1 -20 17 -65 0 1 5 -6 12 0 -14 -5 18 26 A(3,:)=A(2,:)-A(3,:) A = 1 2 2 -6 1 0 1 -20 17 -65 0 0 -25 23 -77 0 -14 -5 18 26 A(4,:)=A(2,:)1*4*A(2,:) A(4,:)=A(2,:)1*4*A(2,:) ↑ {Error: Invalid expression. Check for missing multiplication operator, missing or unbalanced delimiters, or other syntax error. To construct matrices, use brackets instead of parentheses. } A(4,:)=A(2,:)+14*A(2,:) A = 1 2 2 -6 1 0 1 -20 17 -65 0 0 -25 23 -77 0 15 -300 255 -975 A(4,:)=A(2,:)-14*A(2,:) A = 1 2 2 -6 1 0 1 -20 17 -65 0 0 -25 23 -77 0 -13 260 -221 845 A = A = ↑ {Error: Invalid expression. Check for missing or extra characters. } A = Aog A = 2 -1 4 -5 19 -1 3 -2 -1 -18 0 1 5 -6 12 3 -8 1 0 29 A(1,:)=A(1,:)+A(2,:) A = 1 2 2 -6 1 -1 3 -2 -1 -18 0 1 5 -6 12 3 -8 1 0 29 A(2,:)=A(2,:)+A(1,:) A = 1 2 2 -6 1 0 5 0 -7 -17 0 1 5 -6 12 3 -8 1 0 29 A(4,:)=A(4,:)-3*A(1,:) A = 1 2 2 -6 1 0 5 0 -7 -17 0 1 5 -6 12 0 -14 -5 18 26 Aog2=a {Undefined function or variable 'a'. } Aog2=A Aog2 = 1 2 2 -6 1 0 5 0 -7 -17 0 1 5 -6 12 0 -14 -5 18 26 A(2,:)=A(2,:)-4*A(3,:) A = 1 2 2 -6 1 0 1 -20 17 -65 0 1 5 -6 12 0 -14 -5 18 26 A(3,:)=A(2,:)-A(3,:) A = 1 2 2 -6 1 0 1 -20 17 -65 0 0 -25 23 -77 0 -14 -5 18 26 A(4,:)=A(4,:)+14*A(2,:) A = 1 2 2 -6 1 0 1 -20 17 -65 0 0 -25 23 -77 0 0 -285 256 -884 Aog3=A Aog3 = 1 2 2 -6 1 0 1 -20 17 -65 0 0 -25 23 -77 0 0 -285 256 -884 A(3,:)=A(3,:)/-25 A = 1 2 2 -6 1 0 1 -20 17 -65 0 0 1 -23/25 77/25 0 0 -285 256 -884 A(4,:)=A(4,:)+285*A(3,:) A = 1 2 2 -6 1 0 1 -20 17 -65 0 0 1 -23/25 77/25 0 0 0 -31/5 -31/5 A(4,:)=A(4,:)*5 A = 1 2 2 -6 1 0 1 -20 17 -65 0 0 1 -23/25 77/25 0 0 0 -31 -31 A(4,:)=A(4,:)/-31 A = 1 2 2 -6 1 0 1 -20 17 -65 0 0 1 -23/25 77/25 0 0 0 1 1 Aref=A Aref = 1 2 2 -6 1 0 1 -20 17 -65 0 0 1 -23/25 77/25 0 0 0 1 1 A(1,:)=A(1,:)-2*A(2,:) A = 1 0 42 -40 131 0 1 -20 17 -65 0 0 1 -23/25 77/25 0 0 0 1 1 A(1,:)=A(1,:)-42*A(3,:) A = 1 0 0 -34/25 41/25 0 1 -20 17 -65 0 0 1 -23/25 77/25 0 0 0 1 1 A(2,:)=A(2,:)-20*A(3,:) A = 1 0 0 -34/25 41/25 0 1 -40 177/5 -633/5 0 0 1 -23/25 77/25 0 0 0 1 1 A(2,:)=A(2,:)+20*A(3,:) A = 1 0 0 -34/25 41/25 0 1 -20 17 -65 0 0 1 -23/25 77/25 0 0 0 1 1 A(2,:)=A(2,:)+20*A(3,:) A = 1 0 0 -34/25 41/25 0 1 0 -7/5 -17/5 0 0 1 -23/25 77/25 0 0 0 1 1 A(1,:)=A(1,:)+(34/25)*A(5,:) {Index in position 1 exceeds array bounds (must not exceed 4). } A(1,:)=A(1,:)+(34/25)*A(4,:) A = 1 0 0 * 3 0 1 0 -7/5 -17/5 0 0 1 -23/25 77/25 0 0 0 1 1 A(2,:)=A(2,:)+(7/5)*A(4,:) A = 1 0 0 * 3 0 1 0 * -2 0 0 1 -23/25 77/25 0 0 0 1 1 A(3,:)=A(3,:)+(23/5)*A(4,:) A = 1 0 0 * 3 0 1 0 * -2 0 0 1 92/25 192/25 0 0 0 1 1 A(3,:)=A(3,:)-(23/5)*A(4,:) A = 1 0 0 * 3 0 1 0 * -2 0 0 1 -23/25 77/25 0 0 0 1 1 A(3,:)=A(3,:)-(23/5)*A(4,:) A = 1 0 0 * 3 0 1 0 * -2 0 0 1 -138/25 -38/25 0 0 0 1 1 A(3,:)=A(3,:)+(23/5)*A(4,:) A = 1 0 0 * 3 0 1 0 * -2 0 0 1 -23/25 77/25 0 0 0 1 1 A(3,:)=A(3,:)-(23/5)*A(4,:) A = 1 0 0 * 3 0 1 0 * -2 0 0 1 -138/25 -38/25 0 0 0 1 1 A(3,:)=A(3,:)+(23/5)*A(4,:) A = 1 0 0 * 3 0 1 0 * -2 0 0 1 -23/25 77/25 0 0 0 1 1 A(3,:)=A(3,:)+(23/25)*A(4,:) A = 1 0 0 * 3 0 1 0 * -2 0 0 1 * 4 0 0 0 1 1 Aans=A Aans = 1 0 0 * 3 0 1 0 * -2 0 0 1 * 4 0 0 0 1 1 diary off rref(Aog) ans = 1 0 0 0 3 0 1 0 0 -2 0 0 1 0 4 0 0 0 1 1 diary off