% GEEXAMPLE  bigger, random example of Gaussian elimination; see GEBUELER
% Examples:
%   >> geexample
%   >> geexample, while norm(ans)<1e-11, geexample, end, norm(ans)

A = randn(8,8);
x = randn(8,1);
b = A * x;               % now A x = b is random system with known solution

Ab = [A b]               % augmented matrix for this system

Uc = gebueler(Ab)        % do Gaussian elimination

X = Uc(:,1:8) \ Uc(:,9)  % does back-substitution; X should be identical to x

x - X                    % should be zero, but we would accept approx 10^-15 ...
