Chiến lược chọn trụ 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; om Khai báo ma trận AA: ma trận mở rộng hệ Ax = b C AA:=matrix(2,3,[0.003,59.14,59.17,5.291,-6.130,46.78]); m21:= AA[2,1]/AA[1,1]; AA:=addrow(AA,1,2,-m21); Zo Biến đổi sơ cấp hàng ne Khử cột 1: phần tử trụ khác 0, tính hệ số khử tiêu Si Giải lùi nh AA[2,1]:=0; > evalm(AA); Vi en Ý nghóa giải thuật khử Gauss: Phần tử khử phải triệt 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 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 phương pháp trụ tối đại C:=matrix(2,3,[30.00,591400.,591700.,5.291,-6.130,46.78]); SinhVienZone.com https://fb.com/sinhvienzonevn m21:=C[2,1]/C[1,1]; C:=addrow(C,1,2,-m21); C[2,1]:=0; evalm(C); x:=backsub(C); Si nh Vi en Zo ne C om > > > > > > SinhVienZone.com https://fb.com/sinhvienzonevn ...m21:=C [2, 1]/C[1,1]; C:=addrow(C,1 ,2, -m21); C [2, 1]:=0; evalm(C); x:=backsub(C); Si nh Vi en Zo ne C om > > > > > > SinhVienZone. com https://fb .com/ sinhvienzonevn