% Instability of Gaussian
% elimination w/o pivoting
A = [2^(-55), 1; 1 , 1];
b = [1;2];
x1 = A\b;
x2 = GEpivot(A,b);
[L,U] = lufak(A);
z = L\b; x3 = U\z;
format short; disp([x1,x2,x3]);