Metode Numerik Eliminasi Gauss MATLAB Script

Dalam postingan kali ini, saya akan memberikan Script M-File untuk melakukan Eliminasi Gauss pada Sistem Persamaan Linear dalam perhitungan analisis numerik atau metode numerik. Berikut ini coding script yang dapat dijalankan pada program MATLAB dan script ini dapat dituliskan di Blank M-File MATLAB yang selanjutnya dapat dirun dengan nama GAUSS_ELIM:

function x = GAUSS_ELIM(A, b)

A=[2 0 0 0; 1 3 0 0; 3 4 2 0; -1 3 -6 -1];

b=[-10; 4; 2; 8];

n = size(A, 1);

r = zeros(n, 1);

for i = 1 : 1 : n

    r(i) = i;

end

x=zeros(n,1);

for k = 1 : 1 : n

    max = abs(A(r(k), r(k)));

    max_pos = k;   

    for l = k : 1 : n

        if abs(A(r(l), r(k))) > max;

            max=abs(A(r(l), r(k)));

            max_pos = l;               

        end

    end

    temp_r=r;

    r(k)=temp_r(max_pos);

    r(max_pos)=temp_r(k);

    for i = 1:1:n

        if i ~= k

            zeta = A(r(i),k)/A(r(k),k);

            for j = k : 1 : n

                A(r(i),j)=A(r(i),j)-A(r(k),j)*zeta;

            end

            b(r(i)) = b(r(i)) - b(r(k))*zeta;

        end

    end

end

for i = 1 : 1 : n

    x(i) = b(r(i))/A(r(i),i);

end

end