% Detta är en skriptfil (kommandofil)
% Linjär algebra: Ax=b
%  x + 2y + 3z = 1
% 3x +      2z = 2
%      5y -  z = 0

clear all
format rat
A = [1 2 3;3 0 2;0 5 -1] %koefficientmatrisen
b = [1;2;0]  % högerledet
R = rref([A b]) % Vi radreducerar utökade matrisen
disp('Dvs lösningen är')
x = R(1,4) % x, y och z står i högraste kolonnen i R
y = R(2,4)
z = R(3,4)
disp('Kontroll'),pause
disp(['Med dessa värden insatta i ekvationssystemets'...
    'vänsterled så får vi högerledet'])
disp([x+2*y+3*z;3*x+2*z;5*y-z])