function [x] = elimgauss(A,B,n)
%       Program ini bertujuan menyelesaikan Persamaan Aljabar Linier
%       dengan metode Eliminasi Gauss untuk sistem n x n
%       n adalah orde matriks
%       A adalah matriks n x n
%       B adalah matriks n x 1
%       x adalah output matriks n x 1 untuk penyelesaian A x = B

%       Nama File : elimgauss.m
%       Surakarta, Oktober 2005
%   -------------------------------------------------------------------

%   Eliminasi Gauss
for k=1:n
    for i=k+1:n
        m(i,k)=A(i,k)/A(k,k);
        for j=k:n
            A(i,j)=A(i,j)-m(i,k)*A(k,j);
        end
        B(i)=B(i)-m(i,k)*B(k);
    end
end

%   Substitusi Balik
x(n)=B(n)/A(n,n);
for i=n-1:-1:1
    sum=B(i);
    for k=i+1:n
        sum=sum-A(i,k)*x(k);
    end
    x(i)=sum/A(i,i);
end