Chiến lược chọn trụ trong phương pháp khử Gauss Nguyễn Quốc Lân - BM Toán Ứng Dụng, DHBK TP. HCM > with(linalg): > Digits:=4; Khai báo ma trận AA: ma trận mở rộng của hệ Ax = b AA:=matrix(2,3,[0.003,59.14,59.17,5.291,-6.130,46.78]); Khử cột 1: phần tử trụ khác 0, tính hệ số khử m21:= AA[2,1]/AA[1,1]; Biến đổi sơ cấp trên hàng AA:=addrow(AA,1,2,-m21); Ý nghóa của giải thuật khử Gauss: Phần tử được khử phải triệt tiêu AA[2,1]:=0; > evalm(AA); Giải lùi x:=backsub(AA); > B:=matrix(2,3,[0.003,59.14,59.17,5.291,-6.130,46.78]); Hoán vò hàng để tìm ra trụ tối đại BB := swaprow(B,1,2); > m21:= BB[2,1]/BB[1,1]; > BB:=addrow(BB,1,2,-m21); > x:=backsub(BB); Trường hợp thất bại của phương pháp trụ tối đại C:=matrix(2,3,[30.00,591400.,591700.,5.291,-6.130,46.78]); > m21:=C[2,1]/C[1,1]; > C:=addrow(C,1,2,-m21); > C[2,1]:=0; > evalm(C); > x:=backsub(C); >