solve set of linear equations with syms in matlab

main.m
syms x y z
eqns = [2*x + y - z == 3;
        x - y + 3*z == -1;
        3*x + 2*y + z == 4];
sol = solve(eqns, [x y z]); % solve the equations
disp(sol) % print the solutions
175 chars
7 lines

main.m
x = 7/5
y = 1/5
z = -2/5
25 chars
4 lines

main.m
syms x y z
A = [2 1 -1; 1 -1 3; 3 2 1]; % coefficient matrix
B = [3; -1; 4]; % constant vector
X = A \ B; % solve using inverse matrix
disp(X) % print the solutions
165 chars
6 lines

main.m
X =

    7/5
    1/5
   -2/5
29 chars
6 lines

main.m
syms x y z
A = [2 1 -1; 1 -1 3; 3 2 1]; % coefficient matrix
B = [3; -1; 4]; % constant vector
AB = [A B]; % augmented matrix
r = rank(A); % rank of A
rref_AB = rref(AB); % row echelon form of AB
X = zeros(3, 1); % initialize solution vector
X(r) = rref_AB(1:r, end); % solve for variables using back substitution
for i = r-1:-1:1
    X(i) = rref_AB(i, end) - rref_AB(i, i+1:end) * X(i+1:end);
end
disp(X) % print the solutions
428 chars
13 lines

main.m
X =

    7/5
    1/5
   -2/5
29 chars
6 lines