% ENM 3.2
% Denna kod skrevs på övningen 2017-09-07 i SF1545

clear

A = [4 -1 -1 0; -1 4 0 -1; -1 0 4 -1; 0 -1 -1 4]
b = [200; 400; 0; 200]

% Fyra olika sätt att lösa ekvationssystemet
% (generellt rekommenderas \ eller linsolve,
% de är snabbast och mest stabila)
u1 = A\b
u2 = inv(A)*b
u3 = linsolve(A,b)
R = rref([A,b]), u4 = R(:,end)

% Räkna ut residualen (borde bli nära noll)
r1 = A*u1 - b
r2 = A*u2 - b
r3 = A*u3 - b
r4 = A*u4 - b
