n = 5

%% Generate a random matrix and RHS
A = rand(n)

%% Generate a random RHS
b = rand(n,1)

%% Use naive G.E. to solve
x = gauss_elim_scalar(A, b)
%x = gauss_elim_vector(A, b)
%x = gauss_elim_outer(A, b)

%% Check the error
A*x - b