giải phương trình tuyến tính bằng khử Gauss trên Maple (phương pháp tính)
Bạn đang xem bản rút gọn của tài liệu. Xem và tải ngay bản đầy đủ của tài liệu tại đây (15.12 KB, 1 trang )
GIẢI HỆ PHƯƠNG TRÌNH TUYẾN TÍNH BẰNG KHỬ
GAUSS TRÊN MAPLE
> restart:
with(linalg):
> A :=
matrix(4,5,[1,-1,2,-1,-8,2,-2,3,-3,-20,1,1,1,0,-2,1,-1,4,3,4]);
> m21 := A[2,1]/A[1,1];
m31 := A[3,1]/A[1,1];
m41 := A[4,1]/A[1,1];
> A := addrow(A,1,2,-m21);
> A := addrow(A,1,3,-m31);
A := addrow(A,1,4,-m41);
> A := swaprow(A,2,3);
> m43 := A[4,3]/A[3,3];
A := addrow(A,3,4,-m43);
> x := backsub(A);
> B :=
matrix(4,5,[1,-1,2,-1,-8,2,-2,3,-3,-20,1,1,1,0,-2,1,-1,4,3,4]);
> B := gausselim(B);
>
