Примеры решение задач линейной алгебры

ЗАДАЧА. Решить матричное уравнение А • Х= В.

with(linalg);
A:=matrix(2,2,[1,2,3,4]);
A := matrix([[1, 2], [3, 4]])
M:=matrix(2,2,[5,6,7,8]);
M := matrix([[5, 6], [7, 8]])
evalm(inverse(A)&*M);#первый  способ
matrix([[-3, -4], [4, 5]])
linsolve(A,M); #второй способ
matrix([[-3, -4], [4, 5]])

ЗАДАЧА. Решить систему линейных уравнений различными методами:
2*х1- х2+ х3= 0,
3*x1+ 2*x2-5*x3=1,
х1+ 3*х2- 2*x3= 4.

with(linalg);
A:=matrix(3,3,[2,-1,1,3,2,-5,1,3,-2]):b:=vector(3,[0,1,4]) ;
b := vector([0, 1, 4])
linsolve(A,b);#решение линейной  системы Ax*b
vector([13/28, 47/28, 3/4])
evalm(inverse(A)&*b);#решение Ax*b методом обратной матрицы
vector([13/28, 47/28, 3/4])
col((rref(augment(A,b)),4));#решение Ax=b методом Гаусса
vector([13/28, 47/28, 3/4])