function [ A,b ] = BasicGE( A,b )
%function [ A,b ] = BasicGE( A,b )
%   This routine does basic Gaussian Elimination.
%   Inputs:
%       A - an nxn coeficient matrix.
%       b - a n-vector (column vector).
%   Outputs:

M=[A b];
[n,m]=size(M);

for i=1:(n-1)
    for j=(i+1):n
        c=M(j,i)/M(i,i);
        M(j,:)=M(j,:)-c*M(i,:);
    end
end

A=M(:,1:n);
b=M(:,(n+1):m);

end